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Abstract 

Two families of cuprates exhibiting high-Tc superconductivity, YBa2Cu306+j; (YBCO) and 
La23;Sr3;Cu04+5 (LSCO), have been extensively studied by anelastic spectroscopy, by measuring 
the complex dynamic Young's modulus E {lo,T) = E' + iE" of ceramic samples at frequencies 
of 0.5-20 kHz between 1 and 900 K. Results are also presented on oxygen vacancies in the 
ruthenocupratc RuSr2GdCu208-5. The elastic energy loss curves as a function of temperature, 
(a;, T) = E"/E', contain peaks at the temperatures Tm such that UTiTm) — 1, where a;/27r 
is the measuring frequency and r is the relaxation time of any microscopic process coupled to 
strain, like hopping of O atoms, tilting of O octahedra or fluctuations of the hole stripes. By 
measuring such anelastic spectra at different frequencies, it is therefore possible to selectively 
probe the dynamics of the various relaxation processes, precisely determining their characteristic 
times t(T). The reliability of the anelastic experiments and of the assignments to various 
microscopic mechanisms is also discussed. 

The richest anelastic spectra are those of LSCO, and have been studied in the whole range 
of Sr doping (0 < x < 0.2) and at few Ba doping levels. As in the parent perovskite compounds, 
LSCO is made of O octahedra unstable against tilting, which give rise to low symmetry phases. 
The layered coordination of such octahedra and the possibility of achieving a low density of 
pinning defects (interstitial O atoms), makes it possible to observe solitonic tilt waves and fast 
local motion of the octahedra among the several minima of the local tilt potential far below 
the structural transition temperature. The fast local motion is driven by tunneling of the O 
atoms and is enormously enhanced and accelerated by even small doping, demonstrating direct 
coupling between the tilt modes of the octahedra and the hole excitations. 

In addition, it has been possible to probe the dynamics of the stripes into which the holes 
segregate, at liquid He temperature, when they act as walls between domains of antiferromag- 
netically correlated spins (cluster spin glass), and also at higher temperature, when they can 
overcome by thermal activation the pinning barriers provided by Sr^+ dopants. A picture is 
proposed in which the high temperature motion of the stripes involves the formation of kink 
pairs, while at lower temperature only the motion of the existing kinks can occur. 

In YBCO the anelastic spectra are dominated by the motion of the O atoms in the CuOj, 
planes, which are responsible for the doping of the charge carriers (holes) into the supercon- 
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ducting Cu02 planes. The doping level, and therefore the superconducting properties, are 
determined not only by the content x of highly mobile non-stoichiometric O, but also by its 
ordering. The various processes involved in ordering and diffusion of O have different character- 
istic times r, and therefore produce distinct peaks in the anelastic spectra at the temperatures 
for which ujT(Tm) = 1; they have been studied in the whole stoichiometry (0 < x < 1) and 
temperature (50 K < T < 800 K) ranges. It is shown that there are three types of O jumps with 
different rates, depending whether they involve: i) ordered Cu-0 chains in the orthorhombic 
O-I phase {x ~ 1); ii) sparser chains fragments in the O-II and tetragonal phases; Hi) isolated 
O atoms. The first two types of jumps occur over a barrier of ~1.0 eV, whereas the latter has 
a barrier of only 0.11 eV. It is discussed how such widely different barriers for O hopping are 
possible and how the extraordinarily high mobility of isolated O atoms is compatible with the 
slow times for O ordering. In addition to the diffusive jumps, hopping of O between off-center 
positions within the Cu-0 chains is proposed to occur. Finally, an anelastic process has been 
observed, whose intensity increases steeply in the overdoped state, x > 0.9, where all the other 
physical properties remain practically constant. Such a process is attributed to the reorienta- 
tion of small bipolarons on orbitals that do not contribute to the electrical conduction, and can 
be useful for characterizing materials with non-optimal O content, like thin films. 

RuSr2GdCu208-5 is isostructural with YBCO, except for Ru02-5 planes with 5 below few 
percents instead of the widely nonstoichiometric CuOa; planes. In this case, the situation of O 
mobility and ordering is simpler than in YBCO, and it is shown that it can be quantitatively 
explained in terms of hopping of O vacancies whose clastic quadrupoles are weakly interacting 
and start becoming parallel to each other below Tq ~ 470 K. 
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Chapter 1 



Introduction 



The cuprates exhibiting high-Tc superconductivity (HTS) have been receiving enormous at- 
tention in the scientific hterature, due to the great number of interesting physical effects they 
exhibit and to their technological applications such as Superconducting Quantum Interference 
Devices (SQUID) for highly sensitive measurements of magnetic fields, filters for the trans- 
missions with mobile phones, or electric power applications PP. Among the most studied and 
not yet completely understood issues are those connected with nonstoichiometric oxygen, its 
ordering and role in doping, and those connected with spin and charge inhomogeneities on the 
scale of nanometers, generally called stripes (see Ref. |2| for a recent review). 

In HTS cuprates, superconductivity sets in mainly in Cu02 planes doped with holes (or 
electrons in the case of l>fd2~xCexCu04^s) , and doping is due to the charge unbalance from 
aliovalent substitutional cations {e.g. Sr^+ in La2-xSrj;Cu04) and from nonstoichiometric oxy- 
gen (generally excess O^^). The oxygen stoichiometry may be varied over relatively wide ranges, 
but the amount of charge doping depends also on the ordering of nonstoichiometric O atoms, 
which are the most mobile atomic species. For this reason, the detailed knowledge of the dif- 
fusion and ordering mechanisms of oxygen in the cuprates are of great interest, especially in 
materials of the YBa2Cu306+a; (YBCO) family where doping is totally due to oxygen. A host of 
studies have been carried out on the oxygen mobility and ordering, especially with diffraction 
and permeation from gas phase methods, revealing complex phase diagrams and non trivial 
effects. Anelastic spectroscopy is one of the most powerful methods to study these complex 
phenomena, thanks to its ability of selectively measuring different types of hopping rates, from 
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isolated or differently aggregated O atoms, which produce different elastic energy loss peaks in 
the temperature scale. 

The topic of the intrinsic charge and spin inhomogeneities in HTS cuprates is also of great 
interest; in fact, not only the observation that the conducting holes may segregate into fluctu- 
ating stripes is difficult from the experimental point of view and counterintuitive, but it is also 
debated whether it is a phenomenon competing against superconductivity [2] or instead is at 
the basis of HTS [U E] • Rather unexpectedly, the anelastic measurements on cuprates of the 
LSCO family reveal also features attributable to the slow collective dynamics of charge stripes 
and antiferromagnetic domains in the Cu02 planes, and, to my knowledge, they are the only 
experiments where some of these dynamic processes are observable at acoustic frequencies. In 
fact, ac magnetic susceptibility, and dielectric, NMR, ^uSR spectroscopies are dominated by the 
single charge or spin fluctuations, while anelastic spectroscopy is insensitive to them and may 
probe the collective charge and spin motions through their weak coupling to strain. 

Two families of HTS cuprates, YBa2Cu306+x (YBCO) and La2-2:Sr2;Cu04 (LSCO), have 
been extensively studied by anelastic spectroscopy, by measuring the complex dynamic Young's 
modulus of ceramic samples at frequencies of 0.5-20 kHz between 1 and 900 K. In addition, 
some results are presented on the ruthenocuprate compund RuSr2GdCu208 (Ru-1212). All 
the samples have been obtained from a collaboration with the Department of Chemistry and 
Industrial Chemistry of the University of Genova, Italy (M. Ferretti), while several anelastic 
experiments have been done in collaboration with the Physics Department of the University of 
Rome "La Sapienza" (G. Cannelli, R. Cantelli, A. Paolone, F. Trequattrini). 

Several phenomena have been found and studied, among which diffusive hopping, ordering 
and off-center dynamics of O atoms, collective and local tilt dynamics of the octahedra in 
LSCO, and charge stripe fluctuations and depinning. Unless otherwise specified, all the results 
presented here were the first studies of such phenomena by anelastic spectroscopy. 

The Thesis is organized as follows. First is an introduction to anelasticity, limited to the 
concepts that are necessary for interpreting the phenomena discussed later, and emphasizing 
those concepts that are treated little in books on anelasticity; in particular, relaxation between 
energetically inequivalent states, interaction between elastic dipoles in the mean field approx- 
imation, and the relationship between anelastic and other spectroscopies. A brief description 
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follows of the method of measurement and of the samples treatments. A chapter is devoted to 
LSCO, with a short description of its structural and magnetic phase diagram, with the unstable 
tilt modes of the oxygen octahedra, and of the charge and spin stripes. The anelastic measure- 
ments are presented starting with the phase transformations, then interstitial O, followed by 
the newly found relaxational dynamics of the unstable tilts of the octahedra, and finally the 
observations of the thermally activated depinning dynamics of the hole stripes and their motion 
in the cluster spin glass phase, identifiable with the motion of pinned domain walls between 
antiferromagnetic domains. The following chapter is devoted to YBCO, starting with a pre- 
sentation of its complex structural phase diagram due to various types of ordering of O in the 
CuOx planes, and of the main results in literature on the mobility of this oxygen species. Sec- 
tion l5.4l is devoted to the diffusive jumps of oxygen and starts with an overview of the anelastic 
spectra at different stoichiometrics, followed by a discussion of what kind of effects one might 
expect from the interaction among oxygen atoms, at least in the simple Bragg- Williams ap- 
proximation. Then, the three distinct elastic energy loss peaks due to the oxygen diffusive 
jumps are discussed, with emphasis on the extremely fast hopping rate of the isolated O atoms 
in the semiconducting state, and on the role of the charge transfer between Cu-0 chains and 
Cu02 planes in determining various types of O jumps; also the ordering dynamics of oxygen 
is discussed. The chapter on YBCO terminates with a relaxation process identified with the 
hopping of O between off-center position in the Cu-0 chains, which diffraction studies suggest 
to be slightly zig-zag instead of straight, and with a peak which develops for x > 0.85, a doping 
range where all the physical properties are practically constant, and attributed to reorientation 
of pairs of holes (bipolarons) in the apical O atoms. 

The anelastic spectra of LSCO and YBCO reflect the great complexity of the structural, 
magnetic and charge phenomena occurring in these cuprates, and their interpretation is not 
always straightforward; therefore, at the end of the chapters devoted to LSCO and YBCO, a 
summary of the main results is provided, together with brief explanations of how the various 
anelastic effects have been assigned to specific mechanisms. 

Finally, a chapter is devoted to the diffusive hopping of oxygen in the ruthenocuprate 
RuSr2GdCu208-5, where it is more appropriate to talk of oxygen vacancies, since in the RUO2-5 
planes the oxygen stoichiometry is rather stable and close to the maximum. In this case, there- 



10 



fore, there is no complex phase diagram for the oxygen ordering, and the hopping dynamics may 
be described in terms of interaction among the elastic dipoles in the mean-field approximation. 
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Chapter 2 



Anelasticity 



The anelastic spectroscopy consists in the measurement of the complex dynamic comphance, 
or its reciprocal, the elastic stiffness, generally as a function of temperature at fixed frequen- 
cies. It is the mechanical analogue of the dielectric spectroscopy or ac magnetic susceptibility. 
Comprehensive treatments of the theory of anelasticity and of the application of the anelastic 
spectroscopy to the study of solids can be found in the seminal book of C. Zener [7], in the 
classic book by Nowick and Berry ^[ and in the most recent book edited by Schaller, Fantozzi 
and Gremaud The use of tensors for describing the elastic properties of solids can be found 
in books like Ref. 111! I12j . In the present chapter I will only mention what is strictly 
necessary for defining the notation and for the comprehension of the results discussed in the 
Thesis, with emphasis on the issues that are not treated in the above texts. 

2.0.1 Elastic dipole and thermodynamics of the relaxation 

For a perfectly elastic solid, Hooke's law can be written (in matrix instead of tensor notation) 
as 



where Si and ai are components of the strain and stress tensors and Sij the elastic compliance 
matrix, and summation over repeated indices is understood, or equivalently 
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where Cij is the stiffness matrix. In matrix notation, index i = 1,2,3 denote uniaxial strains 
along X, y and z respectively, while i = 4, 5, 6 denote shears of type yz, xz and xy (see Fig. 
I2-Ih ,.b: note that a shear strain is equivalent to two perpendicular uniaxial strains at 45° with 
different sign and equal magnitude, as shown by the gray arrows). Anelasticity results from the 
response to the application of a stress from defects or excitations that can change their state 
and also their contribution to the overall strain, with a characteristic relaxation time r; then, 
the elastic response according to the above equations is accompanied by a retarded anelastic 
response. 

For simplicity, I will refer to a molar concentration c of point defects uniformly distributed 
over the solid, and having at least two possible states; e.g. interstitial atoms that visit sites of 
type 1 and 2, according whether the first neighboring lattice atoms along the x ox y directions 

(Fig. \r^). 



(a) 




^ (b) 





defecl energ ■ 
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Figure 2-1: (a) Uniaxial and (b) shear strains, (c) Interstitial atoms and corresponding elastic 
dipoles under the application of a uniaxial cJi stress and (d) effect of cxi on the defects and 
strain. 

One then defines the concentrations ci and C2 of defects in states 1 and 2, with ci + C2 = c 
and their specific contributions to strain 



(2.3) 
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where a|"^ is the elastic dipole of the defects of type a. The definition "elastic dipole" derives 
from the analogy with the electric or magnetic dipoles of polar or magnetic defects |H], but, being 
a centrosymmetric strain tensor of the 2nd rank and not a vector, it is actually a quadrupole 
jllj representable as an ellipsoid with the principal axes defined by the tensor eigenvalues. This 
is shown in Fig. 12-lb for interstitial atoms causing greater lattice expansion in the direction of 
the nearest neighbor atoms. From Eq. H2.3|) it follows that the anelastic strain is 



^f; = E^A^f • (2-4) 



From now on I will drop the tensor indices unless necessary, and assume that a pure type of 
stress is applied, e.g. uniaxial, and the corresponding components of strain and compliance are 
probed. 

It is easy to show with a thermodynamic argument that the elastic dipole is also minus the 
rate of change of the elastic energy of a defect on application of stress. In fact, the differential 
of the Gibbs free energy g per unit volume (all the extensive variables are expressed per unit 
volume) is: 

dg = —eda — sdT H E^dca (2-5) 

a 

where vq is the molecular volume and s here is the entropy per unit volume; differentiating 
twice one obtains that 

A«.) = |1 = -pL. = -Ifi. (2.6) 
and therefore the application of a stress a changes the defect elastic energy as 

Ea{<7)=E^{0)-voX''a. (2.7) 

The defects occupy the possible states (or energetic levels) according to some distribution 
function and, for the sake of simplicity, let us considered diluted defects whose populations obey 
the Boltzmann distribution function; then 



^-Ea/ksT 

' z 



Z = J^e-^"/^B^. (2.5 
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The application of a stress a (see Fig. 12- lb and d) changes the energies Ea (c) of 

SEa = a—^ = -voX^^^a (2.9) 
oa 

and the new thermal equilibrium requires a repopulation of the states (with a relaxation time 
r) such that 



which results in the anelastic strain 



^an ^ ^ ^(.) = 5] ^ ||iAWA(/5) . (2.11) 

It can be shown that (the case of only two states is trivial) 

pairs a</3 

where Ua = Ca/c, which can also be extended to non-Boltzmann statistics ^Hl- The relaxation 
strength is defined as 

A = ? = 7 E '^(^'"'-^'^0'= E (2.13) 

pairs a</3 pairs a</3 

and this equation tells us that i) the anelastic response is the sum of the partial contributions 
Aa/3 from the repopulation of all the pairs of defect states; ii) Aq,^ is proportional to the square 
of the anisotropy of the elastic dipole (A^"-* — A*-^-*)^ and therefore elastically equivalent states 
(A(") = A*-^^) are not repopulated with respect to each other and do not cause relaxation; in) 
Aai3 is proportional to the defect concentration but also to the depopulation factor rianp, 
meaning that if the two states a and [5 are energetically inequivalent also in the absence of 
stress, then one of the two is less populated and this limits the stress-induced repopulation 
and the relaxation intensity. Note that in case of high density of defects the term n^n^ 
becomes of the type Ua (1 — ria) np (1 — np) since e.g. the jump of a defect from a site a to 
a site /3 is proportional to ria but also requires that site j3 is empty, hence the term (1 — np), 
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and analogously for the /3 — > a transitions. Such expressions are also symmetric in ria and 
1 — Ua, meaning that, e.g. when dealing with the jumps of the O atoms in the Cu02c or Ru02c 
planes of YBCO or Ru-1212 (0 < c < 1) for c —>■ the O atoms are the defects, but for c — > 1 
the O vacancies can rather be considered as the defects. Finally, the l/k^T term comes from 



dn 



^ oc {k^T) ^, which is a consequence of the fact that with increasing temperature all the 



defect states tend to become equiprobable. 

For the case of relaxation between only two states with 



AA = A2 - Ai , AE = E2- El (2.14) 



one obtains 

voc{AXf 



A (T) = '^^r ^ (2.15) 

4sA:BTcosh2(A^/2A;Br) 

which further reduces to the well known expression 

for equivalent states with AE = 0. The term cosh^^ {AE/2kT), coming from the depopulation 
factor nin2 in Eq. 1)2. 13() is generally overlooked, but becomes very important in all situations 
in which AE ~ k^T, since it produces a maximum in A{T) at T = 0.65AE/kB and then 
falls off as exp [—AE/k-QT) for T ^ 0. For example, when dealing with jumps of the O atoms 
in the CuOa; planes of YBCO, oxygen can pass from the isolated to the aggregated state and 
vice versa, which might well differ in energy by several tenths of electronvolt, which means 
thousands of kelvin in the temperature scale (I will often measure the energy in Kelvin, by 
considering E/k-Q instead of E] the conversion is 1 eV= 11600 K); the cosh"^ (^^E/2kT) term 
would cause a reduction of A (T = 500 K) by a factor 0.04 for AE ~ 0.2 eV and a factor 
0.004 for AE ~ 0.3 eV, making certain types of jumps completely unobservable in the anelastic 
relaxation. 

The depopulation factor becomes important at low enough temperature even for relaxation 
between states that ideally are energetically equivalent; in fact, unless dealing with extremely 
low concentrations of impurities in crystals of high perfection, there will always be long range 
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elastic interactions among defects that cause random shifts AE to the defect energies. Such 
shifts have been estimated for the case of 0-H pairs in Nb to be of the order of 100 K for 
impurity concentrations of ~ 1 at% jl4j . This means that, especiahy for disordered sohds hke 
the HTS, relaxation processes below 100 K are very likely affected by the depopulation factor, 
whose effect is of changing the temperature dependence of the relaxation intensity from 
to a nearly constant or increasing function of T. In such cases, I will include the depopulation 
factor in the analysis. 

2.0.2 Relaxation kinetics 

There is no general treatment for the relaxation kinetics and I will limit to the relaxation 
between two states from the start. In this case the rate equations for the defect populations 
are 



with Vij the rate for passing from j to i; thanks to the condition ni + n2 = 1 there is actually 
one independent nj. Setting hi = the equilibrium populations are found as 



which agree with the thermodynamic result if 1/12 e 1^21 satisfy the detailed balance principle 

^ = !l^ = eAW, (2.19) 

Incidentally, Eq. (|2.19j) requires that, for hopping over a saddle point according to the 
Arrhenius law, 




(2.17) 



1^12 + i^21 ' 



(2.18) 




(2.20) 



with the same tq for both states. Defining 



An = n2 — ni 
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the previous equations can be put in the form 



Are — An 

Are = (2.22a) 

r 

= {ui2 + V2i) = exp {-E/k^T) cosh (A^/2A;bT) (2.22b) 



which says that the rate for Are reaching the equihbrium value Are is proportional to its deviation 
from An through the relaxation rate r^^; the latter has been written in terms of the mean 
activation energy E = — {Ei + E2) and therefore contains a factor cosh (A£^/2A;Br) that 
should be taken into account when dealing with asymmetric states. 

On application of a time dependent stress o"(t), there will be an instantaneous elastic re- 
sponse 

e^i (i) = sua (t) (2.23) 
where subscript "U" stands for unrelaxed, and in addition an anelastic strain e^^ = c (niAi + n2X2) 



+ ^AAAn 



where the first term is constant and the time dependent anelastic response 



e^°(t) = ^AAAre(t) (2.24) 

with Are determined by Eq. H2.22a() . 
2.0.3 Dynamic compliance 

Let us calculate the time dependent anelastic response Eq. ()2.24|) on application of a periodic 
stress a = cioe*'^* that modulates AE. It is convenient to refer An to the equilibrium values in 
the absence of stress, defining: 

Are' = An - Are (a = 0) (2.25) 

and for the instantaneous equilibrium An = — tanh {/S.E /Ik^T) 

Are' = An - An (d = 0) = ^ , a = ^- ■ -a (2.26) 

^ ^ dAE da 2kBT cosh^ {AE 2kBT) ^ ' 
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By substituting into Eq. ()2.22a() 



An' 



IkBTcosh"^ {AE/2kBT) 1 + ilot 
and using Eqs. (|2.2l-{l2.24j) the dynamic compliance can be written as 



(2.27) 



s{uj) 



I ■ II 

S — IS 



a 



1 + 



A(T) ' 
1 + iujT 



(2.28) 



where A (T) is the same as given by Eq. ()2.15|) . The real and imaginary parts of s{uj) are 
therefore 



s' (w,r) = s\] 



1 + 



A(T) 
1 + [lotY 



(2.29) 



/(w,r) = suA(r) 



1 + (wr)^ 



where the s" presents the well known Debye peak of amplitude (T) at the condition for 
maximum relaxation 

WT = 1 , (2.30) 

while s' a step of amplitude A (T) . In the lot ^ 1 limit the defects are too slow to follow 
the periodic stress and corresponds to the elastic limit without defects {s' = su, s" = 0). In 
the LOT <ti 1 limit the defect relaxation is so fast that instantaneously complies to the periodic 
stress, s" = 0, and the compliance is totally relaxed, s' = su (1 + A). 

It is possible to make an analogous derivation for the dynamic stiffness 



c {lo, T) = c' + ic" 



(2.31) 



finding 



cu 1 



A(T) 
1 + (lotY 



c" = cuA (T) 



LOT 



1 + (lotY 



(2.32) 
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2.1 Elastic energy loss and anelastic spectrum 

The fact that the dynamic comphance s {uj,T) or stiffness c{lo,T) is complex means that, due 
to the retarded anelastic response, strain is out-of-phase with respect to stress by the loss angle 

s" c" 

tan(p = — = — . (2.33) 

s' c 

Generally, it is A (T) ^ 1 so that 

tan — = A (T) ""J . (2.34) 

The tangent of the loss angle can be measured from the dissipation of elastic energy or acoustic 
absorption; in fact, if cr = do cos {tut) and e = Eq cos {ujt — (p), the elastic energy dissipated in 
one vibration cycle is 

AW = / ade = vreoCo sincp 

Jujt=0 

while the maximum elastic energy stored is W = J^^Zq^'^ o'de' = |eoO"o cos (j), where e' = 
Eo cos (j) cos (iot) is the strain component in phase with a. One defines the elastic energy loss 
coefficient as 

1 AW 

Q-^ = —-— = tan^, (2.35) 
zvr W 

which coincides with the reciprocal of the mechanical Q of the sample. 

Generally, the anelastic measurements are made sweeping temperature at nearly fixed fre- 
quency Lo, and the resulting spectrum contains absorption peaks in correspondence with the 
temperatures Tmax for which the condition of maximum relaxation ujt (Tmax) = 1 is verified, 
while c' (s') contain negative (positive) steps. This is shown in Fig. I2-2| assuming that the re- 
laxation rates follow the Arrhenius law, = Tq^ exp {E/k-QT), between states with the same 
energy, so that A (T) oc 1 /T. Figure l2-2b shows how slower processes are peaked at higher 
temperature, and how it is possible to estimate the energy barrier E by measuring at different 
frequencies, exploiting the condition r (Tmax) = at the peak temperatures. 

It is more convenient to analyze the {uj, T) curves rather than the real parts s' {uj, T) or 
c' (w, r), especially if A is very small, because the absorption peaks stand out of a usually small 
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Figure 2-2: (a) Contribution of a Debye relaxation process to the imaginary and real parts of 
the dynamic stiffness at two measuring frequencies, (b) Anelastic spectrum with two processes 
having the relaxation rates plotted in the upper part; measured at the two frequencies wi and 

background (at least for w < 10 MHz, when anharmonic effects and sound wave diffraction 
at the boundaries are small); instead, even for a perfectly elastic solid, the real part or 
cu has a temperature dependence due to anharmonic effects that should be subtracted in 
order to analyze the anelastic contribution. In addition, the elastic moduli are affected by the 
sample porosity, by microcracks, and internal stresses due to anisotropic thermal 
expansion building up during thermal cycling, which may result in anomalies and hystereses 
on varying temperature. 

Figure presents an example of anelastic spectrum measured on LSCO at three resonance 
frequencies during the same run; the stiffness which is measured is the Young's modulus E. 
Note that the scale of the real part is almost 10 times larger than that of the absorption. Note 
also that the amplitudes of the steps in the real part are larger than 2A deduced from the 
amplitudes A of the corresponding absorption peaks, as one expects from Eq. (|2.29|) . This is 
due to the fact that the relaxation processes are broadened by distributions of relaxation times; 
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the elementary peaks are shifted with respect to each other, so that the resulting peak resulting 
step in E' is the sum of the elementary steps. 
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Figure 2-3: Example of anelastic spectrum (real and absorption parts) of LSCO measured at 
three vibration frequencies during the same run. 

The pre-exponential factor in the Arrhenius law is the relaxation rate extrapolated to 
infinite temperature, and its value deduced from experiments at low temperature should not 
be taken too seriously, since one should also take into account the temperature dependence of 
all the quantities that affect the jump rate, including the vibration entropies. Still, as a rule of 
thumb. To ^ 10^^^ s, of the order of magnitude of the local vibrations promoting the atomic 
jump, is indicative of point defect relaxation, while tq > 10^^^ s is indicative of extended defects 
or collective motions. 

2.2 Elastic energy loss and spectral density of degrees of free- 
dom coupled to strain 

The dynamic susceptibility x is defined as the ratio between a response dr and the excitation 
force df: x = fjj in the elastic case, the compliance s = is the elastic susceptibility with 
f = (J and r = e. The fluctuation-dissipation theorem correlates the imaginary part of a 
susceptibility x" with the spectral density J (tj, T) (Fourier transform of the autocorrelation 
function) of the spontaneous fluctuations of r; for the elastic case and in the classical limit it 



22 



can be written |151 116j : 

s"{co,T) = ^ J dte-*(e(t)e(0)) = ^J(u;,r) (2.36) 

where (...) denotes the thermal average and V is the sample volume. This important equa- 
tion tells us that the elastic energy absorption at angular frequency uj is proportional to the 
corresponding Fourier component of the spontaneous strain fluctuations, which are due to any 
motion or excitation coupled to strain e (as e.g. in Eq. 1)2. 4jl ). For a process with relaxation 
time r, also called of the diffusive or pseudodiffusive type, the strain autocorrelation function 
is {e (t) e (0)) =< fe^ > exp (— |t| /r), and its Fourier transform is 

Jiu;)= <6e^> ] (2.37) 

which, introduced into Eq. (|2.36|) yields the usual Debye formula (|2.29|) with the ther- 
modynamic factor ()2.16|) . This formulation of s" {io,T) and therefore of {io,T) will be 
useful when comparing anelastic and NQR data in Sec. 14.91 and in general when comparing an 
anelastic spectrum with those from other spectroscopies {e.g. neutrons or infrared). In fact, the 
intensities of the transmitted or diffracted radiations are proportional to the spectral densities 
of the fluctuations of the physical quantities x coupled to those radiations (nuclei positions for 
neutron scattering, electric polarization for light, etc.). A J (uj) of the form of Eq. (|2.37j) is 
also called a "central peak" for the following reason. Motions x (t) with resonance frequency 
ujQ and mean lifetime r (vibrations, excitations with energy hujQ, tunneling) have a correla- 
tion function {x (t) x (0)) oc cos exp (— \t\ /r), and their spectral density is a lorentzian peak 
J (uj) = J")i_)_^-2 centered in and with width t~^; Eq. 1)2.37(1 can be obtained by setting 
loq = 0, and therefore is a peak centered at the origin of the frequency or energy scale. 

This is schematically shown in Fig. I2-4| where the physical quantity x{t), e.g. atomic 
positions coupled to strain e and to the electric field gradient producing NQR relaxation (see 
Sec. 14. 9() . displays resonant and pseudodiffusive types of motion. The spectral density of the 
fluctuations of x contains both the peak centered at loq > 10^^ s~^ (local vibration in a potential 
minimum) and the central peak of width (hopping between two potential minima with rate 
r^^). At acoustic frequencies the resonant process has negligible spectral weight and only the 
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Figure 2-4: Spectral densities J (w) of x {t) for resonant (blue) and relaxational or diffusive 
(red) processes. The imaginary elastic compliance or energy loss is related to J (w) through the 
fluctuation-dissipation theorem. 

central peak is observed; s" can be obtained from Eq. (|2.36j) . while the NQR relaxation rate is 
given directly by J (co) [T7| . 



2.3 Distributions of relaxation times 

Up to now we considered anelastic relaxation from equivalent non-interacting elastic dipoles, 
all having the same relaxation time r, as is the case of very diluted solid solutions. More 
common are the situations with distributions of relaxation times ^(t), due to static disorder 
or interactions among the dipoles, which result in broader peaks in s" {uj,T). 

In what follows we are mainly concerned with the imaginary or absorption part of the 
compliance or modulus, because it is easier to analyze than the real part. In fact, the imaginary 
part s" {uj,T) associated with a relaxation process with characteristic time r is peaked around 
the temperature at which ujt (T) = 1 holds, and is superimposed over a background contribution 
•^bkg ("^i ^) ■which is usually small and slowly varying over the temperature range of the peak in 
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■s" {uj,T). Instead, the contribution to s' {oj^T) extends well above the temperature for which 
(T) = 1, and is superimposed to elastic compliance s^j (T), which is large and temperature 
dependent. This is at variance with the dielectric case, where the equivalent of s'^j (T) is close 
to the vacuum permittivity eo- 

2.3.1 Uniform distribution of activation energies or Inr 

The simplest distribution of relaxation times with a clear physical meaning is a uniform dis- 
tribution in the activation energy E: g {E) = {E2 — Ei)~^ for Ei < E < E2. Setting 
Inr = Inro + E/T, this is equivalent to a uniform distribution in Inr: Inro + Ei/T < Inr < 
lnTo + ^2/T: 

g (E) dE = = ^^^"[^ = g (In r) d (In r) (2.38) 

E2 — El m T2 — m Ti 

The imaginary part of the compliance can be easily integrated over such a distribution, 
where we include a relaxation strength cx 1/T [see Eq. (|2.1tj|l ]: 

dE As 2u;Toe^/^ 



{E2-E1) T i + (^roe^/T)2 

= , [arctan fwroe^^/T'j _ Cretan (toToe^'/^)] . (2.39) 

[E2 — El) I \ J V /J 

This distribution, however, usually provides a poor fit to very broadened processes, since it 
develops a plateau at the maximum. It has been used to describe glassy processes, by letting 
T2 (T) diverge at the freezing temperature, and assuming a constant relaxation strength. The 
omission of the 1/T term in the relaxation strength transforms the broad plateau into a linear 
temperature increase from the freezing temperature up to the maximum at ujt2 (Tm) ~ 1, which 
describes reasonably well, for example, the dielectric maximum of some relaxor ferroelectrics 
[TH] , possibly by smoothing the cutoff at T2 P^ . 
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2.3.2 Fuoss-Kirkwood distribution 

Fuoss and Kirkwood 20 showed that, when the imaginary part of a susceptibihty x" rnay be 
written as 



X"(a;)= / dilnr) gilnr)- 
Jo 1 



LOT 



(2.40) 



then the distribution function g of the relaxation times may be expressed by analytical contin- 
uation in terms of x"- 



gilnr) = - 



IT 



X 



(2.41) 



where = 0J~^ at the maximum of x" is the mean value of r, when g (In r) is even in In r/xm = 
Inr — Inr^ and therefore x" is even in InijOTm- In this manner, it is possible to associate a 
distribution function g (In r) to any analytical form of broadened peak in s" (w, T). A commonly 
used expression is 



s" (c^) 



a (a;r)" 

?2a 



a 



(2.42) 



1 + 2 cosh [a In (ljt)] 

with 0<Q!<1;q; = 1 corresponds to a Debye peak, whereas a < 1 broadens the peak by a 
factor in the Inr and therefore T"^ scale (and lowers its amplitude by a); the corresponding 
distribution, called Fuoss-Kirkwood distribution [S], can be calculated through eq. 1)2. 41|) : 



5FK (In r) = - 



vr 



a 



a 



-+-. 



a 



2 cosh [— Q In (t jTm) -|- iaf ] 2 cosh [—a In (t jTm) — ^"f ] 
cosh [a In r/Tm] cos ^ 



TT cos^ ^ + sinh [a In r/ r^] 
and is normalized to 1; g-pY. (2;) can also be well approximated with a Lorentzian 



(2.43) 



5FK ix) 



a 



^cos^f x2 + (u;/2)^ 



with w = 3.38 X (. 



-1 



1) 



(2.44) 



for a down to 0.4. If Inr = Inroexp [E/k-QT), and the effect of a distribution in the values of 
To is neglected in comparison to those of E, then g-pK (In ''") can be attributed to a distribution 
in E, roughly Lorentzian with a temperature dependent full width at half maximum we = 
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3.38 X (a-i - 1) kBT. 



2.3.3 Cole-Cole distribution 

The simple expression 

^ 1 // ^ (^r)"sin(§a) ^ 1 sin (fa) 

l + (ia;r)°' 1 + (u^t)^" + 2 (l^t)" cos (f a) 2 cosh [a In (i^r)] + cos (f q) ^ ' ^ 

is very popular in the analysis of the dielectric susceptibility. The corresponding distribution 
function is 

1 sin (vtq) . 

gcc In t) = - ■ . . . / ./^ . (2.46) 

ZTT cosh [a in (r/rmjj + cos (vraj 

2.3.4 Other expressions 

Other expressions for xi'^) ^-nd corresponding distributions in g{lnT) are used, especially in 
the dielectric literature, for example the Havriliak-Negami one, 

[1 + [lUJT) \ ' 

To my knowledge there are no particular physical reasons for preferring one over the others, 
except the empirical fact that different processes may be better interpolated by different ex- 
pressions. For example, a process whose x" i^^T) is not even in 1/T will be better interpolated 
by an expression of the Havriliak-Negami type, which is not even in In lot, but also a simpler 
expression, like 

[UT) + (UJT) ' 

may be used |21j . The latter is a generalization of the Fuoss-Kirkwood expression, where 
rescaling in 1/T by a is adopted for the low-T/high-r region and rescaling by /? in the high-T 
region. An even more flexible expression was proposed by Jonscher j22j . 

^" = , ,a ' (2-49) 

(WTi) + {UJT2) 

where ri (T) and T2 (T) are two relaxation times, possibly following the Arrhenius law, which 
describe the relaxation of the system at long and short times, respectively. 
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Figure 2-5: Normalized x" vs Inw according to the Fuoss-Kirkwood and Cole-Cole expressions 
for different values of a. 

2.4 Interacting elastic dipoles 

The treatment of interacting elastic dipoles is of particular importance in the case of oxygen 
in the CuO^ planes of YBCO, where the concentration of dipoles 0.1 < x < 1 is by no means 
small, and their interaction is so strong to give rise to a complicated phase diagram. Under 
such conditions, any attempt at describing the anelastic relaxation from the hopping of the O 
atoms should somehow take into account their mutual interactions, which are both of electronic 
and elastic origin. The description of the structural phase diagram of YBCO, however, requires 
rather sophisticated models with asymmetric interactions (different along the a and b axes) at 
least up to the next-nearest neighbors (so-called ASYNNNI models [2'A\ ) . Such models can be 
solved only with Monte Carlo techniques, and are generally adopted to reproduce the YBCO 
structural phase diagram, in only few cases to evaluate the oxygen diffusion coefficient |241 125j , 
and no attempt exists to evaluate the dynamic elastic compliance due to oxygen hopping. The 
numerical results on the tracer diffusion coefficient |24[ I25j are of little help in analyzing the 
anelastic data. 

Wipf and coworkers [^H] carried out an analysis of the high temperature anelastic measure- 
ments of YBCO where the interaction between the O atoms is treated in the Bragg- Williams or 
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mean-field approximation; in that manner, the complexity of the YBCO phase diagram cannot 
be obtained, since only a single ordering phase transformation is reproduced, identified with 
the tetragonal to orthorhombic one, but useful expressions of the dynamic compliance can be 
obtained. This treatment is particularly appropriate to the analysis of anelastic relaxation from 
oxygen jumps in the Ru02-(5 planes of Ru-1212, which can indeed be considered as a diluted 
solution of O vacancies with long range elastic interactions (see Sect I6.2j) . 

A treatment of the dynamics of interacting elastic dipoles had also been carried out previ- 
ously by Dattagupta [271 128j in the mean field approximation, assuming that the actual elastic 
field that a particular dipole senses is substituted with a mean stress field auF due to all the 
other dipoles, the same for all dipoles. Such a treatment has also been reviewed in connection 
with the Snoek relaxation of O interstitial atoms in bcc metals j29j . Dattagupta's model does 
not explicitly take into account the concentration of relaxing dipoles, since it assumes that there 
is one dipole per cell, which is able to change between three possible orientations; it is the elastic 
analogous of the Curie- Weiss theory of interacting magnetic dipoles, and the result is essentially 
the same. If the elastic interaction E = Acmf among the dipoles on the average favors parallel 
orientations of the major axis of A, AE = E± — i^y > 0, a cooperative alignment of the dipoles 
occurs, which results in an increase of both the relaxation strength A (T) and relaxation time 
r by a factor (1 — Tq/T)^^ . The effect of the cooperative motion can be roughly viewed as 
the coordinated motion of several dipoles instead of independent dipoles, which results in a 
larger effective elastic dipole, but also in a slower reorientation time. When the anisotropic 
component of the interaction energy, A£^, is of the order of the thermal energy k-^T, the dipoles 
start freezing into a fixed orientation, resulting in a ferroelastic ordering transition, similar to 
the ferromagnetic or ferroelectric one. The expression of s" above Tq becomes 
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,11 



OJT 



S' 
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2 



(2.50) 



A 



A (T) Tc 



t{T) 



(2.51) 



1 - Tc/T T-Tc 



T = 



1 - Tc/T 



and, assuming A (T) = Aq/T can be put in the form 



s 



Aq cot 



s' 



r (1 - TjTf + (wr)2 



(2.52) 
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where r is the relaxation time for non-interacting dipoles. Such a treatment might be adapted 
to only two orientations, as is the case of oxygen in YBCO, and the expressions of the elastic 
compliance and relaxation time might be extended below Tq', the dependence on the dipole 
concentration c might be taken into account as an interaction energy AE ~ Actmf oc (AA)^ c, 
giving rise to k^Tc ~ AE oc c, as also suggested in However, the result would be 

correct only in the limit of low concentration, and in what follows I will adopt Wipf 's approach 

m- 

2.4.1 Ordering transition 

Let us consider the CuO^ plane of YBCO, where, like the interstitial atoms in Fig. 12-lb . the 
O atoms can occupy sites of type 1 and 2, corresponding to the positions generally called 0(1) 
and 0(5) (see Fig. 15-11) : in the completely ordered phase only sites of type 1 are occupied. The 
starting point of Wipf 's analysis is the equilibrium condition for the O atoms between the 
sites of type 1 and 2: 

/ii (ci,c,r) = /i2 (c2,c,r) (2.53) 

where fia is the chemical potential of the O atoms in the sublattice a and c = | is the fraction 
of sites populated in the CuO^ plane (0 < c < 0.5 for YBCO) with 

ci + C2 = c. (2.54) 

The chemical potentials are evaluated from the free energy F = U — TS, as 

OF 

and in the Bragg- Williams approximation is F = E — kBTlnAf, with E the average energy of 
the dipoles and M the number of ways a state with that energy may be obtained, disregarding 
the interactions among dipoles. Therefore, if the crystal contains N cells, M is the number of 
ways A''i O atoms can be put in sites of type 1 and O atoms in the other N sites of type 
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2: 

N{N-1) ... (N-N^ + 1)N{N- 1) ... {N - + 1) tT .9 .fi^ 

0=1 ^ ' 

and using Stirling's formula InA^! ~ A^lnA^ and Cq = Na/N 

2 

- In = N^[co, In c« + (1 - c«) In (1 - Ca)] (2.57) 

a=l 

The mean energy is 

N 

E <^-^/3 (2-58) 

i,j<i=l 

where in the mean-field approximation j28j the interaction energy between two dipoles of type 
a and /5 in the cells i and j is set 



I E\\ if a = 3 
<= " (2.59) 
' Ex ifa//3 



yielding a mean energy per unit cell 
where 

AS = E|| - E± (2.61) 

and AiiJ < favors parallel orientations. The free energy per unit cell can then be written as 

2 

F = E- kBTlnM = /^E^cl + E^c^ + k^T ^ [c„ lnc„ + (1 - c^) In (1 - c„)] (2.62) 

cc 0=1 

and the chemical potentials 

^la = ^ = 2AEca + kBT In ( ) (2.63) 

OCrv \ 1 Cry / 
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or, adopting the notation of |26j . 

;.. = |£ = -a(c„-|)+.BTln(^) (2.64) 

where the energy is referred to that of the completely disordered state, Cq = |, and a = 
—2AE > is the difference in energy between states 1 and 2, supposed to be weakly dependent 
on temperature. By equating /ii and ^2 we obtain 



" (ci - C2) = In 



kBT 



ci (1 - C2) 

C2 (1 - Cl) 



(2.65) 



having solutions different from the trivial one ci = C2 = ^ only below the critical temperature 

kBTc = a^(l-^); (2.66) 



2 V 2 

below Tc, the equation can be solved numerically for ci, as shown in Appendix A. 
2.4.2 Coupling to stress and relaxation strength 

On application of a stress aki, the elastic energy of each O atom of type a changes by 

-voX'kl^m , (2.67) 

where for tetragonal elastic dipoles A^^ = = Ai , Ayy* = X^x = A2, while the z component 
has no interest since it is the same for both types of sites, so that the anelastic strain due to 
the oxygen dipoles is (compare with Eq. (12.241) 1 

£xx = -Syy = ;^ (ci - C2) (Ai - A2) . (2.68) 



The change in elastic energy, given by Eq. (|2.9|) . results in a change of the populations from 
Cq to c'q, (ct) such that 

^J'l {c'i,c, T) - voX^^^Uki = fJ-2 (4, c, T) - vqA^^Vh (2.69) 
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and expanding to first order in — Ca and using the fact that (c'^ — ci) = - 

[■^kl -^kl ) '^kl _ (Ai - A2) (o-^^x - (Tyy) 



(C2 - C2) 



Ci - Cl 



dci dc2 



9ci (9c2 



which results in an anelastic strain 



eT. = -eZ = ^ [(c'l - c'2) - (ci - C2)] (Ai - A^) 



or relaxation of the compliances Ssijki = e^J /aki'- 

vq (Ai - 



VQ (Ai - A2) {OxX - CTyy) 

dfJ-i _|_ dfi2 
dci ' dc2 



5sii = SS22 



-5s 
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dfii _|_ (9/J2 
dci dc2 



FliT) 



vq (Ai - A2)' 



where the factor Fi (T) is computed from eq. (12.641) : 

'1 



Fi(r) 



c(2-c) ^ c(2-c) 



Cl (1 - Cl) C2 (1 - C2) 



For T > Tc it is Cl = C2 = c/2 and therefore 



Fi{T) 



T 



T-Tc 



{T > Tc) 



(2.70) 



(2.71) 



(2.72) 



(2.73) 



(2.74) 



and 6s becomes the Curie- Weiss susceptibility; below Tc, -Fi (T) falls off to zero, as shown in 
Fie. E 



2.4.3 Relaxation rate 

Following l^n] , the rate equation for the instantaneous out-of-equilibrium concentrations c" and 



Cn = c — c'l' is 



dt 



1 



Tq-i exp i-Eo/kBT) k (1 - c'l') c^^/^'^bT _ ^// _ ^//) ^-AE/2k^T] (2.75) 
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2^ L 



:'2' (1 - c'l') e"«-^)/^BT _ ^// (1 _ ^/^/) ^-a{c'{-c)/k^T 
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where = Tq^ exp {—Eo/k-QT) is the relaxation rate when sites 1 and 2 are equivalent (ab- 
sence of stress and disordered state), the two terms are the rates from 1 to 2 and vice versa, 
taking into account the probability that the neighboring site is empty, while the ^^^^I'^^'bT 
factors account for the non-equivalence introduced by the ordering transition, where the differ- 
ence in the site energies is given by the difference in chemical potential, without counting the 
entropic term: 

AE = E2-Ei = a {c'l - c'^) = 2a (c'/ - c) . (2.76) 

On setting dci/dt = one finds again eq. (|2.65j) . which demonstrates the consistence between 
the thermodynamic-statistical and the kinetic analyses. 

In order to find the dynamic response to a small applied stress, we set c'/ = ci + 5c and 
C2 = C2 — 6c into the rate equation, where 5c (t) is the small response to the applied stress, and 
keep terms to the first order in 5c; after some algebra one obtains 



d5c 6c 

dt Teff 



(2.77) 



with 



Teff = T (^) = (T) . (2.78) 

y'CiC2 (1 - C2) (1 - Ci) 

For T > Tq the factor before -Fi becomes unity, and one finds that 

Teff = ^ ^^-^^^ 

which represents the critical slowing down on approaching the phase transformation. The 
functions Fi and F2 that enhance the relaxation strength and time near Tq are shown in Fig. 
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Chapter 3 

Experimental 



3.1 The samples 

The YBaaCusOe+x (YBCO), La2-xSr^Cu04 (LSCO) and Nd2Cu04+5 (NCO) samples have 
been prepared by M. Ferretti at the Department of Chemistry and Industrial Chemistry of 
the University of Genova, Italy. The powders of the starting oxides were mixed in appropriate 
amounts and calcined, typically 1233 K in air for 6 h for YBCO, 1323 K for 18 h for LSCO 
and 90 h for NCO. The powders were then pressed into bars 5 x 5 x 50 mm^ and sintered in 
fluent O2 (typically YBCO at 1253 K for 22 h and LSCO at 1323 K for 18 h), and further 
oxygenated in fluent O2 at lower temperature for the case of YBCO. The powders were checked 
with X-ray diffraction for impurity phases from decomposition or incomplete solid state re- 
action, and sometimes intermediate grinding and sintering treatments were made. Especially 
RuSr2GdCu208 (Ru-1212) required a long processing route [SO] before the final oxygenation at 
1343 K in fluent O2 for one week, followed by cooling at 50 K/h. 

The ingots were finally cut with a circular diamond saw to the dimensions suitable for the 
anelastic measurements: length of 40 — 50 mm, thickness of 0.2 — 1.5 mm, width of 4 — 5 mm. 
The samples were rather porous (from 10% up to 50% for Ru-1212) and some of the oil used 
during cutting remains in the pores and it is almost impossible to completely remove it with 
solvents (toluene and acetone) in a ultrasound bath. The complete removal of the oil can be 
accomplished by burning it, heating the sample at high temperature. The residual oil, which 
can also be adsorbed from the pumping system ^31^, has a freezing transition around 220 K 
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giving rise to an absorption peak and anomalies in the elastic moduli (see also Sect 15. 7| . 

Further characterization of the samples was made by measuring the resistivity and super- 
conducting transition with a four-probe technique with a closed cycle cryocooler from 300 down 
to 12 K. Checks were made that resistivity p did not depend on the amplitude and direction of 
the applied voltage. The oxygen stoichiometry in YBa2Cu306+a; was estimated by combined 
analysis of the p (T) curves (compared with literature data) and of the anelastic spectra; in some 
cases series of samples have been prepared with different stoichiometrics, checked by iodometric 
titration and from the cell parameters from X-ray diffraction. 

3.2 Anelastic measurements with resonating samples 

The equation of the vibration of a solid near one of its resonances is that of a damped harmonic 
oscillator with resonance frequency |H] /o = ot\/ ^/ P = 2w "^^ere M is an effective elastic 
modulus and a a geometrical factor depending on the vibration mode. Then, the equation of 
motion for strain under application of a stress cioe*'^* with u; ~ is 

CToe*'^* = {M' + iM") e + pe = pujl (1 + i tan c^) e - pu?e (3.1) 

where the real and imaginary parts of the modulus have been introduced. 

All measurements presented here have been made by exciting the samples on their 1st, 3rd 
and 5th flexural modes. Thanks to the sample dimensions, with length 30 — 800 times larger 
than the thickness, and to the low vibration amplitude, always within the linear regime, the 
approximation of simple bending of thin bars is a good one. Under such approximation the 
modulus involved in the flexural vibrations is the Young's modulus that relates uniaxial 
stress and strain of a long sample of uniform cross section. In the case of pure extensional 
vibrations, where stress and strain obey a = Ee all over the sample, the resonance frequency 
of the n-th mode is given by /„ = (n/2/) ^ E j p where / and p are the sample length and 
density (neglecting a correction factor that takes into account the finite size of the cross-section 
with respect to the length and proportional to the square of Poisson's ratio |S21)- In the case 
of flexure, strain is inhomogeneous along the sample thickness, passing from extension on the 
convex side to compression on the concave one, as shown in Fig. 13-11 The effective modulus 
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(the real part) is therefore E, but reduced of a geometrical factor taking into account the sample 
cross section; for the fundamental mode of a thin bar of thickness h, it is |H] 

h 



/ = i.03^/b7p; (3-2) 

the frequencies of third and fifth modes are respectively 5.40 and 13.3 times larger. 



nodes = sin(f()t4) 

Figure 3-1: Electrostatic excitation of the flexural vibrations of a thin bar. Strain is inhomo- 
geneous along the sample thickness: from a uniaxial expansion (blue) to a compression (red). 

We never tried to estimate the absolute value of E because of irregularities in the samples 
shapes and above all due to the large sample porosity (from 10 to 50%), but it is possible to 
measure the variation of E (T) with respect to a reference temperature Tq through 

E{T) ^ f (T) 

E{n) /2(To)' 

where the temperature dependence of p (T) is neglected in comparison to that of E (T). 

The elastic energy loss can be measured with two methods. In the free decay method, 
the excitation is switched off and, according to Eq. 1)3. 1|) with ao = 0, the vibration amplitude 
falls off as 

Therefore, interpolating with a straight line the logarithm of the vibration amplitude A (t) = 
eoex.p {—(f)u>()t/2), log [A (t) /A(0)] = —^Q~^u>Qt, it is possible to deduce tancf) = Q~^. This 
method is useful when the decay constant (5/('^/o) is long enough to let the amplitude decay 
to be recorded. 

For larger dissipations or frequencies the forced oscillation method can be used, where 
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the vibration amplitude is measured on sweeping frequency near resonance. From Eq. p.lf) . 
the vibration amplitude is 



(wg — u)'^) + cJq tan^ 4> 

and by least square fitting of yl^^ (oj^) one obtains the resonance frequency ujq/2tt and the 
elastic energy loss coefficient Q~^. 



3.3 The anelastic relaxation setup 

The anelastic measurements were based on resonant sample techniques developed in the Lab- 
oratory of Solid State Acoustics at the Istituto di Acustica "O. M. Corbino" of CNR, starting 
from Prof. P. G. Bordoni up to Profs. G. Cannelli and R. Cantelli. The sample is suspended 
on thin thermocouple wires in correspondence with the nodal lines for the flexural vibrations; 
generally the nodes at 0.225/ from the sample ends (/ = sample length) are used, since they are 
very close to nodes of both the 1st and 5th flexural modes, whose resonance frequencies are in 
the ratio 1:13.3. It is therefore possible to measure these two frequencies during the same run. 

As schematically shown in Fig. 13- H an electrode is kept very close the sample surface and the 
sample, if not conducting, is covered with silver paint in correspondence with the electrode and 
suspension wires, in order to have electrical continuity and proper grounding (through one of the 
thermocouple wires). An excitation ac voltage 14xc is applied to the electrode with frequency 
//2 and amplitude up to 200 V; this induces a potential of opposite sign on the surface of the 
sample and the resulting electrostatic force oc Vg^c excites vibrations of the sample at frequency 
/ (100 Hz-lOO kHz depending on the sample shape and vibration mode). These vibrations 
modulate the sample/electrode capacitance, which is part of a high frequency (~ 10 MHz) 
circuit, and therefore modulate its resonant frequency, which is in turn demodulated with a 
coupled resonating circuit with tunable capacitance. This revelator, called vibrometer, was 
constructed in the Solid State Acoustics Laboratory by Profs. Nuovo and Bordoni. I improved 
the data acquisition by introducing a lock-in amplifier that measures the signal amplitude and 
can detect low vibration levels at high noise level, and developed the acquisition software. Now 
it is routine to perform measurements at several frequencies during the same run sometimes in 
semiautomatic manner. This is very important whenever the status of the sample is modified 
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by the measurement itself, for example due to oxygen loss in vacuum at high temperature, and 
a subsequent run at a different frequency could not be compared with the previous one for 
making an analysis at different frequencies. If the sample is properly mounted and its intrinsic 
dissipation is Q^^ ~ 10~^ or larger, it is possible to measure several flexural and torsional 
modes with the same mounting; Fig. 12-31 is an example with 1st, 3rd and 5th flexural modes 
vibrating at 1.3, 7.2 and 18 kHz respectively. 

There are two separate inserts for low (1-300 K) and high temperature (100-950 K) mea- 
surements, evacuated by a diffusion pump; He exchange gas is used for thermalization, and an 
adsorption pump is used at low temperature, in order to minimize the adsorption by the sample 
of residual gases when the diffusion pump is closed (being cooled with liquid N2, the adsorption 
pump does not condense the He exchange gas). 

The low temperature insert has a cylindrical brass cell that can be closed both with a conical 
coupling with vacuum grease to the flange or with In wire for measurements below 2 K. The cell 
is equipped with heater and silicon diode for temperature control and a wound copper tube for 
the flow of LHe. The whole insert may be closed with a vacuum tight glass dewar and works as 
a LN bath cryostat and LHe flow or bath cryostat: below 5 K the dewar is filled with LHe and 
may be pumped in order to lower temperature down to 1.1 K. The thermocouple is Au-0.003Fe 
versus chromel. 

The high temperature insert is made of stainless steel and ceramic insulating parts and is 
closed with a quartz tube. Temperature is regulated with a tubular furnace or dewar with LN. 

3.4 Reliability of the anelastic spectroscopy measurements 

Sometimes, there is scepticism on interpretations of anelastic relaxation processes outside mech- 
anisms based on atomic defects, dislocations, domain boundaries and grain boundaries; in ad- 
dition, it is often felt that the anelastic spectroscopy should be particularly sensitive to the 
polycrystalline nature of ceramic samples and to grain boundaries, making the interpretation 
of the anelastic spectra very unreliable. In this section few arguments are put forward demon- 
strating that such doubts are generally unjustified. 

Regarding the effect of impurity phases, generally found at the grain boundaries, it should 
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be remarked that anelastic spectroscopy probes the whole volume of the sample; therefore, 
if defects motions or phase transformations occur in an impurity phase, their effect will be 
weighted with the molar fraction of the phase (and another factor of order unity taking into 
account the different moduli of bulk and impurity phase |33p. Therefore, spurious phases 
may affect the anelastic spectra exactly as they affect other bulk measurements like specific 
heat or neutron spectroscopy, and much less than techniques sensitive to the grain surfaces, 
like resistivity, dielectric susceptibility or optical spectroscopies. In the samples considered 
here the molar fraction of impurity phases did never exceed few percents, as checked by X-ray 
diffraction, which means that their possible contribution to the anelastic spectrum was rather 
small. Of course, if the impurity phase presents a phase transformation, even small amount 
of it may produce noticeable effects on the anelastic spectrum, as for example occurs with 
Lai-i-'Ca^MnOs samples containing Mn304 below the detection limit with standard powder 
X-ray diffraction 33 . The only cases considered here of little characterized anomalies looking 
like phase transformations are those in YBCO around 240 K (Sec. 15. 7() and 120-160 K (Sec. 
15.5. 1|) . but even in these cases, their peculiar dependence on oxygen doping strongly suggest 
mechanisms involving the CuOx planes of YBCO. 

A convincing indication against the contribution of impurity phases and also against instru- 
mental contributions to the spectra examined in the present study is the strong dependence 
of these spectra on doping and the fact that they are completely different for different fami- 
lies of cuprates. This is shown in Fig. 13-21 comparing the Q^^ (T) curves of ceramic samples 
of La2~xSr2.Cu04, YBa2Cu306+2,' and Nd2-a;Ce2,Cu04+5 at various dopings, measured in the 
same conditions exciting the mode at lower frequency, between 0.8 and 2 kHz. All the samples 
were prepared in the same laboratory by M. Ferretti and coworkers, with the same equipments 
and procedures, except for sintering times and temperatures, but the spectra are completely 
different from each other. The fact that there is no single feature that is shared among all 
the spectra excludes the influence of systematic instrumental effects. The only possible instru- 
mental contribution is from freezing of adsorbed oils around 220 K, as discussed in Sec. 15.71 
but the curve of Ndi.95Ceo.o5Cu04 sets the upper limit of this possible contribution to a value 
too small to affect the spectra of LSCO and YBCO; in addition, a sample of AI2O3 with 50% 
porosity, also shown in Fig. 13-21 does not present any anomaly at that temperature, in spite of 
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Figure 3-2: Elastic energy loss versus tempertaure of ceramic samples of LSCO, YBCO, NCCO 
at various dopings and of porous alumina, measured in the same conditions at frequencies 
between 0.8 and 2 kHz. 

being much more susceptible to oil uptake. The influence of gases adsorbed by the samples is 
discussed in Sec. 15.5.11 The lack of peaks common to all the curves in Fig. 13-21 excludes the 
influence of a spurious phase like copper oxides, certainly absent in the almost flat and very low 
Q^^ (T) curve of Ndi.95Ceo.o5Cu04. The same observation is made by varying doping within 
the same cuprate family, as shown in Fig. 13-31 for LSCO. In this figure, there is a progressive 
change of the spectrum with increasing x, indicating that all the features depend on doping of 
La2-a;Sr2.Cu04 and therefore not on spurious phases. If they contributed to one of the peaks 
labeled T or S or to the maximum at LHe temperature, they should exist also at x = 0.2, where 
instead the elastic energy loss is 2 — 3 orders of magnitude smaller (except for the step at the 
structural transformation from cubic to tetragonal, shifted down to 80 K at such doping). 
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Figure 3-3: Anelastic spectra of La2-a;Srj;Cu04 at various doping levels x in %; the vertical 
coloured arrows indicate the structural transitions from cubic to tetragonal. 

In Fig. 13-31 the only instance in which there is no clear smooth evolution of the curves with 
doping is the double peak at 80 K for x = 0, which should not be identified with peak S in the 
Sr-doped samples (see Sec. I4.11.3() . The curve with x = was measured after a rather strong 
out gassing treatment that possibly introduced O vacancies in the Cu02 planes (see Sec. I4.8.3() . 
and the double peak might be related to such vacancies. The peak at 50 K also appears as some 
doping dependent relaxation in La2-xSra;Cu04, but it is labeled with a question mark since we 
were not able to find a plausible mechanism for it, and it will not be discussed in this Thesis. 

Finally, the possibility that the motion of grain boundaries itself may contribute to the 
present measurements can be excluded; this type of relaxation, in fact, requires the diffusion 
of the cations, which occurs at temperatures comparable to the sintering temperatures, above 
1200 K. Instead, all the thermally activated relaxation processes observed here at high temper- 
ature are clearly related to diffusion of nonstoichiometric oxygen. 
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3.5 The UHV system for sample treatments 



The oxygenation and some of the outgassing treatments were made in a Ultra High Vacuum 
(UHV) system reahzed by Ing. Dalla Bella at RIAL Vacuum (Parma, Italy) after my project. 
It consists of a main spherical chamber with Bayard- Alpert head and Residual Gas Analyzer, 
pumped by turbomolecular, ionic and Ti sublimation pumps. Two CF63 flanges on opposite 
sides connect on one side, through an all-metal gate valve, another chamber connected with 
the gas inlet line, capacitive heads and a quartz tube flanged CF63 where the sample is put 
for the treatments. On the other side, the spherical chamber is connected to a chamber for 
introducing the sample without breaking vacuum in the rest of the system, equipped with a 
magnetically coupled rotary-linear feedthrough with a tray for depositing the sample. The trail 
has the possibility of some lateral movement and is surrounded by stainless steel wire, so that 
with a rotary movement it is possible to retrieve also fragile samples, generally wrapped with 
Pt wire for protection and for avoiding direct contact with the quartz tube. The sample is 
heated up to 1100 °C by a horizontal tubular furnace mounted on wheels. 

The gas inlet line consists of a multiway valve connecting to various pure gas bottles (O2, 
H2, D2, N2), and a bakeablc all metal section with inlet and outlet needle valves, a small 
(17.3 cm^) and a large (200 cm^) calibrated volumes and capacitance head, in order to be able 
to admit a known amount of gas. 

The base vacuum before any treatment was in the 10~^ mbar range, but could be improved 
with baking to the 10"^*^ mbar range in particular cases. 
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Chapter 4 



LSCO 

La2-.x(Sr/Ba)2.Cu04 (LSCO or LBCO) is the first high-Tc superconductor discovered by Bed- 
norz and Miiller. Having Tq < 40 K, it is not particularly attractive for applications, but has the 
simplest structure among the superconducting cuprates and is probably the best characterized. 
Doping may be achieved both through excess oxygen in La2Cu04+5 and by partial substitution 
of La^^ with Sr^+ or Ba^+; in the latter case, if excess oxygen is completely removed, one 
does not have the complications due to the ordering of the nonstoichiometric oxygen, which 
characterize the other HTS cuprates. 

4.1 Structural phase diagram 

LSCO is formed by layers of CuOg octahedra intercalated by (La/Sr) atoms, as shown in Fig. 
14-11 The substitution x trivalent La ions with divalent Sr or Ba introduces p = x holes/unit 
cell in the Cu02 layers, which form a conducting band. Interstitial oxygen (Oj) can also be 
introduced in tetrahedral coordination with four apical O atoms and four La atoms; each 0^~ 
provides two conducting holes. A small amount of excess O (5 ~ 10^'^) is present in as- 
prepared La2Cu04+5 and can be increased to few percent by equilibrating in O2 at moderate 
temperatures •31 ; up to 12% Oj can be introduced by electrochemical oxidation jHSlEni- By 
increasing the Sr doping, the equilibrium concentration of Oi decreases, and electrochemical 
oxidation is necessary to obtain La2-aSr^Cu04+5. 
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Figure 4-1: Structure of LSCO. Left panel: ideal HTT structure; the O atoms (not shown) 
are at the vertices of the octahedra; also shown is an interstitial O atom and its tetrahedric 
coordination with the nearest neighbor apical O atoms. Right paneh the two more stable 
tilt patterns of the octahedra, giving rise to the LTO and LTT structures; in the latter case, 
adjacent layers of the octahedra are rotated by 90°, so that the overah structure is tetragonal. 

4.1.1 Tolerance factor and Low- Temperature Orthorhombic (LTO) phase 

When decreasing temperature, the equilibrium bond lengths in the (La/Sr)0 layers decrease 
faster than the equilibrium CuO bond length within the Cu02 planes, and the resulting lattice 
mismatch is relieved by a buckling of the Cu02 planes jHZl below a temperature Tj(x). This 
phenomenon is typical of the perovskite structure, formed by a three-dimensional network of 
BOg octahedra intercalated by A cations. In the ideal cubic case, which is the high temperature 
structure of most perovskites, the ratio of the A-0 and B-0 bond lengths is /a-o /^b-o = V^, 
and it is usual to define a tolerance or Goldschmidt factor 



which is 1 in the ideal cubic case, and, when becomes < 1, indicates that the size of the 
BOe octahedra is too large for the equilibrium A-0 bond length, resulting in rotations of 
the octahedra. In many cases, the transitions of perovskites from cubic to lower symmetries 



t = 



^A-O 



(4.1) 
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Figure 4-2: Perovskite structure with evidenced the A-0 and B-0 bonds, whose lengths in the 
ideal cubic structure are in the ratio \/2. 



involving rotations of the octahedra are understood in terms of t decreasing with temperature. 

The same reasoning can be applied to the perovskite layers of La2-xSr2,.Cu04, defining the 
tolerance factor |37j 



Since the octahedra are relatively rigid units, the buckling results in a collective tilting of 
the octahedra, and the structure changes from high-temperature tetragonal (HTT) to low- 
temperature orthorhombic (LTO, see Figs. 14-11 and l4-4b '). 

Doping reduces the mismatch between LaO blocks and Cu02 planes in two ways i) substi- 
tution of La^"^ with larger Sr'^^ or insertion of interstitial O expand the lattice, and therefore 
relieve the compressive stress on the Cu02 planes; ii) doping holes in the Cu02 planes removes 
charge from the CuO antibonds, therefore shortening them j35U39j . Therefore, the HTT phase 
is stabilized by doping, and Tt (x) is an almost linearly decreasing function of x, as appears 
from the LSCO phase diagram in Fig. 14-31 In the same figure is also reported the average 
tilt angle 9 in the LTO phase from diffraction measurements, again a decreasing function of 
doping; the decrease seems to be much more regular at a local level, as found by extracting the 
pair-distribution functions (PDF) from neutron diffraction jHEl, and prosecutes into the HTT 
phase, which therefore should consist of disordered tilted instead of untilted octahedra. 
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Figure 4-3: (a) Structural (blue), magnetic (gray) and electric (green) phase diagram of LSCO. 
AF = antiferoomagnetic, SG = spin glass, CSG = cluster spin glass, SC = superconducting, 
(b) average and local tilt angle of LTO phase versus doping (Ref. (3B) ) . 

4.1.2 Other tilt patterns and the Low-Temperature tetragonal (LTT) phase 

Actually, other tilt patterns besides the LTO one are possible, all describable in terms of a 
rotation axis of the octahedra within the ab plane, and therefore in terms of two rotation angles 
Qi and Q2 about two orthogonal axes within ab. Such axes are generally chosen at 45° with 
the Cu-0 bond directions (the direction of the a axes of the LTO cell, while the a and b axes 
of the HTT cell is parallel to the Cu-0 bonds) HHIEI], as shown in Fig. H^ . 

The two variants of the LTO phase are then described by (Qi,0) and (0,(52)) while the 
LTT phase by \Qi\ = \Q2\- The intermediate cases Qi / Q2 7^ are also possible, and produce 
intermediate phases, generally precursors to the LTT one j4HI43| lll]. Tilt patterns intermediate 
between LTT and LTO are present also within the twin boundaries in the LTO phase. Such 
twins walls have been observed to be nucleation sites for the LTT phase |45| . 

The LTT structure is stable only at low temperature, near the doping x = 1/8, and if there 
is sufficient disorder in the ionic sizes in the La sublattice. The latter is obtained by substituting 
La with Ba instead of Sr, which has a still larger radius JISI (the radii of 12-fold coordinated 
Ba^^, Sr^^ and La'^^ are 1.61, 1,44 and 1.36 A), or by doping with Sr^^ and substituting part 
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Figure 4-4: (a) The possible tilts of the CuOg octahedra in LSCO; (b) calculated local potential 
for an octahedron in LBCO, with eight minima corresponding to LTT and LTO tilts (Ref. |42j ) . 

of La'^'^ with the larger Nd^+. 
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Figure 4-5: Low-temperature region of the phase diagram of LSCO and LBCO, with the differ- 
ent regions of stability of the LTT structure below Tf. The open circles are Tt from Ref. |4Uj . 
while the closed one from the measurement of Fig. 14-101 . 



Figure shows the region of stability of the LTT phase for La2-x-/\^xCu04 with M = Sr, 
Ba, deduced from various types of measurements [4Ul I47j . It appears also that the LTT phase 
is associated with a depression of the superconducting Tc. 

It is important to note that in the LTO phase all the O atoms in the Cu02 planes are 
equivalent, while in the LTT phase one can distinguish between the O atoms on the rotation 
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axes, and therefore remaining on the plane, and those shifted out of the plane; this provides 
a modulation of the potential felt by the charge carriers, which results in a mutual 
stabilization of the LTT modulation and of the static hole stripes (see Sec. I4.4j) . whose spacing 
becomes commensurate with the lattice spacing at x = |. 

The local potential felt by each octahedron free to tilt has been theoretically found to 
have eight minima in correspondence of the four possible LTO and four LTT tilts separated 
by barriers E/kB ~ 300 - 500 K EH] (see Fig. ^^). 

4.2 Electric phase diagram 

The electric and magnetic phase diagram of LSCO can be considered as representative of the 
other cuprates, except for complications arising from the O nonstoichiometry in the latter. 

As anticipated above, the concentration of holes doped in the Cu02 planes of La2-a;Sr2:Cu04_t-5 
is p = X + 25, neglecting clustering of Oi (see Sec. 14. 8|) . The transport properties (conductivity, 
Hall coefficient, dielectric permittivity) of lightly doped La2--xSr2:Cu04+5 can be understood in 
terms of conventional semiconductor physics [491 1501 1ST] : La2Cu04 has static dielectric constant 
e = 310, and the hole effective mass is ruh = 2me. The holes are thermally ionized from the 
acceptors, to which are bound with an energy i?t(Sr) = 10 — 20 meV for the case of Sr dopants 
and Ei,{0) ~ 31 meV for Oi. Conduction is of band-type at high temperature, namely a oc p (x 
exp{—Eb/kT) at T > 70 K, and variable-range-hopping below 50 K (cr oc exp[— (Tq/T)^/'^]). 

For p > 0.05 the planes start to superconduct below Tc {p), which has a maximum versus 
doping at popt — 0.15 (see Figs. 14-31 and l4-5() : the system is called underdoped, optimally doped 
and overdoped, depending on the value of p with respect to Popt- Underdoped cuprates exhibit 
various anomalies, partly interpreted in terms of opening of pseudogaps in the charge or spin 
excitations and partly in terms of charge stripes; the latter will be dealt with in some detail in 
Sect 14.41 In overdoped cuprates, instead, a uniform metallic state sets in and becomes so stable 
that superconductivity eventually disappears. 

In cuprates of the LSCO family, a depression of superconductivity occurs in correspondence 
with the formation of the LTT phase, and this is understood in terms of locking of the charge 
stripes to the LTT lattice modulation, as explained in Sect 14.41 
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4.3 Magnetic phase diagram 



The study of the superconducting and magnetic phase diagram of the Cu02 planes of the 
superconducting cuprates is a complex and fascinating subject (for a review see e.g. JZl); but 
I will mention only the issues relevant to the anelastic measurements. In the absence of doping, 
the Cu02 planes are semiconducting, with Cu in the Cu^'^ oxidation state having spin s = ^; 
these spins order antiferromagnetically (AF) below the Neel temperature Tn = 315 K, with 
the staggered magnetization within the ab plane, and mainly parallel to b. There is also a small 
component of the spins pointing out of the planes, due to the small tilt of the Cu-0 bases of 
the octahedra; this weak canting produces a ferromagnetic component that dominates the low 
frequency magnetic susceptibility [52] • Doping holes causes some Cu atoms to pass into the 
Cu^^ state with s = 0, and this disturbs the AF order; in La2-xSr^Cu04 Tn (x) drops to K 
already at the critical doping Xc ~ 0.02 (see Fig. lA-'A^ . 

A model of how the holes disturb the AF order has been proposed by Gooding et al. |53j . 
starting from the hypothesis that at low temperature and low doping the holes are localized 
near the Sr dopants. The ground state for one isolated hole would be doubly degenerate with 
the hole circulating either clockwise or anticlockwise over the four Cu atoms nearest neighbors 
to the Sr atom. The hole motion couples to the transverse fluctuations of the Cu spins and 
produces a spiraling distortion of the AF order within the ab plane; at this point it should be 
noted that the model assumes the staggered magnetization of the hole-free plane along c, while 
in fact it is along b, but the model may help in focusing some mechanisms responsible for the 
magnetic phase diagram of the HTS cuprates. The ground states with disordered distributions 
of Sr impurities would consist of AF correlated domains delimited by the Sr atoms and with the 
in-plane AF order parameter randomly oriented, resulting in a (cluster) spin-glass state. The 
domains would be separated by narrow domain walls with disordered spins and ferromagnetic 
character which connect the Sr atoms. The hole mobility is much higher in FM rather than 
AFM regions (the hopping of a hole in a AF domain requires also a spin flip, if destruction of 
AF order has to be avoided), and therefore with increasing temperature and doping the holes 
move along these domain walls, which would therefore correspond to the charge stripes of 
the next Section. 

Whatever model is chosen, experiments probing the local spin fluctuations, like NQR [TTHSU] 
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and ^SR [HI], indicate that the spin degrees of freedom associated with the doped holes are 
different from the in-plane Cu^"*" spin degrees of freedom that order themselves below Tn, and 
the localization of the doped holes allows the associated spins to progressively slow down and 
freeze |171 l5Uj . For x < Xc one has long range AF order below Tn, and the doped spins freeze 
into a spin glass (SG) state below Tj (x) linearly increasing with x (see Fig. l4-c{|) . For 
X > Xc there is no long range AF order and approaching Tg ~ 0.2 K/x AF correlations develop 
within domains separated by hole-rich walls, and with easy axes uncorrelated between different 
clusters, giving rise to a cluster spin-glass (CSG) state (see Fig. 1 4-3(1 . 

More recent neutron scattering experiments [HH] of the magnetic correlations in La2_xSr2,.Cu04 
for X < Xc suggests a different picture of the spin glass phase, with the 3D AF ordered phase 
coexisting below ~ 30 K with domains of the stripe phase observed for x > Xc (see also next 
Section). It has been proposed that the hole localization starting around 150 K involves an 
electronic phase separation into regions with xi ~ and X2 ~ 0.02, and the volume fraction 
of the X2 = 0.02 phase changes as a function of the Sr doping, in order to achieve the average 

X. 

4A Hole stripes 

A phenomenon that seems to be common to many superconducting cuprates, and has been 
attracting enormous interest, is the segregation of the conducting holes into stripes, 
while maintaining very good conductivity or even superconductivity [Ml EZl EHl EHI ^1 12] . 
The literature on the subject is vast, and I will deal only with those aspects related to the 
anelasticity. On the theoretical side, it is debated whether these charge stripes compete against 
superconductivity ;3^ or on the contrary they are an essential ingredient of HTS In both 

cases, an important issue is the dynamics of the transverse fluctuations of these stripes, which 
has been modeled in terms of collective pinning from the doping impurities by C. Morals Smith 
|60j . To my knowledge, the only experimental results on the low frequency stripe fluctuations 
are the anelastic measurements presented here IHH ESI IMl ■ 

The first indications of segregation of the charge carriers into stripes came from measure- 
ments of the correlation length ^ (x) for the AF order in LSCO, deduced from Tn (x) jSfij and 
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Figure 4-6: Left: ordered stripes with periodicity d commensurate with the lattice at x = |. 
The holes (Cu^"*" with spin s = 0) are in red, while the violet sites along the hole stripes do not 
have a definite spin direction, since they are also walls between different AF domains. Right: 
possible situation if the holes were uniformly distributed. 

from NQR experiments the observation is that the increase of on cooling is limited to 
a length / ~ a/x (a is the lattice constant). Since the AF correlations can develop only over 
regions free holes (which have s = instead of / should represent the size of domains free of 
holes and, if the holes were uniformly distributed over the Cu02 plane, their separation should 
scale as / ~ a/^/x, while I ^ a/x indicates that the holes are in one-dimensional walls of fixed 
hole density separating hole- free domains or stripes (see also Fig. I4-6|) . Note that these charge 
stripes are not charge-density waves, having a much sharper modulation and allowing the Cu^^ 
spins to form AF domains between in the charge-poor regions. 

Several other indirect indications of the existence of the charge stripes have been found 
jl7l I67j . including the observation of inhomogeneous Cu-0 bond lengths in the Cu02 planes 
with probes of the local structure like EXAFS jl| and pair-distribution functions (PDF) from 
neutron diffraction [HS] . The strongest evidence of the existence of parallel magnetic and charge 
stripes comes from magnetic inelastic neutron diffraction, which reveals one-dimensional dy- 
namic charge and magnetic correlations with a spacing d incommensurate with the lattice 
parameter and decreasing with doping as d oc 1/x |69j . In addition, the direction of the mod- 
ulation changes from diagonal to parallel with respect to the Cu-0 bonds when increasing 
doping above x = 0.055, which also separates the semiconducting from the superconducting 
region |591 170j . These correlations are observable also as static structural modulations in the 
region of the phase diagram x ~ g with the LTT phase stabilized by partial substitution of La 
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with Nd [SniETj. In fact, at x = ^ the stripe spacing becomes commensurate with the lattice, 
and the LTT tilt pattern provides a modulation to which the stripes are locked. The situation 
for X = ^ deduced from neutron diffraction is represented in the left hand of Fig. 14-61 The hole 
stripes act as antiphase boundaries between regions where the Cu^+ spins have AF correlation; 
varying doping does not modify the hole density within a charge stripe, which remains 0.5, but 
only the stripe separation (hence d oc 1/x). It is also found that, on cooling, the formation of 
the charge stripes precedes that of the AF spin stripes (56j . 

4.5 Nd2-a;Ce^Cu04+<5 

The structure of La2Cu04, also called T structure [JJ, is one of three possible structures of 
A2BO4. Nd2-xCe2;Cu04, instead, presents the so-called T' structure, where the cations and 
the O atoms of the Cu02 planes are in the same positions as in the T structures, while the O 
atoms in the (Nd/Ce)0 blocks correspond to the interstitial sites of La2Cu04+5; on the other 
hand, the interstitial positions in the T' structure correspond to the apical O atoms in the T 
structure. The correspondence between the oxygen positions in the two structures is shown by 
the arrows in Fig. 14-71 The result of this shift in the oxygen positions is that there are no short 
Cu-0 bonds along the c axis and therefore no CuOg octahedra, but only flat Cu02 planes. 

It should be mentioned that Nd2Cu04 supports only electron doping (substituting Nd^"*" 
with Ce^+), at variance with all the other superconducting cuprates. 

4.6 The anelastic spectrum 

The anelastic spectrum of LSCO contains several relaxation processes whose intensity and ap- 
pearance strongly depend on the type and level of doping. Figure present the elastic energy 
loss and Young's modulus between 1 and 800 K of stoichiometric semiconducting La2Cu04 and 
Nd2Cu04. It should be stressed that as-prepared La2Cu04+5 has 6 ~ 0.003 that drastically 
modifies the anelastic spectrum, and the result of Fig. 14-81 is obtained after accurate outgassing 
of the excess oxygen. Nd2Cu04 appears like a normal solid without defects or excitations: the 
absorption is low and the elastic modulus decreases regularly by less than 20% between and 
800 K. Stoichiometric La2Cu04 is also free of defects, in principle, but its anelastic spectrum 



54 




Figure 4-7: Comparison between the T (LSCO) and T ' (NCCO) structures. The arrows 
indicate the correspondence between the O interstitial positions (white) in one structure and 
O in the NdO or LaO blocks in the other one. 

presents extremely intense anomalies. The main difference between the two compounds is that 
Nd2Cu04 has flat Cu02 planes while La2Cu04 has CuOg octahedra unstable against tilting 
(see Sec. 14.11 and Fig. I4-7|) . In fact, almost all the anelastic processes in La2Cu04 are due to 
some type of motion of the octahedra. Starting from high temperature we find (the numbers are 
those in Fig. 14- 7() : 1) the transformation from HTT to LTO structure; 2) motion of the domain 
walls between the two LTO variants (orthorhombic b axis along x or y); 3) octahedra tilt waves 
of solitonic type with thermally activated dynamics; 4) local tilts with dynamics governed by 
tunneling and interaction with the charge excitations; 5) thermally activated fluctuations of the 
charge stripes whose interaction with the lattice is mainly mediated by the octahedral tilts. 

Apart from the well known structural HTT/LTO transformation, all the other processes 
have been revealed by the anelastic experiments presented in the next Sections. 
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Figure 4-8: Elastic energy loss and Young's modulus of La2Cu04 (500 Hz) and Nd2Cu04 (800 
Hz) from 1 to 800 K. 

4.7 Structural phase transitions 

4.7.1 HTT/LTO transformation and determination of the Sr content from 

Tt 

The acoustic anomalies connected with this transformation have been studied by other authors 
|72| I73j and, in view of the polycrystalline nature of our samples, I did not attempt any quan- 
titative analysis of the huge elastic anomaly and the accompanying rise in acoustic losses (#1 
in Fig. l4-8j) . I only observe, together with Lee, Lew and Nowick [221, that only the region 
at T > Tt may be analyzed in terms of Landau free energy, as in Ref. [ZSl, because below 
Tt the motion of the walls between the two possible LTO variants is predominant both in the 
imaginary and real moduli. The motion of the walls is very sensitive to defects like Oi and 
Sr, and is difficult to be modeled; their effect is to mask in most cases the cusps or kinks that 
would otherwise be observed in the moduli, making difficult even to establish where Tt exactly 
is. 

In the present Thesis the main interest in analyzing the HTT/LTO transition is for an 
accurate determination of the actual concentration of Sr and an estimate of its homogeneity. 
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In fact, the transition temperature Tt is strongly dependent on doping [T^ 

rt(x) = (1 -x/0.217) X 535 K , (4.3) 

and the width of the modulus step provides an upper limit to the possible inhomogeneities in 
X over a sample, or, at least, by comparing widths of different E (T) curves it is possible to tell 
whether a particular sample appears to be more inhomogeneous than the others. 




Figure 4-9: (a) Derivative of the relative variation of the Young's modulus of La2-xSra;Cu04 at 
the HTT/LTO transition at various values of x. (b) Sr content deduced from Tt versus nominal 

X. 

As mentioned above, there is no model available for fitting the (T) and E (T) curve; 
as a uniform criterion to extract Tt from the anelastic spectra, I chose to identify Tt with 
the inflection point of E[T). In practice I fitted the derivative of the relative variation of 
the Young's modulus, d\n[E / E [{))]/ dT , with lorentzian peaks, as in Fig. l4-9b . obtaining 
Tt from the temperature of the maximum, and a transition width AT^ from the peak width. 
The method seems to work, since the curve for x = has a maximum exactly at 535 K and 
the fitting Lorentzian at 536 K, as expected from Eq. (|4.3jl (it should be mentioned that 
sometimes distorted shapes of E (T) have been measured on La2Cu04+5, possibly connected 
with the presence of Oi). The Sr concentrations estimated in this manner for the samples 
investigated here are plotted against the nominal x in Fig. l4-9b . and show a good agreement. 
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The relationship between transition width and inhomogeneity of the Sr concentration cannot 
be simply deduced by translating the peak width AT into Ax, since there is a considerable 
intrinsic width of the transition (critical softening above Tt and domain wall motion below). 
This is particularly evident by observing that the curve for x = 0, where Ax = 0, is about twice 
broader than those for x ~ 0.015. It is seen, however, that transition ^2 is broader than the 
others at similar doping, suggesting homogeneity problems with that sample. 

4.7.2 LTO/LTT transformation 

As anticipated in Sec. I4.1.2( extended domains of LTT phase can be found only in LBCO 
or Nd-substituted LSCO near x = |. The signature of the LTO/LTT transformation in the 
anelastic spectrum has been identified as a small acoustic absorption peak and stiffening below 
the transition temperature Tt 47 . Figure H- lUI presents the Young's modulus curves normalized 
to the value extrapolated to T = K for samples where La is substituted with 3% and 6% Sr 
and Ba. 
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Figure 4-10: Young's modulus of LSCO and LBCO at x = 0.03 and 0.06. Tg indicates the CSC 
freezing temperature, and Tt the onset of the LTO/LTT transition. 
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The only clear indication of formation of LTT phase is for the 6% Ba sample, below Tt 
~ 40 K; such a stiffening is compatible with the ultrasonic measurements on LBCO with 
0.10 < X < 0.16 [321 and with the available phase diagram |401 175j . as shown in Fig. 14-51 

Also LSCO at X = 0.06 exhibits similar weak anomalies, which suggest a tendency to 
the formation of LTT domains also in LSCO, in accordance with high-resolution diffraction 
experiments jZS,- In LSCO, however, such domains should be either confined to the twin 
boundaries or fluctuating, extremely small, and without long range correlation; instead, in 
LBCO and LSCO co-doped with Nd the diffraction experiments reveal a stable phase with 
long range order, which can also provide a pinning potential for the stripes. The anelastic 
experiments do not provide any direct information on the extension or topology of the LTT 
domains, and the elastic anomalies are likely connected with the domain boundaries. It is 
therefore possible that narrow domains of minority LTT phase in LSCO, with a high perimeter 
to area ratio, and extended LTT domains in LBCO produce elastic anomalies of comparable 
amplitude. 

4.8 Interstitial oxygen 

In literature extensive reports can be found on La2Cu04+5 that is heavily doped by electrochem- 
ical methods, generally enough to become superconducting. The interest in these researches 
is on electrical and structural phase separation and various types of ordering of Oi that occur 
under those conditions. My research is instead focused on the low concentration limit. Under 
normal conditions, the amount of excess oxygen in La2Cu04_|_5 is very small: by electrochemical 
reduction it has been estimated 6 ~ 0.005 in the as-prepared state |77| I78j : many estimates of 
6 ~ 0.01 — 0.02 have also been reported, but, in view of the results presented in Sec. 14.8.21 I 
think that 6 < 0.005 is a more reliable evaluation. 

Figure B-llb presents the effect of outgassing La2Cu04+5 at progressively higher temper- 
atures in vacuum starting from the as prepared state. Initially, two peaks labelled 01 and 
02 are present at 225 and 300 K (at 300 Hz). Outgassing in vacuum reduces the intensity of 
these peaks in favor of another peak, labeled T, at 150 K. It is obvious to assign 01 and 02 
to interstitial O, which is known to be present in small quantities in the sites indicated in Fig. 
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14-11 and is lost in vacuum at high temperature jTS]. The hypothesis is confirmed by fohowing 
the O2 partial pressure during heating at 6 K/min in dynamic vacuum in a UHV system j34j : 
PO2 starts increasing above 550 K, indicating the loss of interstitial O, as also explained in Sec. 

The phenomenology is confirmed in Fig. 14-llb with another sample equilibrated at 900 K 
in partial pressures of O2 ranging from 0.1 to 820 torr WT. Each oxygenation was preceded by 
an outgassing at 700-750 °C, in order to obtain an initial reference state; the sample was then 
equilibrated at 620 °C for 90-120 min in a static atmosphere of pure O2 and rapidly cooled. 



n no 



■ir^:f d in ^'a^nnn up Is 




U.1U - 



USA - 



■ 1 H ?Af.) C 




2 Kb] 






. «qiiiihrnhd nrnOfi K i"! / 






■ T // 








- f\ 1 









ma 



3C3 



4J1CI 



Figure 4-11: Dependence of the anelastic spectrum of La2Cu04+5 on 6 (a) outgassing an as 
prepared sample at progressively higher temperatures in vacuum (300 Hz); (b) equilibrating at 
900 K at different O pressures (800 Hz). 

Figure 14-121 presents a fit to the Q^^ (T) curves after equilibration in 10 torr with two 
Fuoss-Kirkwood peaks, Eq. ()2.42|) . excluding the region at the right of the dashed lines, where 
contributions from the structural transition and the DW motion are important. The fitting 
parameters are tq = 2 x 10~^^ s , Eoi/k-Q = 6000 K and a = 0.4 for peak 01 and tq = 
8 X 10-1^ s , Eo2/kB = 7000 K and q = 0.5 for peak 02. 

The most likely explanation for peak 01 and 02 is that Ol is due to hopping of isolated 
Of and 02 to pairs of O;, possibly peroxyde complexes with lower oxidation state, as 
suggested by the dependence of 5 on the oxygen pressure (see Sec. I4.8.2() . This hypothesis 
explains the fact that peak 01, expected to increase linearly with 5, is rapidly overwhelmed by 
02, which should increase with 5^. The parameter of peak 01 is of the order of magnitude 
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Figure 4-12: Anelastic spectrum of La2Cu04+5 in the region of the peaks due to hopping of 
isolated (01) and paired (02) interstitial O atoms. The curves are a fit of the data at the left 
of the dashed lines. 

of the local phonon modes promoting the atomic jumps, and the resulting rate (T) agrees 
with the observation that the excess oxygen becomes immobile below 140-150 K 81 : in that 
temperature range the estimated mean time r between subsequent jumps passes from 5 min 
to more than 1 hour. It should be mentioned that experiments based on the analysis of the 
time and temperature dependence of the superconducting properties, which in turn depend on 
the ordering of Oi [HJ, provide activation enthalpies different from Eqi. Those experiments, 
however, determine the characteristic time of a complex process, like the aggregation of oxygen 
into domains and the possible ordering within the domains, while peak Oldirectly probes the 
hopping rate for single Oi defects. 

Two features of peaks 01 and 02 seem at first at odds with a mechanism based on jumps of 
Oi: the fact that they are very broad and that they shift in temperature with increasing doping. 
The broad width can be justified by the strong interaction between excess oxygen and 
the tilts of the surrounding octahedra, as discussed in the next Section. The network 
of these octahedra with their easy tilting modes provides a disordered environment for excess 
oxygen and mediates the interaction between the interstitial atoms much more effectively than 
in a regular lattice. In fact, at lower values of 6 peak 01 is narrower {e.g. a = 0.7 for the 
intermediate curve of Fig. 14-llb , [80 J. This could also explain the shift of peak 01 to higher 
temperature with increasing 6 in terms of correlated and therefore slower dynamics of Oi. 
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These results can be compared with similar measurements made by Kappesser et al. j82j . 
where only one thermally activated peak was observed, which was attributed to oxygen diffusion. 
Since that peak was found after electrochemical doping to 5 = 0.013 and 5 = 0.035, it should 
correspond to peak 02, with 01 possibly masked as in the highest curve of Fig. 14-llb . The 
peak temperature is close to that of 02 at a similar frequency, but the attempt frequency and 
activation enthalpy were estimated to be rather low: Tq"^ = 6.5 x 10^ s and H/k-Q = 3600 K. 

4.8.1 The structure of the interstitial oxygen defect 

The interstitial sites 0(i) for Oi have been determined by neutron diffraction |83| I84[ I36j 
and are shown in Figs. 14-11 and 14- 13b . In the HTT structure and in the absence of ordering 
into superstructures they are at the center of regular tetrahedra formed by the neighboring 
apical oxygen atoms, 0(a), with 0(a)-0(i) distances in the range 1.6 ^ 2.2 A, and lar ger 
regular tetrahedra formed by La atoms. Inspection of Fig. I4-I3b shows that all the interstitial 
sites for oxygen in the ideal HTT structure are equivalent from the elastic point of view. In 
fact, by symmetry the elastic dipole of Oj in the red site must have principal axes z \\ c and x 
and y within the ab plane, parallel to the horizontal edges of the tetrahedron; the tetrahedron 
corresponding to the neighboring green site is rotated by 90° about the c axis with respect to 
the red one, but it may be also be obtained from it by inversion. This means that also the 
elastic dipole of Oj in the green sites may be obtained by inversion of the red one, and therefore 
the principal axes x and y must be equivalent; then, the strain due to Oi is isotropic within the 
ah plane and there is no change of the elastic dipole for jumps of Oj in the ideal HTT structure. 

In the LTO structure, where the anelastic relaxation peaks are observed, the octahedra 
are tilted about the a axis (of the LTO representation). Figure I4-I3b shows a top view of 
three interstitial sites (red and green circles) and the surrounding apical O atoms, with their 
displacements with respect to the tetragonal structure. The black circles are 0(a) below the 
plane of 0(i), with the corresponding octahedra shaded in gray, while the white circles are 
0(a) above it. The sites are labeled as 1 (red) or 2 (green) depending whether the 0(i)-0(a) 
distances are shorter along (llO) (or x) or along (110) (or y). After the insertion of an O atom, 
the lattice will expand in different ways along the directions of the shorter bonds with respect 
to the longer bonds. In the figure is represented the case in which the 0(a) which are closer 
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Figure 4-13: (a) two interstitial O atoms in the HTT structure; b) interstitial sites for O in the 
LTO structure labeled as 1 or 2 depending whether the shorter distance with the surrounding 
apical O atoms is along (llO) or (110). The white (black) circles are the apical 0(a) atoms 
above (below) the plane of the 0(i) sites; the arrows show their shift with respect to the ideal 
tetragonal HTT structure. The octahedra below 0(i) are shadowed. The elastic dipoles due to 
the interstitial O atom are represented as ellipsoids with major axis along the shorter 0(i)-0(a) 
distances. 



are more displaced outwards (but in the case of the formation of a short peroxide bond the 
opposite can be true): the elastic dipole, represented as an ellipse, has the major axis along 
the direction of the shorter bonds and reorients by 90° when the O atom jumps to a site of 
different type, modulating the Eq strain and giving rise to anelastic relaxation of the cqq elastic 
constant. 

The anisotropy of the elastic dipole in the ab plane arises from the deviation of 
the LTO structure with respect to the HTT one, and therefore is expected to depend 
on temperature and doping. In fact, the tilt angles of the octahedra decrease continuously 
with doping, at least in La2-xSr2.Cu04+5 85_; in addition, the orthorhombicity of La2Cu04+5, 
defined as (6 — a) /a, varies from 0.0165 at 5 = to 0.0119 in the oxygen-rich phase [361 186j . 

From Fig. 14- 13b it also appears that only jumps along a contribute to anelastic relaxation, 
because they involve a rotation by 90° of the anisotropic elastic dipole, while jumps in the b 
direction do not. The hopping rates in the two directions may also be different, as discussed in 
more detail in Ref. |8()j . but up to now there are no experimental indications in this sense. 
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4.8.2 Dependence of S on the O2 pressure at high temperature 

Under normal conditions, the amount of excess oxygen in La2Cu04+5 is very small: by elec- 
trochemical reduction it has been estimated 6 ~ 0.005 in the as-prepared state [7711751, but no 
systematic data exist of the equilibrium values of 5 as a function of temperature and pressure 
(the p6T phase diagram). For this reason, I studied the absorption of interstitial O and the 
formation of O vacancies in La2Cu04-|-(5 in equilibrium with O2 pressures between 10"^ and 10'^ 
torr or in vacuum at high temperature, in the UHV system described in Sec. IH.51 ,34^. Such a 
study could be done in a more straightforward manner with a thermobalance; therefore I will 
omit the complications arising from measuring the amount of O absorbed or lost in a closed 
volume [351 ■will concentrate on the results and the implications on the aggregation and 
valence state of excess oxygen in La2Cu04+5. 

Two different temperature regions are found: below about 600 °C the equilibrium down 
to very low pressures of O2 involves only interstitial oxygen atoms with negligible influence 
of O vacancies (Vq), while above 700 °C large amount of vacancies are created. 

I determined the equilibrium values of 5 as a function of pressure at 550 °C, and in view of 
the very low amount of oxygen absorbed {5 < 0.01), it seems appropriate to initially suppose 
that Oi dissolves as O^" ions, with the charge balance guaranteed by the change of the valence 
of Cu from Cu'^+ to Cu^"'' (La'^"^ does not support mixed valence). The reaction involved in the 
equilibrium between gaseous O2 and interstitial can then be written as 

^02 + 2Cu2+ < — > of- + 2Cu^+ (4.4) 

and at equilibrium the concentrations must satisfy the mass action law 

roM rcu3+i' ^ ^ 

K = L^lJJ= J 4.5 

where K (T) is the reaction constant and the square brackets represent the molar concentrations, 
assumed to be <C 1. Introducing the neutrality condition, [Cu^^] = 2 [Of^] = 25, and the fact 
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that for (5 <C 1 it is [Cu^^] ~ 1, one obtains 

(5ocpi/^ (4.6) 




Figure 4-14: The experimental 5 after equihbration at 550 °C in a pressure p of O2 is within 
the gray region; the continuous hne is a fit with Eq. ( I4.1UI ) and a = 0.003, b = 0.02. 

The experimental S (p) relationship is the gray region in Fig. 14- 141 and shows deviation from 
6 oc p^^^ already at 5 ~ 0.002. A transition to p^^"' laws with n = 4, 3, 2 could result from the 
formation of peroxide species like Oj~ or neutral interstitial O; for example, the formation of 
neutral Oj involves the reaction < — >0i with the equilibrium condition K = S/p^^"^ instead 

of Eq. mni 

Peroxide species are indeed found in La2Cu04+5 [HSIIHII; but at much higher doping levels; 
on the other hand, the anelastic spectra of Fig. 14-lll clearlv show that interstitial O forms two 
distinct species already at 6 <^ 0.01; there are also indications that the second species consists 
of pairs of interstitial O atoms stabilized by partially covalent bonds with one of the neighboring 
apical O atoms. These peroxide complexes can be indicated as 02j~ in the equations for the 
chemical equilibrium, since involve two interstitial O atoms, and the reaction for their formation 
is 

O2 + 2Cu2+ < — > O^r + 2Cu3+. (4.7) 
65 



We introduce the notation ci = [O^ ] , C2 = 2 [Ogj ], h = [Cu^"*"] , which yield 6 = ci + C2 and 
the neutrahty condition 

h = 2ci+C2. (4.8) 

The equihbria between gaseous O2 and (Eq. I4.5jl and between gaseous O2 and Oaf can 
be written as 

fe/2)tf ^ 



which yield 



1/6 



(2 + ^pi/2) 



2/3 



2i^' 1/2 

C2 = -J^P ' Ci 



and 



1/6. (l + ft^'^/') 



<5 (p) = ci + C2 = apV^-^ (4.10) 

with a = K^/^and b = At low pressure, O2 dissolves only as isolated O- ions, with 

6 ~ , while at p > b^'^ also the peroxide pairs are formed and 6 tends to {IK')"^^^ p^^^ . 

The case without formation of the peroxide pairs is reobtained setting 6 = (or K' = 0). The 
continuous line in Fig. ITTll is Eq. (|I?TU)) with a = 3.0 x 10^^ torr^^/^, b = 0.02^^ torr"^/^, 
which yield a perfect fit except between 0.2 and 0.4 torr, where it is higher than the experiment. 
It can be concluded that around 550 °C the equilibrium content of excess oxygen in La2Cu04+5 
depends on the O2 pressure as expected from the dissolution as interstitial divalent O ions 
{5 ~ p^^^), but with the formation of peroxide complexes already at (5 ~ 2 x 10 
the concentration of such complexes presumably increases with decreasing temperature, in 
agreement with the anelasticity results. 
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4.8.3 Coexistence of interstitial and vacancy defects 

When heating La2Cu04_|_5 above 700 °C in a closed volume, the evolution of O2 becomes 
massive, indicating that O vacancies (Vq) start forming (up to 6 = —0.08 for T = 750 °C). 
Such vacancies can be reversibly filled by increasing , and no sample decomposition occurs: 
after 9 cycles of outgassing followed by oxygenation, a sample was still superconducting at 
34 K, with a resistivity at 5 K lower than 1 /iOhm cm, and the anelastic spectrum was still 
reproducible [HI] . 

The kinetics for the formation and filling of Vq has been measured by cycling in temper- 
ature and pressure and recording the p{t) curves jMj i and is much slower than that of Oi. If 
reoxygenation is done at the relatively low temperature of 200 °C, a small and fast decrease of 
p (t) in a closed volume indicates that 6 ~ 0.002 oxygen enters as Oi, but the filling of the Vq 
requires equilibration times of years [H^. The first consequence is that Oi and Vq may coexist; 
the second is that they must be in the Cu02 planes, otherwise they would be immediately 
filled by the much more mobile Oi. The coexistence of Vq and Oi is confirmed by the fact that 
traces of the relaxation peaks due to hopping of Oi are seen both in the anelastic and NQR 
spectra [7^1 IHH| of samples which have been outgassed above 700 °C and cooled in vacuum; 
such samples should have up to 8% of O vacancies in the Cu02 plane. This would also be in 
accordance with the observation by neutron diffraction that vacancies are formed in the Cu02 
plane and not at the apical positions jg?j . 

It can be concluded that if one tries to obtain perfectly stoichiometric La2Cu04 with a 
long anneal at temperatures exceeding 700 °C in flowing inert gas, then a large amount of O 
vacancies is introduced in the Cu02 plane, unless the inert gas contains O2 as an impurity, with 
a partial pressure larger than 0.1 torr. 

The absorption kinetics at room temperature could not be followed, but it seems that 
interstitial oxygen can enter in a nearly stoichiometric sample even at room tem- 
perature. In fact, a sample that was carefully reduced in vacuum and II2 atmosphere displayed 
intense peaks due to interstitial oxygen hopping after the sample had been kept for 10 months 
in a glass tube closed with a rubber bung and containing silica gel for absorbing humidity. 
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4.9 Collective tilt motion of the oxygen octahedra 

In this section I will deal with the peak around 150 K labeled as (3) in Fig. 14-81 and T in 
Fig. 14-111 This relaxation process is thermally activated, as shown in Fig. I4-I5h .. and has 
the particularity of presenting the maximum intensity in the perfectly stoichiometric 
material. This fact is demonstrated in Fig. 14- 15b . showing the effect of Oj and substitu- 
tional Sr. The relaxation must be assigned to some intrinsic lattice mechanism, in view of the 
extremely high intensity that it may reach in a stoichiometric and defect-free sample; the am- 
plitude of the process may be better appreciated by noting that the maximum Q-i is only 0.08 
in Fig. 14- 15b . because of the peak broadening, but the corresponding relative modulus defect is 
/S.E / E = 0.25. An intrinsic thermally activated relaxation in La2Cu04 is certainly connected 
with the unstable tilt modes of the octahedra (see Sec. 14. ip 79^, but it must also be 
different from the motion of twin walls between LTO variants 72 , which should be responsible 
for the broad relaxation maximum below the structural transformation at Tt , marked "DW" 
in Fig. 14-121 and "2" in Fig. 14-81 It must also be a collective motion of the octahedra, 
since it is completely blocked hy 5 < 0.001 of Oi: peak T in the curve (5 = 8 x 10"^ oflHSb) 
is completely absent (the rise above 150 K is due to peak 01 from Oi hopping). Substitu- 
tional Sr is instead much less effective than Oi in blocking this collective tilt motion, since peak 
T is reduced of only ^ by x = 0.019, a concentration 20 times larger than the concentration 
of Oi that completely inhibits the peak; the evolution of peak T can be followed for x up to 
at least 0.08 (see Fig. I4-23|) . This fact can be easily explained by looking at Fig. 14-11 or H'-13b 
and considering that an 0^~ pushes outwards 4 apical O atoms and therefore firmly blocks the 
tilts of the corresponding octahedra, which in turn may block several surrounding octahedra. 
Instead, substitutional Sr^+ constitutes a much weaker disturbance, than interstitial O^^, both 
from the steric point of view and because it brings only one instead of two holes. 

A possible picture of the excitations involving the planes of CuOg octahedra has been 
provided by Markiewicz j90j . In a model of rigid octahedra free to tilt as in the LTO structure 
(Fig. I4-1|1 . the tilt of an octahedron propagates all over the plane (Fig. I4-I6fe i'). but if the 
rotation axes pass through the CuO bonds as in the LTT case, the correlations are only in rows 
of octahedra along directions perpendicular to the tilt axes, since they share the O atoms which 
are displaced from the CuO plane, while adjacent rows are little correlated, since the share O 
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Figure 4-15: (a) peak T attributed to the tilt waves of the octahedra in stoichiometric La2Cu04 
measured at three frequencies, (b) Blocking effect of 8 x 10"'* excess O and 0.019 substitutional 
Sr on peak T and on the relaxation at LHe temperatures. 

atoms remaining in the CuO plane (Fig. I4-I6b '). Markiewicz therefore observed that a system 
with LTT ground state is reduced to a one-dimensional array of rows of octahedra, 

where all the octahedra belonging to the same row are tilted according to the same pattern. The 
resulting equations of motion are non-linear and admit solitonic solutions, which correspond 
to one-dimensional propagating walls between domains with different tilt patterns 
j90j . as schematically represented in Fig. I4-I6b ). Although extended X-ray absorption fine 
structure (EXAFS) and atomic pair distribution function [HHl (PDF) measurements of 
La2_xSra;Cu04 exclude a prevalent LTT local tilt at small x, the model by Markiewicz 90 of a 
LTO structure arising from a LTT ground state (dynamic Jahn- Teller phase) seems to provide a 
good framework for analyzing the dynamics of the connected octahedra without charge doping. 

Additional evidence that peak T is connected with the oxygen octahedra and more insight 
on the nature of motion involved came from the collaboration with Rigamonti and Corti, who 
made Nuclear Quadrupolar Resonance (NQR) experiments on the same samples measured 
by us '91 . Details on the NQR principles and technique can be found in Refs. |9Hll7j and will 
be omitted here for reasons of space. I will only mention that in those experiments the spins of 
the ^^^La nuclei can be driven out of equilibrium by radio frequency pulses, and the effective 
relaxation rate Wq for reaching equilibrium can be measured. Such a rate is determined in 
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Figure 4-16: Tilt patterns of the LTO (a) and LTT (b) phases; (c) possible representation of tilt 
waves, where only the apical atoms (filled circles) and their average positions (empty circles) 
are represented. 

the present case by the fluctuations of electric field gradient (EFG) at the La nuclei, and in 
particular by the fluctuating displacements of the surrounding apical O atoms. Under 
such conditions, there is correspondence between Wq {uj,T) and Q^^ {lu,T), since (Hj 



where J is the spectral density of the displacements x of the apical O atoms. On the other 
hand, as shown in Sec. 12.21 also i^iT) oc J, since the same displacements, for nearly rigid 
octahedra, are coupled to the in-plane shear. Therefore, if such displacements produce peak T 
in the elastic energy loss, they must also produce a peak in Wq (w, T) . This has been indeed 
observed [2], as shown in Fig. 14-171 

The relaxation rate t~^, determined from the condition ujt = 1 at the maxima of the Q^^ (T) 
and Wq (T) curves, is also plotted in logarithmic scale against the reciprocal of temperature. 
The three points can be closely fitted by a straight line = TQ^exp{—Ecs/kBT), tq = 
1.9 X 10~^^ s and Ecs/kB = 2800 K, clearly indicating that the same process is observed by 
both techniques. Both sets of curves have been fitted (continuous lines) with tq = 1.7 x 10"^^ s 
integrating over a gaussian distribution of Seff/^B with mean value 2800 K and width of ~ 300 K 




(4.11) 



inn- 
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Figure 4-17: Peak T after subtraction of the background measured in the anelastic Q and 
^^^La NQR relaxation rate. Also shown is the Arrhenius plot deduced from the three 
maxima. 

In order to gain some insight on the nature of the barrier E^fi, one should develop the 
statistical mechanics of the one-dimensional model of correlated tilts of Markiewicz j^H] , which 
is a rather difficult task. However, one can apply the one-dimensional models of non-linear 
lattice dynamics (92) . which have been used to describe the correlated dynamics of off-centre 
atoms in perovskites |93j . In these models the i-th atom moves in a double-potential of the form 
V = —axf + bxf + cxjXj+i, with minima separated by an energy barrier E = (a'^/Ab) and 
with a bilinear coupling to the neighbors. The coupling constant c takes into account a cluster 
average over configurations. The resulting equation of motion has solitonic solutions similar to 
those found for the octahedra in La2Cu04,|nni and it has been shown j92| l9,Sj that the spectral 
density J (uj) of the displacements x contains a resonant peak at the frequency of vibration 
in each minimum, and a central peak due to the jumps to the other minimum with 
mean hopping frequency t^^ = Tq^ exp (Ees/T) with effective barrier E^s ~ 1.75Ey^2c/a, 
increased with respect to the local barrier E by the interaction between the octahedra, and 
To = ^ \/c/a, where v is the average velocity of propagation of the soliton-like excitation through 
the atoms spaced by d. 

The central peak (see Sec. 12. 2() is therefore responsible for the peak in the NQR relaxation 
rate and in the elastic energy loss; the effective barrier E'eff/^B = 2800 K is compatible with the 
theoretically estimated |42l I48j local barrier of ~ 500 K (see also Sec. I4.1.2[l assuming a mean 
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coupling constant cq ~ lOo- The width of the gaussian distribution of c is cr = O.IScq for the 
NQR and 0.25co for the anelastic data and is justified by the distribution in size and shape of 
the regions where the octahedra clusters build up the cooperative dynamics; such regions are 
delimited by Oi, which blocks the relaxational dynamics. An asymmetry energy AE = 280 K, 
10 times smaller than E'efj , has been assumed to reproduce the increase of the intensity of the 
anelastic peak T, and may also be justified as follows. For nuclear quadrupolar relaxation to 
occur it is sufficient that a solitonic front passes near a ^^^La nucleus producing a fluctuation 
of the atomic positions; for anelastic relaxation to occur, instead, it is necessary that a strain 
fluctuation occurs, and the simple propagation of a tilt wave is not sufficient: referring to Fig. 
I4-I6h . for example, it is necessary that also the size of an LTT region changes with respect to 
the neighboring LTO region, since LTT and LTO regions have different strain (in other words 
it must be AA 7^ in Eq. (|2.15)) ): such a process involves energetically inequivalent states and 
therefore AE 7^ 0. 

The proposed physical picture is therefore that tilt waves can arise, propagate with a velocity 
determined by the coupling c between rows of octahedra, and disappear within regions free of 
defects that can block the tilt degrees of freedom (Oi in the first place, substitutional Sr and 
twin boundaries). The creation and annihilation of such excitations would result in strain 
fluctuations and therefore anelastic relaxation. On the other hand, such tilt waves could extend 
over more than few tens of lattice spacings only in accurately outgassed and defect-free crystals 
and would have little or no correlation along the c axis, making difficult their observation 
by diffraction techniques. The lack of correlation among different planes would be 
the main difference between these tilt solitons and the usual twin boundaries and, 
to my knowledge, La2Cu04 is a unique case exhibiting this type of lattice excitations. The 
conditions necessary for their existence are: i) the existence of unstable lattice modes 
describable by a one-dimensional non-linear equation of motion (therefore having 
solitons as possible solutions); ii) a low density of pinning centers, in order to allow the 
formation and propagation of such excitations to occur. La2Cu04 satisfies both conditions, 
since it has planar arrays of octahedra with little coupling between different planes, and with 
the further possibility of decoupling each plane into rows [201; addition, the level of defects and 
impurities can be lowered enough to let these solitonic tilt waves to develop. I tried to observe 
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similar effects in the isostructural La2Ni04+5, but in that case it is impossible to lower the 
content of Oi below acceptable levels. Another isostructural material with planes of octahedra 
close to an instability is Sr2Ru04+5, but it might present problems with excess oxygen jMj . 
and above all the instability is rotational about the c axis |95j . In that case the rotation of an 
octahedron is strongly coupled to the rotations of all the octahedra in the plane, excluding the 
possibility of a one-dimensional equation of motion; it is also intuitive that when a rotation 
of a single octahedron determines the rotations of the neighboring octahedra in both x and y 
directions, excitations like these cannot arise (see also Fig. l6-4() . Indeed, such rotations would 
be coupled to the cge elastic constant, which however stiffens normally below 220 K [HS]- Also 
in RuSr2GdCu208 the octahedra are tilted about the c axis, and in fact its low temperature 
anelastic spectrum does not show any anomaly of this type Wl\ . 

4.10 Tunneling driven local motion of the oxygen octahedra 

This section is devoted to peak ^4 in Fig. 14-81 also shown in Fig. 14- 15b and indicated as "local 
tunneling", and studied in Refs. [791 198j . For this relaxation process, which I label here as 
LT, it is possible to make considerations similar to those made in the previous section for peak 
T. Its intensity, again better appreciated from a modulus defect up to 10% (Fig. 14- 18b ,). is 
extraordinarily high if one considers that at 10 K only small atomic displacements that driven 
by tunneling may occur. Therefore, a great fraction of the lattice must participate in the 
motion, and, again, the unstable tilt modes of the octahedra are the first candidates, but 
on a local scale (single octahedra or even oxygen atoms). 

A confirmation from this hypothesis comes from the effect of Oil with increasing 5 the 
relaxation amplitude is depressed (Fig. I4-I8b ). since the tilts of the octahedra surrounding Oi 
are fixed to accommodate the interstitial. Still, the effect is much less marked than for the 
collective motion, as can be appreciated from Fig. I4-I5b . where in the curve for 6 = 8 x 10~^ 
peak T is completely absent, but peak LT is only halved in intensity. It also appears that the 
peak shifts to lower temperature with increasing doping, and I will show now that this effect 
has to be attributed to the holes doped by Oj. 

The effect of hole doping is better studied by substitution of La with Sr, which we have 
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Figure 4-18: Effect of excess O on the relaxation processes LT attributed to O tunneling (a) in 
a Sr-free sample and (b) in Sr-doped LSCO, where peak LT is accompanied by the relaxation 
of domain walls in the CSG phase (CSG). 

seen disturbs the lattice much less than Oi, and reveals extremely interesting effects of lattice- 
charge coupling, also providing evidence that such a local motion is dominated by atomic 
tunneling. The introduction of holes above ^ 0.01, however, introduces another relaxation 
process attributed to the motion of the domain walls (hole stripes) in the cluster spin glass 
(CSG) phase, as discussed in Sec. 14.11.11 which has to be distinguished from peak LT. Figure 
14- 18b shows the low temperature anelastic spectra of La2-a;Sri.Cu04 at x = and x = 0.022, 
where the Sr-doped sample presents an additional step around 10 K in correspondence with 
the freezing into the CSG phase. The different nature between the tunneling process, LT, and 
the CSG absorption is put in evidence by the fact that excess O depresses peak LT (in the 
as-prepared state) in both samples, but does not affect the CSG step at all. This fact allows 
the two processes to be distinguished from each other. 

Figure 14-191 shows the (T) curves of outgassed LSCO at four doping levels between 
< X < 0.028. What appears from the first three curves at x = 0, 0.007 and 0.015 is that Sr 
doping has two extremely marked effects on peak LT: i) a shift to lower temperature, ii) a 
narrowing and in) an increase of the intensity. Already at a doping level as low as x = 0.015 
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the maximum of peak LT is below our experimental temperature window, and therefore it is 
impossible to say whether the intensity is still higher than that, already considerable, of the 
curve for x = 0.007. For x = 0.028 even the tail of the relaxation is no more visible. The 
effect of shift to lower temperature is also observable with doping with Oi (Fig. I4-I8b l. but is 
obscured by the concomitant depression due to blocking. 




Figure 4-19: Left: elastic energy loss coefficient of slightly doped LSCO at LHe temperature. 
The arrows indicate peak LT; the steps in correspondence with the vertical bars at Tg [x) are 
due to the onset of the CSG phase. Right: peak LT for x = 0.007 measured at three frequencie, 
after subtraction of a constant backgroung. 

With Sr doping the blocking effect is reduced and it appears that doping has a very 
strong effect on the low temperature relaxation.- it increases its intensity and shifts 
it to lower temperature, which implies an acceleration of the pseudodiffusive dynamics 
of the octahedra. Unfortunately the relaxation could be followed only up to the very small 
doping level oi x = 0.015, and it is impossible to extrapolate to the superconducting region of 
the phase diagram. 

The peak for x = 0.007 is sufficiently well defined to be analyzed, and Fig. I4-I9b shows 
the spectral density J{uj,T) = ■^Q^^ (T) measured on the three vibration modes in double 
logarithmic scale, after subtraction of a common constant background. It appears that all the 
curves collapse together in the high temperature region, as expected from a process that is a 
superposition of Debye relaxations with some distributions g (r) and h (A) of relaxation times 
and strengths. In fact, at high temperature one can discard (wt)^ in the denominator of Eq. 
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(tTTSU) and set J (w, T) = ^Q'^ {T)c^T J drg (r) / dA /i (A) r which is independent of uj. The 
J (uj,T) obtained for x = does not satisfy the same condition, because it is described by a 
broad distribution of r (T) which vary with temperature at a sfower rate than in the case x > 0, 
so that the condition ujt ^ 1 is not satisfied in the high temperature side of the peak for aU 
the elementary relaxations. Therefore, from the temperature position and shape of peak LT 
we deduce that doping changes the mean relaxation rate r^^ from a slowly varying 
function of temperature in the undoped case to a function which increases faster 
with temperature. The high temperature side of J {uj,T) indicates that, for x = 0.007, the 
rate increases faster than above 6 K. 

The above analysis demonstrates that peak LT is due to atomic displacements changing 
at an average rate r~^, and both the low temperature of the maximum and the fact that 
oc T" indicate that the rate is not due to overbarrier hopping with the Arrhenius law, 
but is determined by quantum effects. In the standard tunneling model, the changes of 
average atomic positions with rate are associated with transitions of the tunnel system (TS) 
between its eigenstates |99| . The eigenstates arise from the overlap of the atomic wave functions 
(of a single oxygen or possibly of a whole octahedron in our case) over two or more minima 
of the multiwell potential, like that of Fig. l4-4b : the transitions between the eigenstates are 
promoted by the interaction between the TS and the various excitations of the solid, generally 
consisting of emission and absorption of phonons and scattering of the conduction electrons 
j99j . An important difference between the relaxation process LT and those due to TS in glasses 
is that the latter are characterized by an extremely broad distribution of parameters, mainly 
the tunneling energy t and the asymmetry between the minima of the double-well potential; 
instead, in La2-a;Sr2;Cu04 the geometry of the tunnel systems is much better defined, being some 
particularly unstable configurations of the octahedra. As a consequence, the TLS relaxation 
in glasses produces a plateau in the acoustic absorption and a linear term in the specific heat 
as a function of temperature, whereas in La2-2:Sr^Cu04 we observe a well defined peak in the 
absorption, and no linear contribution to the specific heat has been reported. 



76 



4.10.1 Interaction between the tilts of the octahedra and the hole excitations 

I argue now that the marked acceleration of the local fluctuations with doping (narrowing and 
shift to lower temperature of peak LT) is the manifestation of a direct coupling between 
the tilts of the octahedra and the holes, similarly to the TS in metals, whose dynamics is 
dominated by the interaction with the conduction electrons . The transition rate of a TS is 
generally of the form 

r-i(r)oct2[/ph(r) + /ei(r)] , (4.12) 

where t is the tunneling matrix element and /ph (T) and /d (T) contain the interaction between 
the TS and the phonons and electrons, whose excitation spectra depend on temperature and 
hole density. As shown in Sec. 14.1.11 and Fig. l4-3b . doping reduces the tilt instability and tilt 
angle of the octahedra, and hence makes the minima of the multiwell potential closer to each 
other and reduces the potential barriers between them, resulting in a larger t Such changes 
are however gradual, and a doping as low as x = 0.007 can hardly account for the drastic rise 
of only through an increase of t. The greater contribution has to be attributed to the 
increased interaction with the doped holes, corresponding to the term /d (T) , which is null in 
the perfectly undoped material, but immediately becomes » /ph- Models for the interaction 
between tunneling systems and electrons have been developed and successfully adopted, but the 
electronic excitation spectrum of the cuprate superconductors is certainly different from that 
of metals, and new models for the interaction between octahedra and holes are needed. In fact, 
the asymptotic behavior of J{uj,T) at high and low temperature for x = 0.007 indicate that 
/el (T) ~ T" with n ~ 3 — 5, completely different from TS in standard metals, where /ei (T) is 
less than linear |99j . 

Unfortunately, with the present experiments it has been possible to follow peak LT only 
up to X = 0.015, but it cannot be excluded that the pseudo-diffusive lattice modes giving 
rise to peak LT modify the electron-phonon coupling at higher doping. In fact, the 
characteristic frequencies of these modes are far too low for having any influence on the electron 
dynamics for x < 0.02, but increase dramatically with doping. In terms of the multiwell lattice 
potential, it is possible that at higher doping the barriers between the minima become small 
enough to give rise to an enhancement of the electron-phonon coupling predicted by some 
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models with anharmonic potentials jlUUj . 
4.10.2 Static and dynamic tilt disorder 

The observation of an increased intensity of peak LT at x > with respect to x = (Fig. 
I4-I9t i) provides an explanation for the apparent discrepancy between the tilt disorder observed 
by anelastic experiments on La2Cu04+5 on one side and EXAFS |85j and atomic PDF |38j 
measurements on La2-a;Sr^Cu04 on the other side. The anelastic spectra of Fig. 14- 181 indicates 
a dynamic tilt disorder increasing with reducing 6 (and therefore doping), while the latter 
techniques see the opposite effect of an increasing tilt disorder with increasing doping through 
X. Actually, the EXAFS spectra and atomic PDF are sensitive to both static and dynamic 
disorder; at x = the instantaneous fraction of octahedra swept by the tilt waves is relatively 
small; according to a crude estimate jSjl a few percents of octahedra instantaneously involved 
in the tilt waves account for the intensity of NQR relaxation (peak T) and their effect on the 
PDF or EXAFS spectra would be undetectable. At higher x, the instantaneous tilt disorder will 
increase due to both the disorder in the La/Sr sublattice and the lattice fluctuations. Instead, 
the anelastic spectroscopy is sensitive only to the dynamic disorder with characteristic frequency 
comparable to the measurement frequency; the introduction of interstitial O certainly increases 
the static tilt disorder but also inhibits the dynamic one, resulting in the depression of peaks 
LT and T observed in Figs. I4-I8l and l4-Ill By introducing substitutional Sr, which disturbs the 
lattice much less than interstitial O, it is possible to observe that doping actually increases 
the fraction of fluctuating octahedra, and not only the static disorder. 

4.11 Charge and spin inhomogeneities on nanometer scale 

4.11.1 The cluster spin glass phase: fluctuations of the hole stripes between 
the pinning points 

The low-temperature anelastic spectrum of LSCO presents a step-like increase of the absorption 
below a doping-dependent temperature Tg{x), as already noted in Figs. 14- 18b and 14-191 1. 
Figure I1-2UI clear Iv shows that the temperature at which the increase of Q^^ occurs is in close 
agreement with the temperature Tg (x) for freezing into the cluster-spin glass (CSG) state (Sec. 
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14.31 and Fig. As also discussed in Refs. |1U11 1^ . Tg is not easy to be determined, 

and is often identified with the temperature Tf (uj) at which the rate of the spin fluctuations 
fahs below the frequency co probed by the experiment, > 10^^ s~^ for neutron scattering j58j . 
10^ - 10^ s-^ for NMR/NQR [102], 10^ s^^ for muon spin relaxation ^ 10^ - 10^ s^^ for 
anelastic spectroscopy and < 0.1 s~^ for magnetization measurements jlUlj . A system with 
glassy dynamics is generally describable with a broad distribution of relaxation times r, and 
the maximum relaxation time tm of such a distribution diverges at Tg, often following the Vogel- 
Fulcher law tm = tq exp ^^^^^ ^ , resulting in a maximum of the dynamic susceptibilities at 
ujtm — 1- By adopting the above criterion, one measures the freezing temperatures Tf (uj) such 



that w ^ ~ To exp 



, and the lower is uj the closer Tj approximates Tg. For the case 



_kB{Tf-T,} 

of freezing into the CSG state in LSCO, it seems that tm diverges fast enough that Tf measured 
by magnetization is a good estimate of Tg and can be well approximated by the expression [Jl] 

Tg (x) ~ (0.2 K) /x . (4.13) 

The vertical bars in Fig. 14-201 are obtained from this expression of Tg (x) and have an excellent 
correlation with the rise of Q^^ (T) [10311611 162j . clearly indicating that the acoustic absorp- 
tion is correlated to the CSG state. A strong indication that the absorption step in the 
CSG state is of magnetic origin within the Cu02 planes is its complete insensitiveness 
to interstitial O, as shown in Fig. 14- 18b . whereas the rest of the anelastic spectrum, involving 
tilts of the octahedra, is drastically affected by the presence of even minute amounts of Oi. 

Regarding the mechanism producing dissipation, it might be similar to the well known 
stress-induced movement of domain walls (DW) in ferromagnetic materials [3 13 • In 
this case, the spin order is antiferromagnetic (AF), and the anisotropic strain might be 
correlated with the easy axis m for the staggered magnetization |103j . The elastic 
energy of domains with different orientations of m would be differently affected by a shear stress, 
and the domains with lower energy would grow at the expenses of those with higher energy. 
According to the prevalent interpretation presented in Sec. 14. 4( the walls between the 
AF clusters should be the charge stripes where the holes are segregated, and therefore 
the acoustic absorption below Tg would be due to the fluctuations of the hole 
stripes. According to this view, the rise of (T) would be determined by the freezing of 
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Figure 4-20: Step in the Q below the temperature Tg for freezing into the cluster spin glass 
state at different dopings. 

the spins into AF domains with some disorder in the directions of the staggered magnetization 
m. According to neutron diffraction and magnetic susceptibility, in La2Cu04 or slightly doped 
La2_a:Sra;Cu04 (x < 0.03) the direction of lii is parallel to the h axis but we suppose 
that some distortion of the spin order must be present, like e.g. in Gooding's model 
|53j ■ otherwise the anisotropic strain would be parallel to h in all domains and there would be 
no anelastic relaxation. Note that the fact that the hole stripes must have fluctuations and 
topological defects at all temperatures is predicted by several authors |104U105lllOB| . but what 
is needed here is disorder in the average direction of m from domain to domain. Within this 
picture, the diffraction experiments see the ordered fraction of stripes, with regular spacing and 
separating antiphase AF domains, while anelastic relaxation is due to the disordered 
fraction of the spin structure. The relaxation of DW pinned by impurities, the Sr dopants 
in our case, is generally characterized by a broad distribution of relaxation times, and this 
would explain why a plateau rather than a peak is observed. The (T) step at Tg would 
than be analogous to the Q^^ (T) step occurring at Tt in correspondence with the structural 
HTT/LTO transformation with consequent formation and motion of twin walls (see the step 
#1 followed by the broad peak #2 in Fig. H^ . 

A possible alternative or concomitant mechanism of relaxation within the CSG phase is the 
motion of kinks on the otherwise straight stripes, in analogy with the motion of kinks in 
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dislocations The dynamics of such kinks in the hole stripes has been theoretically analyzed 
j6Uj . also in the presence of pinning, and seems to correspond to another anelastic relaxation 
process attributed to collective stripe depinning, as described in Sec. 14.11.31 

The above results indicate that the anelastic relaxation below Tg is due to the stripe fluc- 
tuations, when they act as walls between frozen AF domains. On the other hand, it has been 
established that the stripes can be locked by the LTT structure, when they are parallel 
to the LTT modulation (for x > 0.055 |59[ I70j) and commensurate with it (x = | |56[ I57j). 
The time scale of the neutron scattering experiments that established this phenomenology is 
extremely short, and what appears static in those experiments may be fast for low frequency 
anelastic spectroscopy; nonetheless, in case of real locking of the stripes by the LTT structure, 
the absorption step at the CSG transition should be depressed of the fraction of stripes that 
become static. For this reason, we extended the anelastic measurements to La2_x'Ba2^^Cu04 
(LBCO) |.61) . that has a greater tendency to form the LTT phase (see Fig. 14- 5 1) . 
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Figure 4-21: (a) Compliance of LBCO and LSCO (rescaled to overlap with that of LBCO) with 
3% Sr. (b) Doping dependence of the intensity of the step at Tg for LSCO and LBCO, 
compared with the relaxation amplitude expected from dislocations (a = 3) or ferromagnetic 
domain walls (a = 2). 

Figure l4-21b shows both the absorption and real parts of the compliance of LSCO and 
LBCO at X = 0.03. The data of LSCO practically overlap with those of LBCO if multiplied 
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by 1.45, and the onset of the decrease of the compUance coincides with Tg, indicating that it 
is determined by the same mechanism producing the absorption, i.e. the freezing of the stripe 
fluctuations. The factor ~ 1.5 can be ascribed to a greater magnetoelastic coupHng in the Ba 
compound, possibly resulting from the shghtly different atomic sizes and distances. 

Figure ll-21b presents the ampHtude of the Q"^ step at Tg for LSCO at x = 0.025, 0.03 
and 0.06, and of LBCO at x = 0.03 and 0.06 normahzed to the value assumed at x = 0.03. 
These data are compared with the relaxation strengths expected from dislocations, Ad(x), 
and from ferromagnetic DW, Adm (x), strongly pinned by a concentration x of defects in 
a bidimensional geometry. In both cases a power law A oc n (x) I" (x) is expected, with I 
the mean length between pinning points, n the line density and with a = 3, in the case of 
dislocations [S] or a = 2 in the case of FM DW jl07j . The A (x) curves have been calculated 
assuming that the stripes are fixed at the pinning points , obtaining A oc 2x I ^ ~ ^) ' 



as explained in Appendix C. In this manner, the experimental data of LSCO are reasonably 
well reproduced with a < 3, intermediate between the dislocation and DW cases. It should 
be remarked that the condition of strong pinning, namely that the walls are really pinned 
at the dopants and can move only between them |107j . is probably not verified. In fact, as 
discussed in Sec. 14.11.31 Morals Smith et al. |60j estimate that the stripes in LSCO are in 
the weak pinning regime, where the stripe may accommodate into the pinning potential over 
a collective length Lc enclosing several pinning centers. Nonetheless, the dependence of 
AQ^^ for LSCO on x is intermediate between that typical of ferromagnetic DW and that of 
dislocations, and this appears as an indication that indeed the relaxation process below Tg 
and therefore the stripe fluctuations in the CSG phase can be assimilated to the 
motion of line defects pinned by the Sr dopants. 

Comparing now the data on LSCO with those on LBCO in Fig. l4-21b . it turns out that 
when La is substituted with 6% Ba instead of 6% Sr, then AQ~^^ is ~ 7 times smaller than 
expected. This depression has to be attributed to the transformation into the LTT structure 
below Td ~ 40 K, as indicated by the upward inflection in the Young's modulus observed only 
for La2.94Sro.o6Cu04 (Fig. I4-10|) . and is a clear indication of pinning of the DW motion within 
the LTT phase, which causes a modulation of the potential felt by the holes, as explained in Sec. 
14.1.21 As explained in Sec. 14.7.21 there are indications of the formation of LTT domains also 
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in LSCO, but only in LBCO and LSCO co-doped with Nd a stable LTT phase with long range 
order sets in, which can provide a pinning potential for the stripes (see Sec. I4.1.2|) . Only the 
stripes parallel or nearly parallel to this modulation should be clamped, and therefore from the 
reduction factor ~ 7 we deduce that about 87% of the DW are parallel to the direction of the Cu- 
O bonds in the LTT phase (or even more if the transition to the LTT structure is incomplete); 
this stripe orientation is in agreement with the direction of the magnetic modulation observed 
by neutron scattering for x > 0.055 |59j . 

The pinning of the DW of the CSG phase within the LTT structure is well documented for 
X ~ |, when the stripe spacing is commensurate with the lattice spacing j5bl I57j . Figure 
demonstrates that almost complete pinning of the stripes can occur also at x = 0.06, 
far from the condition of commensurability with the lattice; in addition, the anelastic 
spectroscopy provides evidence of pinning over a time scale r > uj^^ ~ 10^'^ s, much longer 
than the time scale of neutron diffraction, demonstrating that the pinned stripes are really 
static. 

4.11.2 Electronic phase separation at x < Xc 

The study of the low-temperature anelastic spectrum of LSCO has been concentrated around 
the critical doping Xc — 0.02 above which no more long ranged AF order exists (Fig. l4-3() . on 
samples with nominal doping x = 0.015, 0.016, 0.018, 0.024, 0.03. The prevalent view is that 
for X < Xc a spin glass (SG) state still occurs, but different from the CSG one both in the 
nature (long range AF order with uncorrelated spin distortions near the dopants instead of AF 
clusters) and doping dependence of the freezing temperature Tf {Tf oc x instead of oc 1/x). 
This differences should appear also in the anelastic spectra, since in this picture there are no 
walls in the SG state; therefore, below Xc there should be no step-like rise of the acoustic 
absorption at Tj, or possibly one expects a sharp peak corresponding to those appearing in 
the NQR and /uSR relaxation rates, in correspondence with the crossover for spin fluctuations 
becoming slower than the experimental frequency 17_ (the condition lot ~ 1). On the contrary, 
the anelastic spectra at x < x^, present the same (T) step as for x ~ 0.02, at 
nearly the same temperature but with smaller amplitude |103l I63j . The temperatures 
of the (T) steps found for both x > Xc and x < Xc are plotted as open circles in the region 
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of low-T and low-x of the phase diagram of LSCO in Fig. 14-221 where also the canonical Tf (x) 
and Tg (x) lines are indicated, and there is no hint of a crossover from Tj (x) to Tg (x). 
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Figure 4-22: Low-T and low-x region of the phase diagram of LSCO. The open circles are Tg 
deduced from the step in the (T) curves. 

In view of the difficulties in preparing samples with x determined to better than 0.1% and 
considering the possibility of inhomogeneities in the Sr distribution, we carried out a rather 
accurate characterization of our samples. The Sr content x was checked from the temperature 
Tt (x) of the HTT/LTO transformation and the homogeneity was checked from the width of 
the transformation appearing in the (T) and E (T) curves, as explained in Sec. 14.7.11 The 
deviations of the resulting x from the nominal value were not important, as shown Fig. l4-9b . 
and no inhomogeneity could be detected, to justify regions with x > 0.02 in samples with 
X < 0.018. The magnetic state of the samples was characterized by means of NQR and SQUID 
measurements performed by M. Corti and A. Rigamonti in Pavia: Tn of the samples with x = 
0.015 and 0.016 was ~ 150 K, while for x = 0.018 it was only possible to obtain an upper limit 
Tn < 60 K, due to the overlap with the intense maximum in the ^^^La relaxation rate usually 
attributed to spin freezing. 

It can be concluded that the similarity of the anelastic spectra for x < 0.02 with those 
for X > 0.02 is a real effect and not due to poor samples. The anelastic data are then in 
perfect agreement with the neutron scattering experiment by Matsuda et al. |55j on the magnetic 
correlations in La2_xSra;Cu04 for x < Xc, where it is found that also at x < Xc the 3D AF 
ordered phase coexists below ~ 30 K with domains of "diagonal" stripe phase (with the hole 
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stripes at 45° with respect to the Cu-0 bonds). According to these authors, the hole locahzation 
starting around 150 K involves an electronic phase separation into regions with xi ~ 
and X2 ~ 0.02, and variations in doping affect only the volume fractions of the two phases. The 
Q^^ step in samples with x < 0.02 jl031 IH^ is therefore due to the fraction with X2 ~ 0.02, so 
explaining the fact that the onset of the absorption occurs at the same temperature appropriate 
for the local doping x = X2, while the decrease of the absorption amplitude reflects the decrease 
in the phase with x = X2- 

4.11.3 Thermally activated depinning of the charge stripes 

The elastic energy loss rise in correspondence with the freezing into the CSG phase discussed 
in the previous section, although connected with the stripe fluctuation, exists only below Tg, 
and for this reason is attributed to magnetoelastic coupling with the AF domains rather than 
to the stripes themselves. The charge stripes, however, do exist also at higher temperature, 
e.g. in LSCO substituted with Nd and 12% Sr the superlattice peaks of charge order appear 
below 70 K [HEl, and various types of phase diagrams with charge stripes at temperature well 
above that of the spin glass state are extensively discussed j^l Carrying on the analogy 
between the charge stripes and other linear defects, like dislocations or domain walls, we tried 
to identify a possible anelastic signature of the stripes that overcome the pinning potential by 
thermal activation, and found that a relaxation process with maximum at ~ 80 K for ~ 1 kHz, 
here labeled peak S, has all the expected characteristics _65 ; in Ref. it is also discussed 
why alternative explanations should be very unlikely. The peak is evidenced with a red arrow 
in Fig. 14-231 where the evolution of the anelastic spectrum of La2-aSr^Cu04 from x = 0.018 
to X = 0.20 is shown. The peak well visible up to x = 0.045 is peak T discussed in Sec. 14.91 
and attributed to the solitonic tilt waves of the octahedra; it is progressively suppressed by the 
disorder introduced by Sr. 

For X > 0.125 the temperature Tt of the HTT/LTO transformation becomes lower than 
250 K and is easily identifiable as a sharp rise of Q^^ (T) and softening of the Young's modulus 
(the latter not shown in Fig. l4-23j) . Peak S has an intensity with non-monotonic dependence 
on X and is visible up to x < 0.15; in the overdoped state, x = 0.20, there is no trace of the 
peak, and the only feature of the anelastic spectrum is the HTT/LTO transformation at 70 K. 
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Figure 4-23: Anelastic spectra of LSCO between 50 and 250 K from the underdoped to the 
overdoped state, measured at ~ 1 kHz. The red arrow indicates the peak attributed to the 
stripe depinning. 

The Qmax (^) data for peak S are reported in the lower part Fig. l4-24b . and it is apparent 
that the intensity of peak S has exactly the doping dependence expected from a 
polaronic relaxation due to the hole stripes: i) it vanishes at x — > 0, where there are no 
holes and therefore no stripes; ii) it vanishes in the overdoped region, x > 0.15, where the 
holes are uniformly distributed as in a normal metallic state; in) it presents an abrupt drop 
(note the logarithmic scale of Fig. I4-24|1 between at x = 0.045 and 0.08, in correspondence 
with the change of the stripe order from parallel to diagonal with respect to the 
LTO lattice modulation j59| I7()j: iv) it has a local maximum around x ~ |, where the 
commensurability between stripe spacing and lattice enhances the stripe-lattice 
coupling, and the sharp minimum exactly at x = | is most likely due to the fact that 
the locking between LTT lattice fluctuations and stripes is so strong that the fluctuations are 
depressed. The agreement between the Q^ax (^) data and what is expected from a relaxation 
process due to fluctuations of pinned hole stripes is certainly remarkable, but there is more. 

Figure 14- 24b shows the anelastic spectrum for x = 0.015 from 40 to 300 K measured at 
three vibration frequencies, which allows peak S together with the neighboring peaks to be 
reliably fitted, in order to extract the various parameters of the relaxation S and particularly 
its activation energy, that can be identified with the effective pinning barrier Ep (the origin 
of the peak around 60 K is not yet identified). The continuous curves are a simultaneous fit at 
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Figure 4-24: (a) Example of fit of the anelastic spectrum of LSCO with x = 0.015, in order 
to deduce the activation energy and intensity of the peak attributed to stripe depinning. (b) 
Doping dependence of the intensity and activation energy of the peak at 80 K (for ~ 1 kHz) 
attributed to the stripe depinning. The blue line is the predicted collective pinning barrier, 
according to jHO] • 

all frequencies with Eq. (|2.48|) for all three peaks. Peak S is rather broad and symmetric, with 
a = 13 = 0.7, To = 7 X 10"^'' s and Ep/k^ = 1690 K. Fits like this can be done up to x = 0.045, 
while for higher doping there is a larger error in extracting the parameters of peak S, since it 
is smaller and masked by the other processes. Nevertheless, the doping dependence of tq and 
Ep/k-Q parameters could be followed over the complete doping range. The dependence of Ep 
on X is particularly interesting, since it closely reflects that of \ogQ^^^, as demonstrated 
in Fig. 14-241 Also the pre-exponential factor tq has a general decreasing trend with increasing 
doping, passing from 1.4 x 10~^^ s to lO"^'' s between x = 0.015 and x = 0.08; for this reason 
the temperature of the peak measured around 1 kHz remains close to 80 K, but it is understood 
that the peak temperature depends on frequency. Values of tq as low as lO^^'' s are certainly 
too long for a polaronic relaxation, but may be understood in terms of an extended defect like 
the stripe, with a strong lattice contribution, as argued in what follows. 

Before discussing the analogies between the doping dependences of relaxation strength and 
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activation energy, the present results are compared with the predictions of the model of Morals 
Smith et al. for the dynamics of the pinned stripes jHO]. From the balance between the 
electrostatic interaction Ec of the hole stripe with the Sr^"^ dopants, which tends to pin the 
stripe into a curved configuration, and the elastic energy ei for bending the stripe, which tends 
to keep it straight, both the collective length Lc for the stripe reconformations and the 
pinning barrier Ep for such reconformations to occur are estimated. It turns out that 

L, = ax-'>/'^{ei/e,f'^ (4.14) 

where £c = e^/ (eoo) ~ 0.1 eV is the Coulomb energy scale with an isotropic dielectric constant 
eo — 30 and a the lattice parameter, while = J ~ 0.1 eV is the elastic energy of the 
stripe determined by the AF exchange constant J. Setting ei = J tacitly assumes that the 
hole stripe fluctuates in a static AF background of spins; this is appropriate when the spin 
fluctuations of the Cu^^ atoms are frozen below the characteristic stripe fluctuation frequency 
under study, as in the CSG state, but certainly not at 80 K for x > Xc- I suppose that the 
spin fluctuations would lower the elastic energy ei of the stripe, and therefore also Lc, possibly 
making a crossover to the strong pinning regime with Lc < I — aj \fx , where I is the mean 
distance between pinning centers. Otherwise, with the above estimate of ei ~ Sc the collective 
pinning length — ax~^l'° is always larger than /, resulting in the weak pinning regime. 
The collective pinning energy is estimated as |60j 

Ep = {eiel) ~ (1300 K) x'^^^ . (4.15) 

In view of the rough nature of the estimate, especially of £c, it is reasonable to expect that the 
dependence of Eq. 1)4. 15() can be adopted, but with the constant factor as an adjusting 
parameter. The blue line in Fig. 14-241 is a fit of Ep (x) of peak S, excluding the 2; ~ | region, 
with Eq. ()4.15|) . and the constant is found to be (e/^c) /k^ = 840 K. It can be concluded 
that the agreement between the general trend of Ep (x) and the predictions for the 
collective pinning energy barrier of stripes pinned by the Sr dopants of Ref. |60| is 
surprisingly good, and provides further support to the interpretation of peak S in terms of 
stripe depinning. It should be noted, however, that the treatment of Ref. ^60j does not include 



88 



the interaction with the lattice, which is at least as important as the electrostatic in interaction 
with the dopants, as argued in the following Section. 

4.11.4 Hole-lattice coupling and octahedral tilts 

Regarding the mechanism producing peak S, it would be tempting to simply identify the hole 
stripe with a dislocation or domain wall, and apply the models of anelastic relaxation from these 
linear defects, as proposed for the relaxation in the CSG state in Sec. 14.11.11 The problem 
with peak S is that at 80 K the hole stripes do not separate domains with different orientations 
of the anisotropic strain, as is the CSG state, since at this temperature the spin fluctuate very 
fast U7j- Therefore, it is not possible to associate an anelastic strain e^^ proportional to the 
area swept by the stripe (the continuum equivalent of the elastic dipole of a point defect), as 
in the above models. Instead, different strain contributions e^"^ should be associated 
with the different configurations that the stripe may assume by thermal activation 
over the pinning barrier, e.g. straight or curved or kinked stripe segment. It is likely that 
e^"^, like the distortion coupled to the presence of holes [HZj, is mainly connected with the 
in-plane shear also associated with the tilts of the octahedra. Then, it is reasonable 
to assume that e^'^ is an increasing function of the local degree of tilting of the octahedra. 
This would partly contribute to the general reduction of the relaxation strength with increasing 
doping, because both the in-plane shear strain and tilt angles 9 of the octahedra are decreasing 
functions of doping (see Fig. l4-3() . 

It has already pointed out how the sharp variations of Qmix with doping correspond to a 
variation in the degree of coupling between lattice and stripe; it has to be explained why also the 
activation energy Ep presents similar features. An explanation can be found by observing that 
the octahedral tilts (the most unstable lattice modes) are responsible for the stripe-strain 
interaction, and therefore determine the relaxation strength, but also cause stripe pinning 
from the lattice, as demonstrated by the well known locking between LTT tilt pattern and 
stripes. Therefore, large tilt amplitudes are associated with both large relaxation strength and 
strong pinning from the lattice; in other words, anomalies in the tilt-stripe coupling result in 
anomalies in both Q^ax ^"^^ ^p- 



89 



4.11.5 The picture of the slow stripe fluctuations after anelastic spectroscopy 

The two anelastic relaxation processes we attribute to the stripe motion are again shown in Fig. 
l4-25b , in a sample with x = 0.03. In a first instance jl03l 1611 165j we assumed a strong pinning 
regime, where the CSG absorption below Tg is due to the motion between pinning points and 
the 80 K relaxation as the thermally activated depinning process. According to the analysis of 
Morais Smith et al. the weak or collective pinning regime would be verified also at very 
low doping, meaning that there is no clear distinction between motion of stripe segments that 
are between or at the Sr dopants; rather, there is a collective length (x) characteristic for 
the stripe displacements that minimize the free energy, as mentioned in Sec. 14. 11. H and 14.11.51 
If this type of motion produces peak S around 80 K, its contribution becomes totally negligible 
at Tg < 15 K. 




Figure 4-25: (a) Features of the anelastic spectrum of Lai.97Sro.o3Cu04 associated with the 
dynamics of the charge stripes: freezing below Tg and thermally activated depinning above 
80 K. (b) Sketch of how the absorption below Tg might be decomposed into the contribution 
from DW fast local motion multiplied by the fraction of domains that are frozen. 

This is schematically illustrated in Fig. l4-25b with double logarithmic scale, where the 
extrapolation of peak S is clearly negligible at Tg. The relaxation appearing below Tg must 
therefore originate from a much faster and probably more local type of stripe motion, indicated 
as a (totally hypothetical) broad maximum. The fact that the elastic energy loss is seen to grow 
only below Tg indicates that, while the dynamics may be associated to the stripes, the mecha- 
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nism producing anelasticity involves the frozen AF domains, as discussed in Sec. 14.11.11 This is 
represented in Fig. l4-25b with the fraction of domains between stripes having a characteristic 
fluctuation frequency r^^ frozen below the measurement frequency m; the resulting Q^^ (T) 
curve is the product of this fraction with the broad absorption that would be hypothetically 
measured if the AF spins where already frozen at higher temperature, allowing the mechanism 
of the stress induced domain wall motion to be operative. The fact that this absorption is broad 
is indicated by the shape of Q^^ (T) below Tg. 

The mechanism governing the dynamics of the domain wall or stripe motion at these tem- 
peratures is necessarily different from the thermally activated overbarrier motion of peak S, and 
must be totally quantum mechanical. To my knowledge there are no theoretical treatments of 
this problem; the stripe relaxation time in the quantum limit r ~ 10~^^ s estimated in Ref. 

based on ohmic dissipation [/ei oc T in Eq. (|4.12|) ]. is far too short for explaining the 
anelastic relaxation below Tg at our frequencies; in fact, the Debye factor <^t/[1 + (ujt)'^] would 
cut the absorption amplitude by a factor of 10"^'', making it unobservable. On the contrary, the 
broad shape of the (T) curve below Tg indicates that a consistent weight of the distribution 
function of the relaxation frequencies is in the kHz range. 
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Figure 4-26: Phase diagram of LSCO with the dynamic ranges of the stripe motions deduced 
from the anelastic experiments. 

Figure 14-261 is the phase diagram of LSCO (see Fig. 14- .Sj) with the dynamic ranges of 
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the stripe motions deduced from the anelastic experiments. The depinning rate curves have 
been calculated from the relaxation time r (T) deduced from peak S; they correspond to the 
mean relaxation rate for the stripes to overcome the pinning potential from the Sr^"*" dopants, 
presumably in the weak pinning regime. As explained in the previous Sections, the relaxation 
rate reflects the abrupt changes of the effective pinning barrier in correspondence with changes 
in the stripe-lattice coupling, like the change of the stripe orientation at x = 0.055 or the 
possible insurgence of LTT fluctuations near a; = |. The black symbols are the temperatures 
for the onset of the incommensurate spin order, deduced from the superlattice peaks in magnetic 
neutron scattering; for x < 0.1 these temperatures correspond to Tg (^x,ll! ~ 10^^ ^~^)' namely 
the onset of the CSG phase on the neutron scattering time scale pHl. The CSG phase on the 
(quasi)static time scale from magnetic susceptibility measurements is the gray region. The 
points labeled "fluctuations of pinned DW" and the two accompanying upper and lower lines 
indicate the rise of the acoustic absorption (10%, inflection, 90% of the step) in the CSG 
phase. The absorption is due to the stress-induced motion of the AF clusters separated by the 
hole stripes, and in these temperature region there is a substantial spectral weight of stripe 
fluctuations with frequencies in the kHz region. 

Before concluding I mention a recent ultrasonic experiment jl08j . where the influence of Nd 
codoping on LSCO has been studied. It is found that an increase of the Nd content causes 
an increase of the stripe pinning barrier deduced from peak S, and eventually the formation of 
LTT phase clamps a fraction of stripes, as expected. 

4.12 Summary of the main results obtained in LSCO 
4.12.1 Structural transformations 

The main feature of the anelastic spectrum of La2_xSra;Cu04 is the HTT/LTO structural 
transformation at Tf (x) (Sec. I4.7.1|) . which produces steps in both the real and imaginary 
modulus, accompanied by some broad feature a little below Tt, due to the motion of the 
twin walls. The analysis of the position and width of this transformation has been useful for 
checking the actual doping level x and homogeneity, especially for x < 0.03. The LTO/LTT 
transformation is much more subtle, but can be identified as a rise of the modulus (see e.g. 



92 



the upper panel of Fig. I4-1U() occurring at the expected Td, as already observed in ultrasound 
experiments j47j . 

4.12.2 Absorption in the CSG state 

The other feature that has an evident relationship with other types of experiments is the step- 
like rise of (T) below the temperature Tg at which the Cu spins are frozen into the Cluster 
Spin Glass state (Sec. I4.11.1() . as shown e.g. in Fig. I4-2UI The occurrence of a clear (T) 
step below Tg convinced us that there must be sufficient coupling between spin degrees of 
freedom and strain to produce visible effects in the anelastic spectra. Considering the prevalent 
picture of the CSG state as due to antiferromagnetically correlated spin clusters separated by 
hole-rich walls (or stripes), it is obvious to assign the mechanism of elastic energy dissipation 
to the movement of the walls between spin clusters. The effect is analogous to the 
Q-i ^ji-j g^gp observed below Tj at the HTT/LTO transformation: in the case of the structural 
transformation the fluctuating anelastic strain is provided by the different orientations of the 
orthorhombic axes of different domains, while below Tg it is due to anisotropic magnetoelastic 
strain coupled with the direction of the staggered magnetization in each spin cluster. Having 
confidence in the correlation between Tg and Q^^ (T) step, and in the homogeneity of the 
samples evaluated from the narrow width of their HTT/LTO transition, it has been possible to 
show that (Sec. I4.11.2|) . contrary to the generally accepted opinion, below the critical doping 
Xc — 0.02 the CSG phase does not disappear in favor of a spin-glass state, but there is a 
phase separation between CSG and AF hole-poor states, as also indicated by detailed neutron 
spectroscopy studies |55j . The anelastic experiments are also in agreement with the observation, 
mainly after neutron spectroscopy, that the hole stripes, or equivalently the spin walls, are 
mostly parallel to each other and parallel to the direction of the Cu-0 bonds for x > 0.55, so 
that they are clamped by the lattice modulation of the LTT phase. In fact, the amplitude of 
the (T) step in a sample transformed into the LTT phase was reduced by almost 90% with 
respect to those in LTO phase, consistent with the hypothesis that the majority of walls is 
clamped by the static LTT lattice modulation and therefore does not contribute to the 
dynamic compliance. 
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4.12.3 Peak S 

Having found such close correlations between the amplitude and temperature of the Q^^ (T) 
step below Tg as a function of doping and the phenomenology of the hole stripes deduced from 
several types of experiments, the next step was to proceed with the analogy between these 
stripes and other linear defects and look for possible thermally activated depinning of the 
charge stripes from the Sr^+ dopants. This kind of process has been theoretically studied 
but never observed experimentally; this is understandable, since the dynamic response of the 
single Cu spins and holes dominates most spectroscopies, like magnetic, NMR, dielectric and 
optical. The anelastic response, instead, is insensitive to the charge and spin dynamics, except 
for those excitations or structures that are coupled to strain, as linear charge structures should 
be. The theoretical estimate t60} for the (collective) pinning energy of the stripes to the Sr 
dopants is Ep/k-Q > 1300 K; therefore, the possible signature of stripe depinning is a thermally 
activated peak with effective barrier ~ Ep. In addition, the intensity of this peak should 
vanish at zero doping, where there are no holes, and for x > 0.2, where the holes are uniformly 
distributed in a metallic state. The analysis of tens of anelastic spectra from samples with 13 
different Sr doping levels, allowed this process to be identified with peak S in Fig. 14-231 besides 
the above requirements, the doping dependences of its intensity and activation energy exhibit 
other outstanding features shown in 14-241 like a drop at x > 0.55, where the stripes have been 
observed to change orientation from parallel to diagonal with respect to the lattice modulation, 
and a local increase at x ~ |, where an instability toward the LTT structure exists, whose 
modulation would be again parallel to the prevalent stripe direction and commensurate with 
their mean spacing. All these features have been explained in Sects. l4.11.3l and l4.11.H in terms 
of coupling of the hole stripes to the lattice through octahedral tilting. 

4.12.4 Peaks Ol and 02 

The assignment of the other relaxation peaks found in LSCO is easily accomplished after the 
identification of the peak due to the presence of interstitial oxygen (Oi) in La2Cu04+5; this 
follows straightaway from the dependence of the intensity of the pair of peaks 01 and 02 in Fig. 
14-1 II on the content 6 of excess oxygen. Peak 01 is therefore assigned to hopping of isolated Oi 
atoms, while peak 02, with slightly higher activation energy and prevailing at higher 6, to Oi 
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pairs. These peaks allow the presence of Oi well below 5 ~ 10~^ to be monitored, as estimated 
from oxygen absorption and evolution experiments (Sec. I4.8.2() : the fact that peak 02 prevails 
over 01 even for 6 <ti 0.01 is explained by observing that the anisotropy of the elastic distortion 
due to Oi (and therefore the intensity of peak 01) is proportional to the small deviation of 
the orthorhombic structure of La2Cu04_|-5 from tetragonal (Fig. I4-13() . while stable Oi pairs 
are highly anisotropic defects (Sec. 14.8.1)1 . Incidentally, peaks 01 and 02 allow a very precise 
determination of the hopping rate of Oi to be made and the existence of stable Oi pairs to be 
assessed. 

4.12.5 Peak T 

Once established how to obtain La2Cu04+5 virtually free from excess oxygen, we found that the 
anelastic spectrum of stoichiometric defect-free La2Cu04 contains two extremely intense ther- 
mally activated relaxation processes, causing strong absorption peaks and modulus softenings 
up to 25% above 150 K (at 1 kHz, #4 in Fig. ]^oi peak T) and 10% above 10 K (#4 in Fig. 
14-81 or peak LT). These relaxation processes are hindered by doping and therefore must be due 
to intrinsic lattice mechanisms, which must have to do with the most unstable lattice modes 
of LSCO: the octahedral tilting giving rise to the LTO and LTT structures (considering the 
low relaxation frequency involved, r^^ i_o 

10^ s"\ by tilting I mean switching among the 
several potential minima of the octahedra, as in Fig. l4-4b . and not continuous tilting). In this 
respect, the comparison with Nd2Cu04 is particularly significant; in Nd2Cu04 the O atoms 
that in La2Cu04 form the apices of the CuOg octahedra are shifted to the interstitial positions 
(Fig. 14-7)1 ■ so that there are no octahedra, and in fact the anelastic spectrum of Nd2Cu04 
is completely flat. The evidence that octahedral tilts are involved comes from the fact that 
peak T is observed also in the ^^^La NQR relaxation rate (Fig. I4-17|1 . which is sensitive to 
fluctuations of the La-0 distances, and therefore to fluctuations of the octahedra. Additional 
clues to the nature of these relaxation processes come from their dependence on doping through 
substitution of La with Sr or through addition of interstitial O; the latter type of doping is more 
effective in blocking the tilts of the octahedra, and in fact 5 < 0.001 is sufficient to completely 
suppress peak T (Fig. I4-I5b ). while the effect of Sr is much weaker (Figs. 14- 15b and 14-2^ ). 
The fact that as little as 6 < 0.001 of Oi completely blocks the octahedral motion giving rise 
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to peak T demonstrates that this a collective type of motion, requiring the coordinated motion 
of hundreds of octahedra, and that can be identified with the solitonic tilt waves predicted 
by Markiewicz (HOI- It should be noted that this type of motion is possible only for planar 
coordination of the octahedra with little correlation among different planes, with the additional 
hypothesis that it is possible to start from a ground state like the LTT tilt pattern that can 
be decomposed in weakly correlated rows of octahedra (Fig. I4-16|) : in this manner Markiewicz 
obtained a one-dimensional non linear equation of motion for the tilts of the octahedra, which 
has solitonic solutions. On the contrary, with rotation patterns like that in RuSr2GdCu208 
(Fig. l6-4j) the rotation of each octahedron is strongly coupled to those of all the neighboring 
octahedra, so that no low energy excitations of this type are possible. 

4.12.6 Peak LT 

The other intrinsic lattice relaxation mechanism in ha2-xSTxCu04^ is peak LT below 10 K 
(Fig. I4-I5b and l4-18() . Contrary to peak T, it is only gradually suppressed by the presence of Oi, 
indicating that the motions involved are of local type, so that octahedra far from the Oj atoms 
are not affected; I imagine that even single O atoms in the Cu02 planes may be involved in such 
a fast relaxation. The most striking feature of peal LT is its evolution with hole doping; in fact, 
the suppression of the intensities of peaks T and LT is mainly due to steric effects, especially 
for Oi {6 ~ 0.001 of Oi suppresses peak T but contributes very little to hole doping). Figure 
I4-I9t i shows that hole doping even enhances the intensity of peak LT, and above all greatly 
shifts to lower temperature and narrows it, indicating an enhancement of the magnitude of the 
fluctuation rate and of its temperature derivative. This phenomenology can be compared 
with the anelastic relaxation due to tunneling of interstitial H in bcc superconducting metals 
or due to two-level systems in insulating and metallic glasses ^ . In such cases the anelastic 
relaxation is due to transitions of the tunnel systems promoted by the interaction with phonons 
and conduction electrons, whose contributions to r^^ (T) are well known; it is also established 
that the direct interaction of the tunneling atoms with the conduction electrons is the main 
responsible for such transitions. It is therefore clear that peak LT is due to single octahedra or 
O atoms tunneling between nearly equivalent potential minima, and that the transition 
rate within such tunnel systems (the relaxation rate t~^) is enhanced by the interaction with 
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the holes. Simple plotting of peak LT in double logarithmic scale shows that r~ oc with 
n ~ 5.4 (Fig. I4-I9b ). which is completely different from the (T) laws known for metallic 
systems; this is reasonable, since the hole excitations in these cuprates are certainly different 
from those in metals, and the analysis of peak LT would provide information on their excitation 
spectrum. 
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Chapter 5 



YBCO 

The family of the YBa2Cu306+x superconductors (often abbreviated in YBCO or Y-123 from 
the cation stoichiometry) is extensively studied because of the particularly high superconducting 
temperature Tc ~ 92 K which make various applications possible, including large scale electric 
devices and wires 1_. On the other hand, YBCO is, among the cuprate superconductors, 
the one where nonstoichiometric oxygen plays the most fundamental role, since it is the only 
responsible for hole doping, but also has a wide stoichiometric range, < x < 1, and high 
mobility. The diffusive oxygen jumps give rise to reorientation of the associated elastic dipole 
and therefore the anelastic spectroscopy is particularly useful in studying these cuprates. 

5.1 Structure and phase diagram 

The structure of YBa2Cu306+x is shown in Fig. 15-11 The ideal YBa2Cu306 has planes of Cu^"*" 
and the Cu02 planes in the Cu^"*" oxidation state; there are no mobile charges, the material 
is semiconducting and the Cu^"*" spins order antiferromagnetically below Tn > 320 K. The 
CuOj: planes may accommodate O atoms up to x < 1 and the stoichiometry can be controlled 
through temperature and O2 partial pressure. The hole doping in the Cu02 planes is due to the 
oxidation from these nonstoichiometric O atoms in the CuO^; planes and depends also on their 
ordering. In fact, the O atoms tend to form parallel Cu-0 chains. The right hand side of Fig. 
15- II shows two cells of the ideal orthorhombic structure of YBa2Cu307 with the excess O atoms 
perfectly ordered in Cu-0 chains along the b axis. In practice only intermediate stoichiometrics 
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0.05 < X < 1 are achievable, and the x — T phase diagram (Fig. I5-2|) contains three main phases: 

i) the orthorhombic-I (O-I) phase with parahel Cu-0 chains along the b axis when x ~ 1; 

ii) the orthorhombic-II (O-II) phase with Cu-0 chains alternately filled and empty when 
X ~ 0.5; Hi) the tetragonal phase at x < 0.3 or high temperatm'e, where the oxygen chain 
fragments are very short and oriented along both the a and b directions. Additional phases, e.g. 
O-III with a pattern of filled-filled-empty chains, may be obtained by carefully equilibrating 
the sample in controlled O2 atmosphere. 



Figure 5-1: Structure of YBa2Cu306+a;- Left: one cell of tetragonal YBa2Cu306; right: two 
cell of orthorhombic YBa2Cu307. The vacant O sites in the CuOx planes are indicated as white 
atoms. 

The structural phase diagram has been mainly determined by neutron diffraction investi- 
gations jllOl I23j . and is shown in Fig. 15-21 in a schematic version. Numerous studies have 
also appeared where the phase diagram is reproduced by Monte Carlo simulations of models 
with different interaction energies of the O atoms up to at least the next nearest neighbors. 
The most successful model is the ASYNNNI (asymmetric next nearest neighbors interactions) 
|23| mil I112j . Such models have also been extended to explain the hole doping in the Cu02 
planes, which depends on oxygen ordering |109j as well as content in the CuOa; planes, and are 
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Figure 5-2: Simplified phase diagram of YBCO (from Ref. |in9j ). 

too numerous to be reviewed here. It will be sufficient to mention that they assume at least 
three interaction energies V2 and V3 for nearest, 2nd and 3rd neighboring O atoms, with 
Vi > and V2 < in order to reproduce the tendency to form parallel Cu-0 chains, and V3 > 
in order to reproduce the O-II phase with alternately filled and empty chains. 

5.2 Oxygen ordering and charge transfer between chains and 
planes 

The mean length of the Cu-0 chains may be estimated from the NQR ^^Cu spectra, where 
different peaks are found for Cu atoms which have zero, one or two neighboring O atom in the 
CuOx plane, and therefore are 2, 3 or 4- fold coordinated, including the apical O atoms jll3j : 
a quantitative analysis of the spectra therefore provides the concentrations uq, rii and 712 of 
such atoms. These are related with the mean chain length and degree of oxygen ordering; in 
fact, for infinitely long chains one has n2 = x, ni = 0, uq = 1 — x, while in the limit of a small 
concentration of isolated O atoms one has 77-2 = 0, ni = 2x, uq = 1 — 2x. The concentration 
and nature of the holes doped by the nonstoichiometric O^^ ions may be studied by X-ray 
absorption spectroscopies jll41 1115] , and it turns out that for x ^ Cu in the empty chains is 
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in the Cu"*" state, while in the Cu02 planes it is in the Cu state, so that the ideal undoped 
condition is 

YBasCusOe = Y^+Ba^+O^" (Cu^+O^")^ Cu^+ ; (5.1) 

the addition of isolated oxygen in the CuO^ planes oxidizes all the pairs of neighboring Cu 
atoms into the Cu^"*' state, but chain fragments with n > 1 consecutive O atoms would require 
the oxidation of n — 1 Cu atoms to the Cu^"*" state, which is not observed jll61 1113j (see also 
Fig. I5-9|) . Therefore it may be supposed that in a first step there are n — 1 O^^ atoms, or n — 1 
holes on the atoms in the chain; such holes may hop among the O atoms of the chain, 
providing a minor contribution to the electrical conductivity. This behavior may be rationalized 
in terms of electrostatic repulsion between the holes jll7j : the repulsion between two holes on 
the same atom (Cu^^ = Cu^"^ + 2/i*) is too strong and they distribute over one and one 
Cu^+. With increasing the chain length n, the electrostatic repulsion between the holes in the 
oxygen chain atoms drives part of them, m, into the Cu02 planes where they delocalize giving 
rise to most of the electric conduction and superconduction. A model for the charge transfer 
from the chains to the superconducting planes requires the knowledge of how many m holes 
are transferred from a chain fragment of length n and has been developed by Uimin jll81lllJ] . 
based also on the experimental observation that for long chains it \s m/n ~ 0.7 jll9j . For short 
chains of length n = 2, 3 and 4, it has been estimated that m = 0, 1, 2 holes are transferred to 
the planes, and the transfer proceeds with increasing n up to the optimal value m = 0.7n. For 
the interpretation of the anelastic spectra, it should be noted that the migration of an isolated 
O atom in the CuO^ plane does not change the oxidation state of the crystal, while already 
the formation of an oxygen pair causes the oxidation of the two neighboring Cu atoms, and 
the joining to longer chains involves charge transfer with the Cu02. Therefore, the jumps of 
isolated O atoms are expected to be much easier than those involving the joining to or coming 
out of a chain fragment, since the latter involve substantial energy changes of the electronic 
system. 



101 



5.3 Diffusive dynamics of oxygen in the CuO^; planes 



While the features of the phase diagram have been extensively studied both experimentally and 
theoretically, the knowledge of the diffusive and ordering dynamics of oxygen is less detailed. 
The macroscopic diffusion has been studied by tracer methods as recently summarized in jl2U| . 
and by in- or out-diffusion in a gas atmosphere in connection with thermogravimetry or electrical 
resistance measurements to monitor the time dependent oxygen content jl21j . The first type 
of experiments provides an average macroscopic diffusion coefficient, which may depend on 
the microstructure and possibly on inhomogeneous oxygen ordering, while the latter type of 
experiments is also heavily affected by the kinetics of the oxygen exchange with the gas phase 
and requires a model for the dependence of the electrical resistivity on the filling of the CuOx 
planes. All these factors heavily influence the apparent diffusion coefficient: for example, the 
apparent activation energy for chemical diffusion may vary from 0.4 eV to 1 eV within the 
same study, depending on microstructure |122j . and even more from experiment to experiment; 
an important role is attributed to the formation of an oxygen rich shell in each grain, that 
obstructs further oxygen in-diffusion during absorption experiments |122j : a surface barrier of 
1.7 eV has been estimated for out-diffusion |123j . The spread of the results on oxygen diffusion 
in YBCO in the literature is even more impressive (several orders of magnitude) in terms of 
the diffusion coefficient or hopping rate at a fixed temperature. The anelastic results have the 
great advantage that, once identified the elastic energy loss peak due to oxygen hopping, from 
the condition lot ~ 1 at the peak temperature, Eq. 1)2. 30() . a relaxation time r very close 
to the oxygen hopping time in the bulk can be measured, even when interpreted by different 
microscopic mechanisms such as vacancy-assisted jl24j . Zener pair relaxation ^125 , interacting 
elastic dipoles in fact, different mechanisms may change the ratio between relaxation and 
hopping time by less than one order of magnitude. Instead, the evaluation of the local barrier 
for hopping depends on the type of model that is assumed, ranging from 0.68 eV for a KWW 
relaxation |126j to 1 — 1.4 eV for Curie- Weiss interactions jl27j . 
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Figure 5-3: Comparison between the chemical diffusion coefficient of O in YBCO from var- 
ious permeation experiments (a) and the hopping time of O deduced from several anelastic 
experiments (b). The figures are from Ref. ( |121j ). 

5.4 Anelastic measurements of the oxygen diffusive jumps 

The anelastic spectroscopy is certainly the best method to study in detail the dynamics of the 
mobile O atoms in the CuOx planes. In fact, as explained in detail in Sec. 12.4.21 an anisotropic 
elastic dipole is associated to each O atom, which reorients by 90° after a jump; it is then 
possible to selectively probe different hopping processes in different environments, through the 
analysis of the distinct maxima in the (T) curve, occurring when the condition lot (T) = 1 
is met for each process. Unfortunately, the present Thesis reports only qualitative results of 
the high temperature peaks, because for technical reasons the measurements were made in high 
vacuum, and oxygen loss occurred during the measurement of the anelastic spectra above 500 K. 
Therefore, the peak shapes were affected by the oxygen loss and could not be reliably analyzed 
in order to extract precise information on the oxygen dynamics. Still, the results clearly show 
the presence of two completely distinct hopping regimes: i) one that involves barriers of ~ 1 eV 
and has been studied by several authors with various techniques, ii) a much faster hopping over 
a barrier about 10 times smaller, that is observed only at the lowest oxygen contents. 
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5.4.1 Anelastic spectra at different oxygen concentrations 

The anelastic spectra of YBCO change completely with varying x between 1 and ~ 0. Figure 
l5-4b presents the anelastic spectra of YBCO measured at ~ 1 kHz at three representative 
values of x: 0.1, 0.5 and 0.9. Each of them contains a different peak; for compatibility with the 
labeling originally given [12111121113011121], they are called respectively P2 (at ~ 80 K), PHI 
(at ~ 550 K) and PH2 (at ~ 750 K). All of them are thermally activated, P2 with an activation 
energy of 0.11 eV, and PHI and PH2 with ~ 1 eV. The upper panel represents the CuOa; 
plane with x increasing from left to right; the oxygen jumps causing the three main anelastic 
relaxation processes are put in evidence. 



P2 PH2 




Figure 5-4: (a) CuO^ plane with x increasing from left to right; the O jumps causing the three 
main anelastic relaxation processes P2, PHI and PH2 plotted in (b) are put in evidence. The 
intensities of these peaks are plotted in (d) as a function of x, and are strongly correlated with 
the probabilities of the corresponding clusters plotted in (c). 

Starting from the superconducting phase with x ~ 0.9, one finds peak PH2, which disappears 
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on lowering x and, according to the pendulum experiment of Xie et al. jl24j . also when re — > 1. 
For this reason those authors suggest a vacancy mechanism for the oxygen diffusion in the O-I 
phase, as represented in Fig. l5-4b . and associate peak PH2 with such a hopping mechanism. 
For 0.2 < X < 0.8 peak PHI is observed, with intensity strongly dependent on thermal history; 
as shown later in Sec. 15.4.51 peak PHI cannot be attributed to hopping of isolated O atoms 
and it must therefore be associated with jumps of oxygen aggregated into chain fragments with 
the exclusion of those in the ordered O-I phase, which produce peak PH2. Finally, at the 
lowest values of x attainable, only peak P2 is observed, which is therefore associated with the 
jumps of isolated O atoms. Additional support to these assignments comes from the comparison 
between the intensities Qmax(^)i) of three elastic energy loss peaks (Fig. I5-4H ) and the 
probabilities p{x,j) of the various oxygen clusters of type j = P2, PHI, PH2 or nothing 
(Fig. l5-4h ) . These cluster probabilities have been reproduced from Ref. [lUyj, where they are 
calculated with Monte Carlo simulations on the ASYNNNI model assuming T = 300 K. The 
clusters are represented as insets in the various plots oi p{x,j) and labeled according to the 
relaxation process j that they cause according to our interpretation; for PHI the 5 clusters in 
the three plots have to be summed together. The Q^ax {^^j) curves have been plotted from 
a large number of measurements on different samples, including EuBa2Cu306+x; and have a 
considerable error both in the intensity (notice the logarithmic vertical scale) and in x, since the 
three peaks are not always as clearly distinguishable as in Fig. l5-4b and due to uncertainties 
in the determination of x; also, the ends of the curves for x < 0.1 and x > 0.9 are a guess. 
Nevertheless, the correlation between Qmax (^) j) Pix,j) from Ref. jl09j is excellent, for all 
three relaxation processes; a discrepancy can be found in the fact that Q^ax (^) -P^) is about an 
order of magnitude smaller than Q^^x (^' -PH2) , while p {x, P2) is slightly larger than p (x, PH2) . 
A possible explanation might be that the pix,j) have been calculated for T = 300 K, while 
peak P2 is measured at lower temperature, where the probability of isolated O atoms should be 
reduced in favor of chain fragments; on the other hand, in the light of discussions in the next 
Sections, I think that it is too much to look for a quantitative agreement between i^ij) 
and these p {x,j), and the comparison should remain on the qualitative level. 

Before discussing in more detail the information obtainable from the three relaxation pro- 
cesses P2, PHI and PH2, I will discuss what kind of anelastic spectra one would expect from 



105 



the information presented in Sec. |2J 



5.4.2 Elastic dipole of oxygen in the CuO^ plane 

It is possible to estimate the elastic dipole associated with an O atom in the CuOx plane by 
the X dependence of the cell parameters, and it has been noted that the elastic dipole is almost 
independent of the oxygen content and ordering |127j . The estimate can be made by calculating 
the lattice parameters ot and ao, bo of the tetragonal and orthorhombic cells in terms of the 
lattice parameter a of the ideal oxygen-free tetragonal cell and elastic dipoles A*-^^ and A*^^^ of 
O in the sites of type 1 and 2. These sites are generally labeled as 0(1) and 0(5), depending 
whether the nearest neighbor Cu atoms are in the y ox x direction respectively (with x \\ a and 
y II 6), but I will use the label 2 instead of 5. Let us indicate the components of the elastic 
dipole associated with oxygen as Ai = X^^x = A^^*, A2 = A^^ = X^yy > Ai. The remaining 
component along z is not interesting, since it does not change after a jump. The strain due to 
the occupation of the ni + n2 = x sites is 



and the cell parameters in the O and T phases, assuming n\ = xq, = Q and n 
= X'j:/2^ are: 



Eij = ni 




(5.2) 



ax = a 1 + XT- (Ai + A2) 



(5.3) 



ao = a [1 + xqXi] , 60 = a [1 + X0A2] 



(5.4) 



from which one deduces 



a = 



aTXQ - 2 (qq + bp) XT 



(5.5) 



XQ - XT 



(A2 - Ai) 



bo - ao 



(5.6) 



XoCL 




(5.7) 
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In the following table I consider the cell parameters measured by neutron diffraction in two 
tetragonal and two orthorhombic samples in Ref. |132j . and two pairs of data taken from Fig. 
5 of Ref. [TTn] . 



XT 


xo 


ax (A) 


ao (A) 


bo (A) 


Ref. 


a (A) 


(A2-A1) 


5(Ai + A2) 


0.18 


0.96 (O-I) 


3.8587 


3.8227 


3.8872 


,132 


3.8596 


0.0174 


-0.0012 


0.18 


0.78 (O-III) 


3.8587 


3.8265 


3.8875 




3.8592 


0.0203 


-0.00073 


0.25 


0.96 (O-I) 


3.8586 


3.8227 


3.8872 




3.8599 


0.0174 


-0.0013 


0.25 


0.78 (O-III) 


3.8586 


3.8265 


3.8875 




3.8591 


0.0203 


-0.00069 


0.095 


0.93 


3.8600 


3.8227 


3.8872 




3.8606 


0.0180 


-0.0016 


0.28 


0.93 


3.862 


3.8227 


3.8872 


mni 


3.8650 


0.0179 


-0.0028 



The sample with xq = 0.78 was detwinned and showing O-III superstructure (one filled chain 
every three) and the authors report negligible occupancy of 0(5) in the ortho phases; therefore 
it should be meaningful to assume nf^^ = also for such a relatively low value of x. From these 
data it can be concluded, in accordance with Ref. |127j . that (A2 — Ai) ~ 0.019 it 0.0015 and 
i (Ai + A2) ^ -0.0018 ± 0.0001. 

An important remark must be made on this type of derivation of the elastic dipole of oxygen 
in the CuO^ planes: it is valid for oxygen aggregated into chains or chain fragments, but 
not necessarily for isolated O atoms; in fact, in Sects. l5.4.6ll^^T!Hl it will be shown that 
around room temperature and below, the concentration of isolated O atoms is very small also 
in highly oxygen deficient samples. As explained in those Sections and in Sec. 15.21 the state of 
aggregated and isolated O atoms is different, due to charge transfer effects dependent on the 
chain length, and it is not obvious that the elastic dipole of an isolated O atom should be the 
same as that of an oxygen belonging to a chain. Therefore, it is not obvious that the value of 
the anisotropy AA = 0.019 found above is appropriate for evaluating the anelastic relaxation 
from oxygen hopping. 

Let us see if it is compatible with the peak intensities reported in Fig 15-41 In Appendix 
B is calculated an estimate of the relaxation strength for the Young's modulus of an isotropic 
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polycrystal; neglecting porosity, it should be 



where c is the concentration of relaxing dipoles. The Young's modulus of ceramic YBa2Cu306+x 
corrected for porosity has been estimated as = 90 GPa and 125 GPa for low and high x, 
respectively [nSlIISll- Setting E = 100 GPa, p ~ 6.3 g/cm^, vq = 175 x 10"^^ cm^ we get 

A 30 K 

A~^c. (5.9) 

If we want to compare with the height of peak PH2 we set T ~ 750 K and c ~ 0.1, since 
according to the vacancy mechanism |124j the intensity is proportional to the O vacancies in 
the O-I phase and not to x ~ 0.9; it turns out Q^ax = Y ~ ^ ^ 10~^, instead of the observed 
10~^ — 10"^ (Fig- I5-4|) . The intensity however, might be enhanced by the proximity to the 
ordering transition jl27j . by a Curie- Weiss like factor; from the YBCO phase diagram (Fig. 15-2(1 
it appears that at x ~ 0.9 the transformation to the tetragonal phase, if any, is at T > 1000 K, 
and therefore T = 750 K in the denominator should be substituted with T — Tq > 250 K, 
which does not change much the situation. It seems therefore that AA = 0.019 is too 
small for explaining the large relaxation strength we observe for peak PH2, but 
these estimates are really rough. In any case, it is difficult to draw quantitative conclusions on 
the elastic dipole of oxygen from the anelastic spectra, because, when measuring peaks PHI 
and PH2 at high temperature, one has oxygen loss during the measurement, while for peak P2 
one needs an independent estimate of the concentration c of free oxygen. 

5.4.3 Interactions among the oxygen atoms in the Bragg- Williams approxi- 
mation 

As explained in detail in Sec. 12. 4| Wipf and coworkers performed aftereffect jl35| I2f)( I127| 
experiments on YBCO and took into account the interactions among elastic dipoles in the 
Bragg- Williams approximation, analogous to the Curie- Weiss approximation above Tq. They 
showed that the anelastic aftereffect (the time dependence of strain after application of a con- 
stant stress) for x ~ 0.4 around 380 K presents an enhancement of the relaxation strength and 
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time by a factor T/ [T — Tq)~ , where Tq is the temperature below which the elastic dipoles 
associated with the O atoms start ordering themselves, namely the temperature of the T/0 
phase transformation. Around 380 K, a critical concentration Xc — 0.45 is found, where both 
the relaxation strength A (x) and relaxation time r (x) are enhanced by at least a factor 10 over 
the values found away from Xc |127j . The value Xc ~ 0.45 is somewhat higher than x ~ 0.3 — 0.4 
deduced from the 0-T line of the experimental phase diagram jl()9| I136j . 

The elastic interaction alone cannot be the origin neither of the tetragonal to orthorhombic 
transition nor of the several orthorhombic phases, which are mainly determined by short range 
interactions of electronic origin, as schematized in the ASYNNNI model. In fact, if one assumes 
that the ordering temperature is given by k^Tc = a| (l — |) [see eq. (|8.3|) ] and tries to fit the 
0-T boundary of YBCO with such Tq (x), one obtains a/k^ = 2500 — 4000 K, as shown in Fig. 
I5-5| On the other hand, Wipf |137j estimated the elastic contribution to a as a^^/k^ ~ 290 K, 
which is an order of magnitude smaller than the value needed to reproduce the observed ordering 
temperatures. 




Kin YE^JIjiIi. ^ 

Figure 5-5: Fit of the tetragonal/orthorhombic boundary with Tq (x). 

Therefore, a treatment of the interaction among the O atoms in the mean field approxima- 
tion, as presented in Sec. 12.41 is not fully consistent if the main contribution to the interaction 
parameter a is of local electronic origin; still, it is useful to reproduce important effects in the 
relaxation strength and rate of the anelastic peaks due to oxygen hopping, as shown by their 
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critical enhancement observed in after-effect experiments |127j . 




Figure 5-6: s" {uj,T) curves in logarithmic (upper) and linear (lower) scales calculated in the 
Bragg- Williams approximation with tq = 10^^'^ s, E/k = 11000 K and a'/27r, x and a as 
indicated in the legends. 

According to the phase diagram of YBCO (Fig. l5-2j) . on varying x the temperature Tq {x) of 
the T/O-I ordering transition crosses at x ~ 0.5 the temperatures of peaks PHI and PH2. Then, 
it would be interesting to evaluate the effect of this ordering transition, at least in the simple 
Bragg- Williams approximation presented in Sec. 12.41 Above Tq, the simple Curie- Weiss-like 
formula (|2.50|1 can be used, but below Tq one has to calculate the expressions p.72|) and (|2.78j) . 
or equivalently the order parameter ^ = ^ (ci — C2). In Appendix A numerical approximations 
are provided for ^ (c, T) ,which might be useful for fitting purposes. As already noted, our 
data are generally not susceptible to reliable analysis, since oxygen loss occurred during the 
measurements, therefore I simply plot in Fig. 15-61 some s" {uj,T) curves calculated according 



110 



to Eqs. H2.72() . (|2.78(1 with the approximations (|8.5(1 . H8.6() . The curves are plotted in linear 
(upper panels) and logarithmic (lower panels) scales as a function of x in (a) and (b) and as a 
function of frequency in (c) and (d); the relaxation time in the absence of interaction has been 
chosen as the typical hopping time of O, r = 10~^^ sexp (11000/T), the interaction parameter 
a/k-Q has been chosen as 3500 K, in order to have the transition in the correct temperature 
range, and also plotted as a dashed line is the reference peak with a = (no interaction) at 
Tp = 520 K (for 1 kHz). The ordering transition on cooling is clearly visible as a narrow peak, 
when Tq > Tp and as a drop of intensity when Tq < Tp. Although the anelastic spectra we 
measured for x < 0.5 contain some narrow peaks clearly due to oxygen ordering transitions 
|138j . there is no clear relationship with the curves calculated in Fig. 15-61 and it is evident that 
a more complete treatment of the short range interactions, hke in the ASYNNNI 
model, is necessary to describe the high-temperature anelastic spectra of YBCO. 

5.4.4 Hopping in the O-I phase - peak PH2 

The information obtainable on PH2 from the present measurements is only limited to the 
observation that it is a process occurring in the O-I phase, distinct from PHI; in fact, the 
experiments are made in vacuum and rapid oxygen loss occurs during the measurement, affecting 
both shape and intensity of the peak. More reliable measurements of this process have been 
made by other authors with the pendulum jl24j . which allows some O2 pressure to be maintained 
around the sample and moreover lowers the peak temperature. The mechanism devised for PH2 
in Fig. 15-41 is the one requiring a neighboring vacancy, as suggested by Xie et al. jl24j . on the 
basis of fact that its intensity vanishes for x ^ 1. That analysis, however, is not fully consistent, 
since it assumes an asymmetry between the energies in sites 0(1) and 0(5), IS.E (x) = E^ — Ei, 
which is maximum for x = 1 and decreases to when the sample becomes tetragonal. Such 
an asymmetry energy was taken into account for reproducing the decrease of the apparent 
activation energy with decreasing x, but it was not considered that it would also reduce the 
intensity of the peak, due to the depopulation factor in Eq. (j2.15|) . thereby accounting for a 
consistent part of the reduction of the peak intensity when x ^ 1. 

Among other measurements of peak PH2 |1251 113{?1 1126j . I would like to signal isothermal 
spectra measured with the pendulum |140j that confirmed the presence of two peaks, the one at 
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higher temperature is associated with the orthorhombic phase and should therefore correspond 
to PH2, while the other one to PHI. This confirmation is important, since it is based on 
stationary and reproducible spectra. 

5.4.5 Other jumps of aggregated oxygen and oxygen ordering - peak PHI 

Figure 15- 7b . (from Ref. JISOJ) shows peak PHI measured in an EuBa2Cu306+x; which has 
anelastic spectra and all the physical properties very similar to YBa2Cu306+x- The peak with 
higher intensity is for x ~ 0.5, and proceeding with the out gassing treatments it decreases in 
intensity and shifts to higher temperature; the effective activation energy, estimated from the 
shift of the peak at higher frequency, also increases from 1.1 to 1.3 eV. While the decrease in 
intensity is easy to understand after the considerations of Sec. 15.4.11 the increase in activation 
energy has no obvious explanation. 
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Figure 5-7: (a) Peak PHI measured in EuBa2Cu306+x starting from x ~ 0.5; with decreasing 
x,the peak decreases in intesity and shifts to higher temperature, (b) After quenching to room 
temperature from 990 K, peak PHI is strongly depressed, but reappears with aging at high 
temperature. 

Figure IS^Tb shows another interesting effect, which excludes the hopping of isolated O 
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atoms as possible mechanism for peak PHI. Curve ^1 has been measured on heating 
after quenching from 990 K to room temperature in few seconds, by pouring water on the 
quartz tube where the sample was heated. After such a quenching the concentration of isolated 
O atoms should be maximum, but peak PHI is strongly depressed, and instead a small peak 
appears near 480 K. The latter does not shift in temperature with increasing frequency and is 
therefore signalling the occurrence of an ordering transition near 480 K, like the frequency 
independent peak at 560 K in the curves of Fig. l5-6t l: it is difficult to say whether the ordering 
occurring above 480 K is simply a lengthening of the Cu-0 chains or the formation of O-II 
domains. It is clear however, that peak PHI is due oxygen jumps within such long chains or 
domains and not to isolated O atoms or even pairs or very short chains, which are promptly 
formed. In fact, the peak develops with thermal cycling above 500 K (curves 2 and 3), and the 
increase of intensity cannot be attributed to oxygen uptake, since the measurements are made 
in high vacuum. 

We also observed additional frequency independent anomalies up to 580 K jl38j . which ap- 
pear only on heating and not on cooling (as for curve 3 in Fig. l5-7b ). but still there is no better 
explanation than some types of ordering transitions. This complex phenomenology demon- 
strates that situation is even more complex than it appears from the diffraction measurements 
of oxygen ordering. 

5.4.6 Isolated oxygen atoms - peak P2 

There is extensive evidence that the effective barrier for the diffusive jumps of the O atoms in 
the CuO^. planes is of the order of 1 eV, as discussed in the previous Sections. It should be 
noted, however, that none of the high temperature (T) peaks presents the dependence on 
X expected from the hopping of a concentration c = x/2 of an interstitial species (if one regards 
oxygen in the CuOx plane as an interstitial species, and takes into account that there are 2 
sites available for each Cu atom). As a matter of fact, none of the behaviors expected from an 
interstitial solution is observed; e.g. the proportionality of the relaxation strength to c in the 
high dilution limit, or to c (1 — c) if one takes into account the filling of a finite number of sites 
jl3j . or to or some other power of c if interstitial complexes contribute. This is due to the 
strong interaction between the O atoms, and the difficulty of probing the two limits x ^ and 
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X ^ 1. The interaction among the O atoms is so strong to prevent the fiUing of the planes above 
c = 0.5 (or X = 1). That concentration corresponds to the O-I phase, with completely filled 
chains, and can be approached but probably never obtained. We attempted at reaching the 
high dilution limit x ^ 0, but this also cannot be done, due to the decomposition of YBCO in 
vacuum at high temperature into the more stable oxides. Bormann and Nolting jl41j found the 
decomposition limit p{T) below which YBa2Cu306+a; transforms into the more stable oxides 
CU2O, BaCu02 and Y2BaCu05 (the latter easily recognizable from the green color), and found 
that the lowest attainable oxygen content is YBa2Cu306.o5. 

The first anelastic measurements on highly outgassed YBCO |128( I142j showed a rather 
unexpected result: all the relaxation processes at high temperature were completely suppressed 
jl42j and a new intense peak, labeled P2, appeared around 60 K, close to a Debye relaxation 
with an activation energy of 0.11 eV |128j (see also Fig. I5-4|) . Since oxygen was certainly 
present in the CuOx planes at the level of several molar percent also after outgassing, the most 
obvious conclusion was that P2 is due to hopping of isolated O atoms. This is striking, since 
the barrier for jumps of isolated O atoms would be 10 times lower than that for 
jumps of oxygen within chains, and also in any other known oxide; such a low barrier is 
rather comparable that for interstitial H in bcc metals. In fact, the same relaxation process 
has been observed also in NQR experiments and attributed to polaron hopping jl43j . Even 
accepting that the O atoms in the nearly empty CuOx planes overcome a barrier as little as 
0.1 eV to perform a diffusive jump, some difficulties arise in conciliating such a high mobility 
with the long time required for oxygen to reach an equilibrium configuration; this led de Brion 
et al. [mj to assign P2 to hopping of isolated O atoms trapped by some defect, e.g. a vacancy 
in the apical oxygen sublattice. It has later been argued |145j how peak P2 is more likely to 
be due to oxygen rather than small polaron hopping, considering that its activation energy is 
independent on concentration, and how it is possible that the extremely fast mobility of isolated 
O atoms may coexist with a slow kinetics for reaching equilibrium. 

Part of the following results appeared in Ref. |145j . We approached the x ^ limit 
by heating the samples in a vacuum of ~ 10^^ mbar at progressively higher temperatures 
(up to 1000 K) for ~ 2 h. The last outgassing treatment at 1000 K was slightly beyond 
the decomposition limit, and traces of superficial green phase had to be removed with emery 
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paper. The oxygen content after the final outgassing treatment was estimated from the c lattice 
parameter by X-ray diffraction, and was found x = 0.05 it 0.05, in agreement with the results 
on the YBCO stability |141j . Figure [5^^ shows the evolution of peak P2, starting with x > 0.4 
(curve 1); where the peak is practically absent, and peak P3 at 90 K (8 kHz) is still visible; the 
peak appears after outgassing up to 810 K (curve 2), with x ~ 0.3, increases after 2 h in vacuum 
at 980 K (curve 3), and starts decreasing with further outgassing (curve 4 after additional 2 h 
at 1000 K) or after long aging at room temperature (curve 5 after 13 months). 




Figure 5-8: (a) development of peak P2 with outgassing and aging at room temperature; (b) 
fits of peak P2 measured at two frequencies in the states corresponding to curves 2 and 5 in 
(a). 

The characteristics of peak P2 other than the intensity are practically independent of x or 
sample history, indicating relaxation of some simple point defect. Figure shows fits to the 
peak measured at two vibration frequencies (0.7 and 9 kHz) in the states with the highest and 
lowest intensity (curves 3 and 5 of Fig. l5-8b '). For peak P2 the Cole-Cole expression 1)2. 45() 
was found to provide a slightly better interpolation than the Fuoss-Kirkwood one, and it has 
been allowed for relaxation between slightly inequivalent sites [Eqs. (|2.15|) and ()2.22b|l with 
/S.E 7^ 0], in order to reproduce the fact that the peak intensity is slightly larger at higher 
temperature, instead of exhibiting 1/T behavior (see the discussion at the end of Sec. I2.().lj) . 
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Therefore peak P2 has been interpolated with 
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plus a linear background and a small contribution from peak P3 around 90 K. The peak may 
be fitted with tq = 1.5 x IQ-^^ g ^ E/ku = 1290 K, a = 0.84 and AE ~ 130 K in all cases, 
except for the final state (curve 5), where the slightly different values tq = 4 x 10^^^ s and 
E/k-Q = 1370 K give a better interpolation. The slight shift of the peak temperature from 
curve 2 to curve 5 is due to changes of the resonance frequency of the sample. It can therefore 
be concluded that the relaxation parameters of peak P2 are independent on x and 
sample conditions within experimental error. 

5.4.7 Isolated and aggregated oxygen atoms 

The first question to answer is how is it possible that the oxygen jumps for x < 0.3 occur 
over a barrier about 10 times smaller than for x > 0.3. The proposed explanation is that 
there are two substantially different types of jumps: fast jumps of isolated O atoms, and slow 
jumps, occurring at high values of x, that involve breaking and formation of chemical bonds 
with a neighboring O atom. Therefore, the real barrier for oxygen hopping in the CuO^ 
planes would be the one of peak P2, 0.11 eV, while the additional 0.9-1.0 eV are 
due to the formation of chemical bonds and, as we shall see later, to the electrostatic 
interaction between neighboring O atoms. This hypothesis is reasonable in view of the 
model for the hole states and charge transfer from chains to planes presented in Sec. 15.21 where 
it has been pointed out that: i) jumps of isolated O atoms do not affect the electronic state; 
ii) the formation (or dissolution) of a pair of nearest neighbor O atoms involves the change of 
a hole of Cu'^^ character into O^^, with a higher energy (but lower than having two holes on 
a same Cu^^ atom); in) jumps to or from longer chains involve charge transfer between chain 
and Cu02 plane. 

Figure 15-91 presents a sketch of these processes. From these considerations it appears that 



116 



Figure 5-9: Effect of different types of O jumps in the CuO^^ plane on the holes. From left to 
right: 1) no change for an isolated O atoms; 2) a Cu^"^ hole is converted into a O^^ hole with 
higher energy; 3) after the formation of a longer chain, a hole is pushed to the Cu02 planes. 



only the jumps of isolated O atoms involve only a simple reorientation of the elastic dipole, as 
in the Snoek effect, while all other types of jumps, besides a higher effective barrier, should 
also present different initial and final energies, and possibly elastic dipoles. The energies for 
creating a hole or changing its state are of the order of several tenths of eV, which means 
thousands of kelvin in the temperature scale. From Sec. 12.0.11 it appears that in such cases 
the depopulation factor makes the relaxation strength strongly increasing with temperature, 
and this should be easily noted as a peak whose intensity is higher when measured at higher 
frequency, instead of being smaller by the 1/T factor. The anelastic data at high temperature 
collected by us do not evidence clearly such a behavior, even though the oxygen loss occurring 
during the measurement in vacuum does not make this observation very reliable. 



Figure 5-10: Jumps of an O atom between terminals of two perpendicular chains of similar 
length leaves the hole state unaffected. 
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A possible explanation for this fact would be that the anelastic relaxation observable 
at high temperature is due to those jumps between Cu-O chains that do not change 
the overall hole state, and therefore whose initial and final electronic energies are the same. 
Figure 15-101 presents the example of a jump between the ends of two perpendicular chains. 
If both chains have n < 5, no hole will be transferred to or from the planes after the jump, 
at least according to Uimin's model |118| IllJ] : also, if their length is similar, then the total 
count of the various types of holes will be the same before and after the jump. Jumps of 
this type, although involve a high activation energy for the temporary change of the electronic 
configuration, do have the same initial and final energy, and therefore the related relaxation 
strength is proportional to 1/T, as usual. It should be noted that this type of jumps is probably 
involved in the motion of twins. Considering that the hopping energy or bandwidth of the holes 
is ~ 0.4 eV |117j . it is likely that the difference between initial and final energy of a jump that 
changes the electronic state is of the order of tenths of eV; already setting IS.E = 0.2 eV in 
Eq. ()2.15|) the depopulation factor cosh~^ {/S.E/2k^T) reduces the relaxation strength by 25 
times at 500 K, a temperature where the peaks due to the chain rearrangements are observed. 
These considerations have been overlooked up to now in all anelastic studies of the oxygen 
jumps, which implicitly assume the equivalence between the electronic energy before and after 
the oxygen jump. A reliable analysis of such processes should take into account the energies 
involved in the rearrangements of the chains and the statistical weights of the various chain 
configurations; this is outside the scope of the present investigation, in view of the impossibility 
of measuring a stable anelastic spectrum in vacuum at high temperature. 

5.4.8 Slow achievement of equilibrium and expected concentration of iso- 
lated oxygen atoms 

The mobility of the isolated O atoms deduced from peak P2 is extremely fast, and corresponds 
to 10^^ jumps per second at room temperature. With such high a mobility one expects an 
almost instantaneous achievement of the equilibrium configuration of the O atoms, contrary to 
observations. In fact, oxygen ordering at room temperature may proceed for months both at 
high and low values of x. A well known example is the fact that if a sample with relatively 
low X is quenched to room temperature, its superconducting temperature is found to increase 



118 



of over 10 K with aging at room temperature over periods of days |146| : this is attributed to 
slow reordering of oxygen at room temperature, with consequent change of the concentration 
of holes injected in the Cu02 planes and change of Tc. Another example is the slow decrease 
of the intensity of peak P2 with aging at room temperature, which can be attributed to a 
slow decrease of the number of the isolated O atoms at room temperature and below. The 
discrepancy between the supposed fast hopping rate of the isolated O atoms and the slow 
reordering kinetics has been put in evidence by de Brion et al. |144j with the following simple 
argument. In the low x limit and assuming that only O pairs are formed, the authors assumed 
that the time for reaching equilibrium, teqj is the time for forming pairs starting from a random 
distribution of isolated (free) O atoms; therefore they obtain teq ~ ''"//^^i where r/ is the mean 
hopping time of the free O atoms and x their concentration. The extrapolation of r from peak 
P2 to room temperature yields Tf ~ 10^^^ s, which would require x ~ 10^^ for obtaining 
teq ~ 1 day and x ~ 5 x 10^^^ for extending teq to 1 year, as we observed. Such concentrations 
are of course unrealistic, being x > 0.1, and this fact induced de Brion et al. to attribute peak 
P2 to hopping of oxygen trapped by some defect, or to the reorientation of O pairs, leaving 
Tf ^ Ta = TQBxp (Ea/k-QT), the usual hopping time for aggregated oxygen with a barrier of 
-Ea ~ 1 eV. A flaw in this argument is the assumption that equilibrium is reached once all the O 
atoms have had the opportunity of forming a pair or chain fragment for the first time. This is 
not true, since the process of reaching equilibrium will require further steps of dissociation and 
aggregation of the chain fragments in the intermediate metastable configurations, and these 
require jumps with the characteristic time of the aggregated oxygen; in conclusion, the limit 
steps are those with r^. 

The picture of only two hopping times, tj = tq exp [Ef /k-QT) with Ef = 0.11 eV for hopping 
of isolated O atoms and = tq exp {Ea/k-^T) with Ea — 1 eV for leaving a chain fragment or 
diffusing between chains (see Fig. 15- lib ), is nevertheless inadequate. In fact, the rate equation 
for the equilibrium between free and aggregated O, in the simplest form and neglecting any 
geometrical factor, would be 



where x is the probability that the jump of the free O atom joins another O atom. By posing 
the condition of stationarity and recalling that Ca + Cf = x and Ta ^ tj, one obtains for the 



dCf Ca Cf 
dt Ta Tf 
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equilibrium concentration of free oxygen 
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XTa + Tf r< 
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(5.13) 



which is of the order of 10^^^ at room temperature with a binding energy = Eg — Ef ^ 0.9 
eV. We find again that the concentration of free O would be negligible at room temperature 
and below, while peak P2 requires cj ~ 10^^ below 100 K. 





Figure 5-11: Potential felt by an isolated O atom (3) on approaching another O atoms (1) in 
site 2. In (a) only the binding energy Ef, of the pair is taken into account; in (b) also the 
electrostatic repulsion is considered. 

The explanation for this apparent inconsistency also explains why peak P2 is not observed 
at X > 0.4. In fact, it is proposed that the electrostatic repulsion between the ions, 
which is effectively screened in the conducting phase by the mobile holes, becomes 
strong in the semiconducting phase, with x < 0.4; the result is an increase of the saddle 
point, as shown in Fig. 15-llb . The electrostatic energy of the two O atoms in the saddle 
point configuration between second and first neighboring positions is x 0.6 eV, where q is 
the effective oxygen charge, q ~ —2, without taking into account any screening. It is an energy 
2.5 larger than Ea, and this demonstrates the importance of screening in determining the short 
range interaction between the O atoms. If the saddle point between sites 3 and 2 is larger than 
that between two free sites 3 and 4, than we have to introduce a "trapping" activation energy Et 
for the pair formation, and Ef, becomes Eg—Et, where Et is not known but can be a large fraction 
of Ea- The new picture is that of a generally metastable concentration c/ of free O atoms that is 
small at room temperature and below, but larger than the equilibrium cj, unless very prolonged 
aging over years occurs after the last excursion to high temperature. These free O atoms jump 
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with an extremely fast mean time Tf and therefore approach continuously other O atoms 
and chain fragments; the mean time for aggregating with other O, however, is the 
much longer Tf = tq exp (Ef / k-^T) , at least in the semiconducting regions. This fact 
maintains the concentration of excess free O rather stable, making perfectly reasonable the 
observation of peak P2. In this respect, it should also be noted that the concentration of free O 
atoms at the peak temperature Tp ~ 65 K is not Cf (Tp), but much closer to the starting room 
temperature value Cj (RT) > cj(RT). In fact, assuming that the starting metastable room 
temperature value Cj (RT) is determined by the longer hopping time Tq, during cooling, c/ (T) 
will be approached with the faster characteristic time r^, but for a measuring run with a cooling 
rate of the order of 1 K/min, c/ will remain frozen when Tt reaches ~ 10 min, therefore above 
cj (T*), where Tt (T*) = 10 min. I do not know of any experimental indication of Ef,, but a 
range of values sufficient to freeze c/ ~ 10""^ necessary for observing peak P2 can be estimated 
setting 10"^ ~ c/ (T*) ~ exp (^^^ft^ and 10 min= Tt (T*) = tq exp (^^^^ with the constraint 
Eb + Et = Ea^l eV. It results that for 0.2 eV < < 0.3 eV one freezes 10'^ < c/ (T*) < lO""^ 
around 250 K < T* < 280 K. 

At last, we note that in the proposed picture, at room temperature and below, the rearrange- 
ment of the aggregated O atoms is extremely slow (tq (290 K) ~ 6 h) and a minor concentration 
of highly mobile O atoms remains, that are bound to diffuse in the semiconducting regions free 
of O chains. Since their mean hopping time tj is much faster than the mean time for joining 
other O atoms, Tt, being r//rj = exp [{Ef — Et) /k-sT] > 10^, they can be considered as effec- 
tively free from the strong interactions giving rise to the ordered phase and described in Sects. 
12.41 and 15.4.31 In fact such interactions involve changes of the great majority of aggregated O, 
which is frozen on the time scale of tj . It is therefore consistent to ignore any ordering temper- 
ature Tq in eq. 1)5. 10() . The introduction of a mean asymmetry AE ~ 130 K between the free 
sites is justified by the fact that the regions free of O pairs and chains, where the isolated O 
atoms diffuse, are certainly small and their site energies are strongly perturbed by strain due 
to the surrounding aggregated O. This effect is particularly evident for defect relaxation with 
small activation energy, giving rise in this case to AE/E = 0.1, but is unnoticed for the high 
temperature relaxation with AE/E = 0.01. The magnitude of AE is typical of the site energy 
shifts due to strain interactions, as discussed in the end of Sec. 12.0.11 
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5.5 Motion of oxygen in off-center positions of the Cu-O chains 



The low temperature anelastic spectrum of well oxygenated YBCO is shown in Fig. 15-121 and 
presents three peaks, which we label PI, P3 and P4, with activation energies of 0.11, 0.16 and 
0.19 eV respectively |147j . Peak PI will be dealt with in the next section. 
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Figure 5-12: Low temperature elastic energy loss of YBCO measured at 1.8 and 23 kHz. The 
continuous lines are a fit as explained in the text with the dashed lines as background. 

The continuous lines in Fig. 15-121 are a fit with the Fuoss-Kirkwood expression 1)2.42(1 for 
peaks P3 and P4 and ()2.48|) for PI plus a background of the form tanh [(T — Tq) /w] (dashed 
lines). Again, in order to reproduce the intensities of the peaks at both frequencies with 
the same set of parameters it is necessary to introduce the possibility that relaxation occurs 
between energetically non-equivalent states. This requirement is normal in such a disordered 
system like YBCO at temperatures below 100 K (see the end of Sec. I2.fl.l|) . It is therefore 
necessary to introduce a correction of the type cosh~^ {/S.E/2k^T) to the relaxation strength 
and cosh""*^ {AE /Ik^T) to the relaxation time (Eqs. ()2.15() and ()2.22b|) ). The parameters used 
for the fit are reported in the table. 
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The pre-exponential factors tq are all compatible with point defects, while the activation 
energies of P3 and P4, 0.15 and 0.21 eV, are very close to those deduced from the condition 
LOT = 1 at the maxima of the peaks measured by several authors in a broad range of frequencies, 
0.16 and 0.19 respectively |147j . The peaks are also considerably broader than pure Debye 
relaxations, and the asymmetry parameters AE are within 10% of the activation energy; they 
do not affect much the other parameters but definitely improve the fit. The evolution of peaks 
P3 and P4 with doping is shown in Fig. 15-131 after an outgassing treatment to x ~ 0.5, as 
deduced from resistivity (absence of superconductivity down to 70 K and increase of resistivity 
by 20 times at room temperature), peak P4 was almost suppressed, while P3 slightly enhanced. 
The broad peak at 220 K seems to be connected with some type of oxygen ordering, as briefly 
discussed in Sec. 15.71 
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Figure 5-13: Elastic energy loss and resistance normalized to room temperature of an YBCO 
sample with x ~ 0.95 (curves 1) and after outgassing to x ~ 0.5 (curves 2-4) [from |147j ]. 

The existence of relaxation processes with low activation energy and disappearing when the 
Cu-0 chains become shorter and more disordered is naturally put in relation with the fact that 
the O atoms in the Cu-O chains actually occupy positions that are 0.15 A ofF-centre 
along the a axis. This fact is deduced from the anomalously large Debye- Waller factor of 
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these atoms in the a direction, from refinements of neutron jl48j and X-ray jl49j diffraction 
spectra. Such anomalous thermal factors are relatively common in perovskite-related materials 
and are indication of off-centre occupation. In addition, a Mossbauer experiment indicated that 
such zig-zag configurations are dynamic with a correlation time > 10^^ s |15U| . The doping 
dependence in Fig. 15-131 indicates that P4 rather than P3 should be associated with these 
off-center jumps, since the outgassing treatment at least halved the number of the off-center 
O atoms in the chains, and P3 did not decrease at all. On the other hand, the reduction of 
the intensity of P4 by much more than a factor ^, might indicate that only the O atoms in 
relatively long and ordered chains are free to jump between off-center position; in 
fact, in a phase of disordered chain fragments, most of the pairs of off-center sites would be 
highly inequivalent, and the O atoms would remain in the sites of smaller energy. 

A model has therefore been developed of how the jumps of the O atoms between the off- 
centre positions might produce anelastic relaxation jl47j . The reasoning goes along the lines 
of Sec. l2.U.lll?ir^ with the additional observation that the relaxing units are pairs of 
neighboring O atoms. In fact, an isolated O atom hopping between the two off-center 
positions would not give rise to anelastic relaxation, since its hopping would simply cause an 
inversion of the off-centre defect, and the centrosymmetric elastic dipole would not change. 
Instead, if we consider pairs of O atoms or chain segments, than we can distinguish between 
pairs of atoms that are on a same side of the chain or on both sides. If we call A and B the 
sites on either side of the chain, we can distinguish between chain segments of type 1, in AA 
or BB sites and having elastic dipole Ai, and segments of type 2 in AB or BA sites with elastic 
dipole A2. An oxygen jump like that on the left in Fig. 15-141 will change the contribution to 
the elastic strain by Ae™ = 2 (Ai — A2) = 2AA, while a jump like that on the right hand will 
let e^"^ unaffected. 

Also the profile of the potential energy is different in the two cases, since in the first case 
relaxation occurs between states differing in local strain 2AA and in energy by 2AE ^ 0, while 
in the second case also the energy must be the same by symmetry (simple inversion of the 
configuration). Therefore, the origin of the consistent asymmetry AE/kB — 210 K found for 
peak P4, may be indicative of the energy inequivalence between the two types of chain segments, 
rather than disorder. 
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Figure 5-14: Possible types of O jumps between the off-center positions A and B of the zig-zag 
Cu-0 chains. 

5.5.1 Possible (anti)ferroelastic and (anti) ferroelectric ordering of the off- 
center oxygen atoms 

Finally, we point out that (anti)ferroelastic correlations may arise between the off-centre O 
atoms, which also involve electric dipoles along the a axis, associated with atoms that 
are out-of-axis with respect to the row of Cu^^ atoms along b. Therefore, such correlations 
would also be of (anti) ferroelectric nature, although the electric field from these dipoles is 
effectively screened in the metallic or superconducting environment, and is therefore hardly 
detectable in experiments. Nonetheless, interactions between neighboring O atoms along the 
chain or even across adjacent chains cannot be excluded and might even result in short-range 
(anti) ferroelectric and (anti)ferroelastic domains. Indeed, the (T) curves shown here are 
measured on cooling, but the measurements on heating after some aging below 90 K exhibit 
clear phase transformations between 120 and 170 K jl51j . indicated with PTl and PT2 in Fig. 
15-151 In Ref. jl51j we discussed how these anomalies in both dissipation and modulus might 
arise from ferroelectric and/or antiferroelectric domains that are formed during low temperature 
aging along the Cu-0 chains. We also proposed a possible evolution of peak P4 with decreasing 
X toward a slower and more correlated dynamics, as suggested by a shift to higher temperature 
and broadening of the peak in some measurements of outgassed samples. Such an increase of 
the correlation between the chain segments would be due to the decreased electric screening 
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when going into the semiconducting state. 




Figure 5-15: Elastic energy loss and resonance frequency ( oc VE) of YBa2Cu306.35 measured 
on cooling and subsequent heating after 20 h at 80 K. 

In general, one should be very cautious in considering phase transitions like these, since 
during several hours of low temperature aging in an ordinary vacuum system without partic- 
ular precautions a porous sample may adsorb residual gases and then present their 
solid/liquid transitions on heating. This fact is exemplified by the measurement made on a 
sample of alumina with ~ 50% porosity (kindly supplied by E. Roncari at CNR-ISTEC, Faenza, 
Italy) presented in Fig. 15-161 without using the adsorption pump during the experiment. 

The measurement on heating present three anomalies around 35 — 45 K, 230 K and 268 K, 
attributable to O2/N2, CO2 and II2O which were adsorbed onto the porous sample, especially 
at LHe temperature, and which pass from the solid to the liquid state and finally desorb on 
heating. It was also possible to monitor an increase of pressure in the system in correspondence 
with these anomalies. The anomaly of Fig. 15-151 however, should be intrinsic of YBCO, since 
it occurs at a temperature where nothing is observed in the alumina dummy sample. 
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Figure 5-16: Elastic energy loss coefficient and Young's modulus of porous alumina measured 
on cooling and subsequent heating. The anomalies on heating correspond the solid/liquid 
transitions of residual gases adsorbed by the sample. 

5.6 Bipolaron reorientation in the overdoped state 

The anelastic relaxation process with the lowest activation energy, E = 0.11 eV, is labelled PI 
and is observed only in the overdoped state: it is absent with x < 0.85 and its intensity rises 
very fast on approaching full oxygenation 21 . This is a really noticeable feature, since most of 
the physical properties of YBCO are practically constant in the range 0.85 < x < 1, including 
the superconducting temperature Tq (x), which presents a maximum at the so-called optimal 
doping X ~ 0.93, but whose variation remains within 4% in that range, as shown in Fig. 15-181 
fromRef. [T^ . 

Figure 15-171 presents the anelastic spectrum of YBCO below 180 K, measured at 22 kHz 
on the same sample in several different states from optimally doped to overdoped, obtained in 
different ways over a period of 6 years. The peaks at 100 K and 120 K (P3 and P4 discussed 
above) remain remarkably stable, and only the intensity of PI at 50 K changes by at least one 
order of magnitude. Curve 1 is the as-prepared state, where x was estimated from the lattice 
parameters as xi ~ 0.93 21 . The curve is closely similar to that measured on an identical 
sample 6 years earlier and therefore demonstrates the stability of samples and anelastic 
spectra. Notice that, in view of the near constancy of the lattice parameters for x > 0.9, 
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Figure 5-17: Anelastic spectrum of YBa2Cu306+a; at various values of x in the overdoped 
region, with peaks PI, P3 and P4, starting from low temperature. 

the above estimate of x has great uncertainty, and should be increased. In fact, after a mild 
annealing at 470 K for 20 h in a vacuum better than 10^^ mbar, some oxygen loss certainly 
occurred, and Tq increased from 90.3 K to 91.9 K, indicating that the sample passed from the 
overdoped to the optimally doped state. 

With the help of the Tq {x) curve from literature |152j . we can set xi = 0.96 and X2 = 0.93, 
as shown in Fig. 15-181 After such a small decrease of the oxygen content, the intensity of 
PI decreased by 3-4 times (curve 2); after 15 days the intensity was further reduced (curve 
3), possibly due to some oxygen reordering. The sample was then oxygenated as fully as 
possible in the UHV system described in Sec. 13.51 by heating to 600 °C, introducing a static 
atmosphere of 1250 mbar O2 and slowly cooling to room temperature (where the O2 pressure 
became 740 mbar); the cooling rate was 0.2 K/min except between 490 and 300 °C, where it was 
reduced to 0.1 K/min. After this full oxygenation the intensity of PI become about twice the 
original one (curve 4), and remained stable during subsequent runs (curve 5). The sample has 
been finally outgassed in the UHV system and equilibrated at 460 °C with a known amount of 
O2, resulting in a final pressure of 9.6 torr O2, which should result in 2; = 0.85 according to the 
phase diagram of Ref. |153j . From the amount of gas absorbed it was estimated x = 0.89 ±0.04, 
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and therefore it can be set xq = 0.87; at this oxygen content no trace of peak PI was present. 
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Figure 5-18: Maximum of Tq (x) inYBa2Cu306+x (after '15^). Superimposed in red are the 
intensities of peak PI corresponding to the curves in Fig. 15-171 at values of x estimated with 
the help of the Tq (x) curve. 

Considering that all the structural parameters remain practically constant within 
the range of x in which PI develops, it is impossible to associate its mechanism with 
structural defects, impurities, defects in the order of the Cu-0 chains, or off-center atoms, as 
also discussed in detail in [2J|. The only quantity that to our knowledge starts increasing in 
the same doping range is the concentration of holes of character p^, deduced from the peaks 
in the X-ray absorption spectra (XAS) of the Cu-K edge |115j . and shown in Fig. 15-191 

These holes, having symmetry p^, must reside on the apical O atoms (see Fig. 15-11) . which 
are the only atoms with consistent hybridization of the 2pz orbitals. The correspondence 
between the growth of the concentration of these holes for x > 0.8 and the growth of peak PI is 
demonstrated in Fig. 15-191 where vQmax of PI is also plotted. In Fig. 15-201 the pz orbitals are 
sketched together with those of O 2px and O 2py character, hybridized with the Cu 3d orbitals in 
the Cu02 planes. Regarding the possible anelastic relaxation processes arising from hopping of 
the holes among these orbitals, the hopping between px and py would of course produce anelastic 
relaxation, but these holes are extremely mobile, in a band giving rise to superconductivity. 
The relaxation rate characterizing the equilibrium between the px and py populations is 
certainly much faster than the thermally activated rate over a barrier of 0.11 eV found for peak 




129 




Figure 5-19: Concentration of holes with prevalent O 2p character, measured by XAS (from 
jllSj ). The concentration oi Pz holes is evidenced in red, and is in good agreement with yQnmx 
of peak PI (right ordinate axis). 

PI. Also the process of change of a hole from Pxjy the Cu02 plane to apical pz can hardly 
explain peak PI, since, as already discussed for peak P2 in Sec. 15.4.71 different electronic states 
are expected to differ in energy in the order of 0.1 eV or more; this would lead to a depression 
of the relaxation strength by a factor cosh^ [IS.E /2k^T) which, already for /S.E = 0.1 eV would 
be 3 X 10^^. This means that the relaxation process corresponding to the Px/y ^ Pz exchange 
is unobservable. 

The hopping of a hole among the pz orbitals of different apical O atoms might well have a 
correlation time like that of PI, since these orbitals are too much separated to form a band, 
and two intermediate p^/y ^ Pz exchanges are necessary. In this case initial and final energy 
are the same and there is no problem of depression of the relaxation strength; rather, anelastic 
relaxation is impossible because of the equivalence of all the Pz orbitals also under the appli- 
cation of any stress. On the other hand, if two holes on neighboring apical O atoms form a 
stable pair, a bipolaron, the resulting elastic dipole has orthorhombic symmetry and reori- 
ents itself together with the hole pair. It is therefore proposed that peak PI is due to the 
reorientation of pairs of holes, or small bipolarons, on neighboring apical O atoms. 

The intensity of PI is then expected to be proportional to the square of the 
concentration of pz holes, and this is verified in Fig. 15-191 where the square root of the 
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Figure 5-20: Px hole in the Cu02 plane and Pz hole on an apical O in YBa2Cu307. 

intensity of PI (closed circles and right ordinate) corresponds very well with the number of 
Pz holes. The activation energy of 0.11 eV would be connected wit the p^/y <-> Pz exchange 
process, necessary for the bipolaron reorientation. 

It is noticeable that the position in temperature and shape of peak PI do not change 
with doping within experimental error. Figure 15-211 presents the experimental data of 
Fig 15-171 after subtraction of a linear background and normalization to the peak intensity, and 
all data fall on the same curves. As an additional check of the reliability of the background 
subtraction, also the difference between curves 4 and 2 is found to fall on the same curves, and 
it does not depend on the choice of the background. 

The continuous lines are a fit with the Jonscher expression, eq. 1)2. 49(1 . with ti = T2 = 
Toexp{E/kBT) , To = 5 X lO^^^ s, E/ks = 945 K, a = 0.33 and /3 = 0.41. The values 
of a and (3 are definitely smaller than 1, indicating considerable broadening with respect to 
the monodispersive Debye case. This fact is not unexpected, especially for a low-temperature 
relaxation process in a rather complex and disordered system like YBCO, where random energy 
shifts due to interaction with various defects, including disorder in the Cu-0 chains, may easily 
be of the order of 0.1 — 0.2 times the activation energy. In the present case, however, the main 
source of disorder are O vacancies in the Cu-0 chains, and their concentration 1 — x decreases of 
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Figure 5-21: Peak PI after linear background subtraction and normalization to the peak in- 
tensity; the symbols are the same of Fig. 15-171 . The continuous line is the fit as explained in 
text. 

about 4 times from curve 1 to curve 4 without resulting in any narrowing of the peak. Neither 
can the dynamic interaction among the pz holes be the source of broadening, since also their 
concentration changes by more than 4 times. The spectrum of relaxation times should 
therefore be ascribed to the interaction with phonons and the other charge carriers, 
which remain practically constant within the range 0.85 < x < 1 (see also Fig. I5-19|) . 

5.7 The anomaly near 240 K 

In YBCO there is a relaxation peak in the absorption around 240 K with the characteristics of 
glassy dynamics, accompanied by a hysteresis in the Young's modulus below the same temper- 
ature |154[ 11551 1156j . This anomaly can easily be masked by the freezing transition at ~ 220 K 
of the pump oils adsorbed in porous samples from the measurement vacuum system, as first 
noticed by Gzowski [^J; in addition, the hysteresis in the elastic moduli may be affected by 
the sample microstructure jl57| I158j , since the anisotropic thermal expansion may cause severe 
stresses from grain to grain during thermal cycling. Nonetheless, an intrinsic effect exists, has a 
strong dependence on the oxygen content, and has been confirmed by several later works. Re- 
garding the possible mechanism, there are different proposals, like an oxygen ordering transition 
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jl54l 11551 1133j . a ferroelectric or ferromagnetic transition |159j or some polaronic mechanism 
[THnilTHT] . 

5.8 Anomalies at the superconducting temperature 

There are various reports of elastic anomalies at the superconducting transition 11621116,^11641 . 
also by means of non traditional techniques jl65j and we also found that sometimes anomalies 
in the Young's modulus appear near Tc. I will not try to analyze such anomalies, considered 
the multitude of effects affecting the elastic moduli of ceramic samples (see also the previous 
Section), and, as observed by Mizubayashi et al. |166j . the presence of the modulus defects 
associated with peaks P3 and P4 in the same temperature range (see Fig. I5-12|) . 

5.9 Summary of the main results obtained in YBCO 

The interpretation of the anelastic spectra of YBCO starts from the realization that they must 
be dominated by the hopping with consequent reorientation of the anisotropic elastic dipole of 
a concentration < x < 1 of nonstoichiometric O atoms in the CuOx planes. One can convince 
himself or herself of this fact, by noting that in LSCO a concentration 6 < 0.01 of interstitial 
O atoms dominates the (T) curves (Fig. 14-llb ). even though these hopping defects are 
in a first approximation isotropic (in a tetragonal cell). The barrier for the jumps of O in the 
CuOx planes of YBCO has been determined as ~ 1 eV from several anelastic relaxation and 
other types of experiments, but the inspection of spectra at x ~ 0.1, 0.4 and 0.9 (Fig. l5-4b ) 
demonstrates that the situation is actually more complex. The complexity should arise from 
the fact that when an O atom jumps into or out of a chain, not only causes a reorientation of its 
elastic dipole, but also a change of the electronic energy, possibly by tenths of eV, as discussed 
in Sec. [5X71 

5.9.1 Peak PH2 

The peak that is observed at the highest values of x is obviously due to hopping of O atoms 
in the O-I phase, where the O atoms are ordered into parallel chains. From the Q^^ (T) 
curves presented here it is impossible to extract any quantitative information, due to the loss 
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of oxygen during the measurements in vacuum, and therefore I accept the picture based on low 
frequency anelastic experiments jl24j that oxygen hopping requires the presence of a vacancy 
in the neighboring chain. In Sec. 15.4.71 however, it is pointed out that the existing analyses 
overlook the change in electronic energy that occurs when an O atom leaves or joins a chain; 
this fact might drastically reduce the number of jumps that actually contribute to the dynamic 
compliance (due to the depopulation factor Eq. (|2.15|) 1 to only those that leave the electronic 
energy almost unchanged, like the jumps occurring in the twin walls (Fig. I5-Ifl|) . 

5.9.2 Peak PHI 

Also in a broad range 0.3 < x < 0.8 there is a peak with an activation energy ~ 1 eV, but its 
temperature is definitively lower than that of PH2, and it is labeled here as PHI (Fig. 14-llb 
and 15- 7b .): the fact that PHI and PH2 are actually distinct peaks is reliability confirmed by low 
frequency isothermal measurements |140j . The difference between PHI and PH2 apparently 
resides in the fact that the latter occurs in an environment of full Cu-0 chains, while PHI 
might involve all the other jumps, as depicted in Fig. l5-4b (save the possible condition that 
only jumps of the type of Fig. 15-101 are observable). 

5.9.3 Peak P2 

The unexpected consequence of reducing x below 0.3 is that peak PHI disappears in favor of 
peak P2 at much lower temperature (Fig. 14-llb and I5-<S|) . This is the only relaxation process 
that may be assigned to the hopping of isolated O atoms in the CuO^^ plane; in fact, quenching 
experiments exclude that PHI may contain contributions from isolated O atoms (Fig. l5-7b l. 
Considering that hopping of isolated O atoms must produce a Q^^ (T) well visible down to 
concentrations well below 0.01, it must be concluded that this peak is P2, with an activation 
energy of only 0.11 eV. To my knowledge, this is the lowest barrier for oxygen hopping ever 
reported, comparable to that for hopping of interstitial H in metals. The main difference 
between jumps of isolated O atoms and those causing peaks PHI and PH2 is that the latter 
involve changes in the chemical or electronic state, as discussed in Sec. 15.4.71 

In Sec. I5.4.8l it is also discussed how the unscreened electrostatic repulsion between O atoms 
in the oxygen-poor semiconducting phase should prevent isolated O atoms to form pairs or join 
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existing chain fragments, so explaining why a small but finite concentration of isolated O atoms 
survives at the low temperature where peak P2 is observed (otherwise the free O atoms would 
immediately aggregate into the energetically favorable chains). This also explains the long times 
(years) necessary for oxygen to reach an equilibrium configuration at room temperature. 

5.9.4 Peak P4 

At relatively high oxygen contents three (T) peaks are observed, with activation energies 
of 0.08 — 0.21 eV (Fig. I5-12|) . Of these, peak P4 presents a doping dependence compatible with 
a process occurring in the O-I phase, and it has therefore been associated with short jumps of 
the O atoms between the sites of the Cu-0 chains, which have been shown to be off-centre with 
respect to the chain axis by diffraction experiments jl481 114ll] . 

5.9.5 Peak PI 

The evolution of peak PI with doping, instead, is very peculiar (Figs. 15-171 and I5-1<S|) : in fact, 
it appears only for x > 0.87 and then grows more than linearly with x in a doping region where 
virtually all the physical properties, including the superconducting critical temperature, remain 
almost constant. The only physical quantity that has been reported to display a similar doping 
dependence is the amount of holes in the Pz orbitals of the apical O atoms |115j . Such orbitals 
have negligible overlap with each other, so that the Pz holes cannot form a band but have to pass 
to orbitals in the Cu02 planes or CuOx chains in order to move; this process may well require 
an activation energy like that of PI, 0.08 eV, and therefore the mechanism proposed for PI is 
the reorientation of pairs olpz holes (bipolarons); in fact, the motion of single holes, having the 
same symmetry of the lattice, would not cause any anelastic relaxation. A confirmation that 
pairs of Pz holes are involved is the proportionality of the intensity of PI to the square of the 
concentration of these holes (Fig. I5-19|) . 
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Chapter 6 

Ru-1212 



6.1 Structure 

The ruthenocuprates are a relatively new class of HTS. p^Tj which has attracted much interest 
for the apparent coexistence of ferromagnetism and superconductivity |168j . RuSr2GdCu208 
(Ru-1212) has a cell similar to that of YBCO, with Gd instead of Y, SrO instead of BaO layers, 
and Ru02 instead of CuOa; planes. Figure puts in evidence the CuOs pyramids, found also 
in YBCO, and RuOg octahedra, absent in YBa2Cu306+x where at most CUO4 squares in the 
be planes exist in correspondence with the chains. From the structural analogy with YBCO, 
if O vacancies may be introduced in Ru-1212, they are expected to be in the Ru02 planes. 
The RuOe octahedra are slightly rotated about the c axis in an antiferrodistortive pattern (see 
Sec. 16.2.11 below), and the rotation angles decreases from 14.4° at T = K to 13.9° at room 
temperature |169: ; a transformation to the symmetric unrotated structure is expected at higher 
temperature, but not yet observed. 

Doping in Ru-1212 is supposed to be due to the mixed valence of Ru^^/^"*", which, from X-ray 
absorption spectroscopy jl7Uj and bond valence sums [T7T] . is found to be 40-50% Ru^+ and the 
remaining fraction Ru^^. The mixed valence of Ru is reflected in mixed valence of Cu, with each 
Ru^"*" considered to inject a hole in one of the Cu02 planes above or below therefore providing 
the holes in the usual superconducting Cu02 planes. However, the deduced concentration of 
holes P = ^ [Ru^^] = 0.2 — 0.25 in the Cu02 planes contrasts with the transport properties, 
which are typical of underdoped cuprates, therefore with p ~ 0.1. A possible explanation for 
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Figure 6-1: Structure of Ru-1212 with the RuOg and CuOs polyhedra put in evidence. 

this incongruence is that most of the holes are trapped by defects or by magnetic order and 
do not contribute to conductivity |17()j . Of course, the presence of O vacancies would affect 
the charge balance and doping, and therefore also in Ru-1212, as in the other cuprates, it is of 
general interest to study their formation, dynamics and possible ordering. 

In spite of extensive experimental investigations, nominally identical samples may be su- 
perconducting or nonsuperconducting, and there is no consensus yet, about the influence of 
the oxygen stoichiometry and sample microstructure. Some studies find that annealing at high 
temperature in vacuum or inert atmosphere causes considerable oxygen deficiency [1671 1172j 
while others find no influence at all |173| I169j even up to 800 °C [mi. It has then proposed 
that the prolonged annealings affect the cation ordering |175j . the grain boundaries jl76j . or 
the microstructure |177j . 

The c lattice parameter is almost perfectly 3 times greater than a, so that there is little 
mismatch between domains oriented perpendicularly with each other; as a consequence, there 
is little driving force for the domain growth, and the domain size is very small, unless prolonged 
aging is made. For this reason, it has been proposed that the reason why superconductivity often 
presents character of granularity is the high number boundaries where the Cu02 planes meet at 
an angle of 90°; prolonged anneals therefore would favor superconductivity due to the growth of 
these microdomains and consequent reduction of granularity, rather than stoichiometry changes 
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Although interesting effects are found also in the anelastic spectrum of Ru-1212 below 
room temperature in the present Thesis only the high temperature jl78j results will be 
considered. 



6.2 Oxygen vacancies 

Figure 1^21 presents a series of measurements in vacuum on a same sample at heating/cooling 
rates of 1.5 — 3 K/min, reaching increasing maximum temperatures. The elastic energy loss 
coefficient is shown in the lower part, while the upper part is the change of the Young's 
modulus relative to the initial value at T = K. A dissipation peak starts developing above 
600 K, and it is due to O vacancies (Vq); in fact, the oxygen loss when heating above that 
temperature was confirmed by the increase of the oxygen partial pressure during heating ramps 
in the UHV system. The peak, which is centered at 690 K at 860 Hz and we label as PI, is com- 
pletely developed after heating up to 850 K and remains stable during further heating/cooling 
runs up to 930 K. The (T) and E (T) curves measured concomitantly at 9.8 kHz (not 
shown in the figure) where shifted to higher temperature, indicating that the peak is thermally 
activated, as expected from jumps of Vq- The relatively fast achievement of an equilibrium 
spectrum with an intense thermally activated peak, without all the complex phenomenology 
found in YBCO, indicates that only a small concentration 6 of Vq can be introduced 
in RuSr2GdCu208^5, in accordance with other results in the literature. The oxygen loss 
estimated from the mass loss or gain of a sample subjected to outgassing or oxygenation in the 
UHV system is (5 ~ 0.02 - 0.03. 

Since the final spectrum (curves 3 to 5) is stable, it is meaningful to analyze it to make 
deductions on the dynamics of the O vacancies. In addition, in view of the good quality of the 
data and rather high values of Q^^, instead of Q^^ {uJi,T) I chose to fit the imaginary part 
of the compliance s" referred to the value sq extrapolated to T = K (the low temperature 
measurements are not reported here): 
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Figure 6-2: Evolution of the anelastic spectrum of Ru-1212 measured in vacuum at progressively 
higher temperatures. 

so avoiding any approximation of ^ 1; I also took into account the variation of the 
resonance frequencies tUi (T) by 6% and therefore s' by 12% within the broad temperature range 
spanned by the peak. In fact, formulas like eq. 1)2. 5U() are valid for both s" and s' due to the 
relaxing defect, whereas s' generally contains an important temperature dependent contribution 
from the elastic constants. The resulting s" curves are plotted in Fig. 16-31 for both the 1st 
and 5th vibration modes. There are in fact two peaks: the main one, PI, at 680 K and a 
smaller one, P2, at 550 K at 0.8 kHz. Figure IFw?K presents a fit with two peaks of the Cole-Cole 
type, eq. (|2.45|) . It is apparent that such a model of relaxation is inadequate for reproducing 
the main peak; in fact, i) the intensity of the main peak decreases faster than 1/T; ii) the 
parameters for the relaxation time are tq = 2.7 X 10^1^ s and E/kB = 21800 K; Hi) even with 
these unphysical values of tq and E it is impossible to reproduce the asymmetrical broadening 
of the peak. All these features indicate an anomalous enhancement of both compliance 
and relaxation time on cooling, as expected from interacting O vacancies in the 
Curie-Weiss-like approach of Sec. 12.41 
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Figure 6-3: Imaginary part of the compliance of Ru-1212 after stabilization of the two peaks. 
The continuous lines are fits with Cole-Cole peaks (a) and including the Curie- Weiss corrections 
(b). 

Therefore, the fitting function is (see Eq. ()2.5U|) ') 



s"iu,,T) = As—^ (6.2) 

1 + [UJT) 

A~s = ^^, ^ ^ ToexpCB/fcB^ 

T-Tc l-Tc/T ^ ' 

for PI and again the Cole-Cole formula (|2.45|) for P2. The resulting fit is the dashed curve in 
Fig. l6-3b . with the following parameters: tq = 1 x 10^^'^ s, E = 1.33 eV, Tq = 460 K for PI and 
To = 2 X lO"^"^ s, E = 1.4 eV, a = 0.63 for P2. The fit to PI is definitely better, even without 
any broadening parameter, and with a perfectly reasonable value of tq; this demonstrates that a 
mean-field scheme of interaction between the Vq elastic dipoles is able to reproduce 
the salient properties of PI. 

In order to further improve the agreement with the experiment, I introduced a peak broad- 
ening according to the Cole-Cole expression, obtaining the continuous lines with the parameters 
To = 1.1 X 10"^^ s,E= 1.46 eV, Tq = 470 K, a = 0.82 for PI and tq = 2 x 10"^"^ s, E = 1.1 eV, 
a = 0.7 for P2. This broadening is not excessive and may be justified by the high density of 
boundaries between the three possible orientations of the c axis (see Sec. 16.1(1 and of antiphase 
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boundaries between domains with opposite directions of the rotation angle of the octahedra 
about the c axis; the parameters of PI, including the onset temperature Tq of orientational 
ordering of the elastic dipoles are little affected, but those of P2 assume more physical values. 
The fact that the fit is not perfect may be ascribed to an inadequacy of the description of 
the dipole interaction in the mean-field approximation, but also to a variation of the various 
physical parameters in the broad temperature range from 500 to 900 K. A slow variation in the 
electronic and/or structural parameters may result in a temperature dependent activation en- 
ergy E. For example, the rotation angle of the RuOg octahedra is a linearly decreasing function 
of temperature within the — 300 K range and probably continues decreasing also above 
room temperature, certainly affecting the hopping parameters of the O vacancies. Incidentally, 
one might expect a structural transformation to non-rotated octahedra above room 
temperature, similarly to the case of LSCO and many perovskites, but the anelastic spectra 
do not show any trace of such a transformation up to 920 K. 

Regarding the minor peak P2, no clear indication of Tc > can be found, also due to the 
fact that the peak hardly emerges from the tail of PI, but it must also be related with Vq- 
Two possible origins of P2 are: i) the presence of Vq also in the Cu02 planes, although with a 
much smaller concentration than in the Ru02 ones; ii) Vq trapped by some defect in the Ru02 
planes, like Cu substituting Ru. 

6.2.1 Rotations of the RuOe octahedra 

The fact that the RuOg octahedra are rotated about the c axes (by about 14° at room tem- 
perature) complicates only slightly the treatment of the anelasticity due to oxygen hopping 
within the Ru02 planes. Figure 1^^^ shows the case of non-rotated octahedra, and the x and y 
axes are chosen parallel to the Ru-0 bonds (at 45° with the standard setting of the a and b cell 
parameters). The elastic dipoles A associated with Vq will therefore have major and minor axes 
parallel to x or y; the ellipsoids represent the principal values of A at two O vacancies, compared 
with the unit circles of no strain at regular O sites (the deviation from unity is exaggerated, 
and actually an anisotropic expansion instead of contraction might occur, as shown for YBCO 
in Sec. I5.4.2|) . Figure l^^b includes the rotations of the octahedra, again exaggerated, and it 
can be seen that the strain ellipsoids at the O vacancies remain oriented along x or y, which are 
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still symmetry planes of the structure. In conclusion, it is possible to neglect the rotation 
of the octahedra in analyzing the anelasticity from the hopping of O vacancies. 




Figure 6-4: RuOg octahedra and elastic dipoles (ellipsoids) of two O vacancies neglecting (a) 
and considerin the rotations of the octahedra about the c axis (b). The dashed unit circles 
represent the absence of strain. 

One might ask why collective dynamic rotations of the octahedra of the type 
observed in La2Cu04 are not found in Ru-1212, even though it has planes of octahedra 
able to rotate without constraints from outside the plane. The reason is that soliton-like rotation 
waves cannot develop for rotations about the c axis of octahedra coupled in the ab plane; in 
fact, it is not possible to decouple rotations along different rows as for the LTT pattern of 
La2Cu04 and therefore it is not possible to derive a one-dimensional equation of motion with 
solitonic solutions. This can be readily checked in Fig. I(>4| where it appears that the rotation 
of a single octahedron determines the rotations of all the octahedra in the same plane. 

6.2.2 Elastic dipole of the oxygen vacancy and comparison with YBCO 

In view of the low concentration 6 of Vq in the Ru02_5 planes, it is natural to think in terms 
of relaxation of Vq possessing an elastic dipole A^, instead of jumps of O atoms with elastic 
dipole A^, but the two treatments are otherwise equivalent, with A*-* = — A^. This is apparent, 
for example, also from the symmetry in c and 1 — c in the equations describing the relaxation 
of interacting dipoles in Sec. 12.41 where c = | for RUO2-5 or c = | for CuO^:. From the 
intensity of peak PI and the estimated concentration 5 = 0.02 — 0.03 of Vq one can evaluate 
the anisotropic part IA2 — Ai|of the elastic dipole of Vq- The estimate is very rough, due to 
the polycrystalline and highly porous (~ 50%) nature of the samples. As already done in Sec. 
15.4.21 and shown in detail in Appendix B, neglecting porosity (which is of the order of 50% in 
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the present samples), the relaxation strength for Vq hopping is 

where the Curie- Weiss correction has been included. The Young's modulus of our sample 
deduced from Eq. (|3.2() was E = 40 GPa, probably due to the large porosity. Inserting 
vo = 170 X 10-24 cm^, c = 0.025, Tq = 470 K, A = 2Q~l^ = 0.025 at 680 K one obtains 
AA ~ 0.4; the error in this estimate is large, due to sample porosity and to uncertainties in 
the determination of c, but is about 20 times larger than the value estimated in Sec. 15.4.21 for 
YBCO from the lattice parameters. A possible justification for the large difference is that the 
assumption that = — A^ would be valid in the absence of electronic effects, which instead 
are important. 

6.2.3 Possible roles of the oxygen vacancies in determining the supercon- 
ducting properties 

The samples used in the present investigation had been annealed for 1 week at 1070 °C in 
flowing O2, in order to achieve superconductivity, but no treatments at high O2 pressure were 
done. Still, the superconducting transition was in two steps: a first step at 45 K, indicative 
of intragrain superconductivity, and a second one at 19 K, when the supercurrents flow also 
between different grains. In these conditions, the anelastic spectrum above room temperature 
is flat until the oxygen loss starts above 600 K; looking at the Q^^ (T) curve 1 just below 
600 K it appears that the initial concentration of O vacancies is at least one order of magnitude 
smaller than after outgassing, i.e. RuSr2GdCu208_5 with 6 < 0.002. This means that even 
without high pressure oxygenation the concentration of O vacancies in Ru-1212 is very small, 
and, also after outgassing, it does not exceed 5 ~ 0.02 — 0.03 (see Sec. 16. 2() . These values of 
5 cannot account for the deficit of conducting holes with respect to the concentration of Ru''"'" 
ions, as discussed in Sec. 16.21 and therefore the present study supports the view jl70j that the 
scarce reproducibility of the superconducting properties from sample to sample are due to the 
microstructure rather than to oxygen stoichiometry. 
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Chapter 7 



Conclusions 

The main findings described in the present Thesis from the analysis of a large number of 
anelastic spectra of HTS cuprates can be divided into three types: i) information on the diffusive 
motion of nonstoichiometric oxygen and its ordering; ii) local and collective motion of oxygen 
between off-center positions; in) low frequency fluctuations of the hole stripes pinned by the 
dopants. 

7.0.4 Oxygen diffusion 

The most relevant results on the oxygen mobility have been obtained in YBa2Cu306+a;, where 
stoichiometry and ordering of oxygen play a fundamental role in determining hole doping. It 
has been demonstrated that, although it is generally considered that there is a barrier of ~ 1 eV 
for oxygen hopping in the CuOa; planes, there are actually three distinct types of jumps: in 
the almost fully oxygenated O-I phase, possibly requiring the presence of an O vacancy in the 
otherwise filled Cu-0 chains; jumps involving more sparse Cu-0 chain fragments; jumps of 
isolated O atoms in the semiconducting state over a barrier of only 0.11 eV. The latter process 
implies an extraordinarily high mobility of the ions in the nearly empty CuOx planes, 
but is not easily identifiable with other experimental techniques, since the O atoms have a 
strong tendency to form chains, and their concentration as isolated ions is very small already at 
room temperature. A discussion of the apparent paradoxes posed by such a short hopping time 
(lO"^-*^ s at room temperature) and the long characteristic times for reordering (years) has been 
provided, based on the hypothesis that the electrostatic repulsion between ions is scarcely 
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screened in the semiconducting states, and rises the barrier for oxygen clustering. An analysis of 
the oxygen jumps involving changes in the electronic state, also due to charge transfer between 
Cu-0 chains and Cu02 planes, demonstrates that the situation is more complex than generally 
supposed in other anelastic spectroscopy studies. 

The ruthenocuprate RuSr2GdCu208 is isostructural with YBa2Cu306+a;) except for having 
Ru02-5 planes completely filled with oxygen, instead of the CuOa; chains at most half-filled, 
and it is of interest to ascertain whether the O vacancies are responsible for the difficulty 
in obtaining reproducible superconducting states. The high temperature anelastic spectra of 
RuSr2GdCu208-5 show that their concentration remains low even after long annealings in 
vacuum. It is also shown that the O vacancies jump over a barrier of 1.4 eV, create an anisotropic 
distortion (elastic dipole) much larger than that of oxygen in the CuOa, planes of YBCO, and it 
is possible to describe the high-temperature anelastic spectrum of Ru-1212 taking into account 
the long range elastic interaction among these elastic dipoles. 

In La2-a;Sr2;Cu04 it is possible to introduce excess oxygen in interstitial positions between 
layers of CuOe octahedra. It has been shown that, although in the ideal tetragonal structure 
oxygen jumps between these interstitial positions do not cause reorientation of the elastic dipole, 
this occurs in the slightly distorted orthorhombic structure, and the elastic energy loss peaks 
associated with jumps of isolated and paired interstitial O atoms have been identified. The 
formation of stable pairs confirms the existence of peroxyde species lowering the hole trans- 
fer to the conducting Cu02 planes, but it is also shown that such pairs form at much lower 
concentrations than formerly believed. 

7.0.5 Off-centre positions and tilts of the octahedra 

In addition to the diffusive motion of nonstoichiometric oxygen, the anelastic spectra demon- 
strate that oxygen may perform smaller jumps between off-center positions both in YBCO and 
LSCO. In the first case, an anelastic relaxation process has been identified with the hopping of 
oxygen between off-center positions on either side of the Cu-0 chains (the zig-zag nature of the 
chains is suggested by diffraction experiments showing anomalously high oxygen thermal fac- 
tors); it is also suggested that an ordering transition of the chain O atoms in (anti)ferroelastic 
and (anti)ferroelectric fashion might be responsible for anomalies observed between 120 and 
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170 K on heating. 

In LSCO off-center hopping is connected with the unstable tilts of the oxygen octahedra, 
responsible for the structural transformations; the anelastic spectra presented here demonstrate 
that at temperatures much lower than that of the tetragonal-to-orthorhombic transformation 
and where also the twin motion is frozen, the dynamic modulus (1 kHz) of stoichiometric and 
defect-free La2Cu04 presents still two hardenings of ~ 25% on cooling below 150 K and of 
10% below 10 K; these magnitudes are huge, if one considers that they are not associated to 
any structural transformation. The first is thermally activated and has been identified as tilt 
waves of solitonic type within the planes of octahedra, while the second has been identified 
as faster and local motion of single octahedra or O atoms, with the dynamics determined by 
quantum tunneling. The fact that octahedral tilts are involved has also been confirmed by 
NQR measurements, but a consistent picture of these lattice motions has been possible thanks 
to a complete analysis of the effect of doping through excess oxygen and substitutional Sr. It 
has therefore been possible to clearly distinguish between collective and local motion and to 
determine that the local tunneling is strongly coupled to hole excitations, with as little as 3% 
doping making its dynamics so fast to bring the maximum of the relaxation process below 1 K 
for frequencies of kHz. This finding suggests that at higher doping the anharmonic potential 
for the O atoms might be of the type that has been proposed to enhance the electron-phonon 
coupling to levels necessary for explaining the high-Tc superconductivity |100j . It has also been 
shown that the spectrum of the hole excitations responsible for the tilt-hole coupling is different 
from the known cases of electron-lattice interaction governing atomic tunnel systems in other 
metals or superconductors. 

7.0.6 Hole stripes and antiferromagnetic clusters 

It has been a surprise to discover that the anelastic spectra of La2-xSr^Cu04 contain clear 
signatures of the freezing of the Cu^"'' spins into the cluster spin glass state at Tg (x): a steep 
rise of the acoustic absorption at Tg. This absorption has been interpreted in terms of stress 
induced motion of the walls between the antiferromagnetic clusters with uncorrelated directions 
of the staggered magnetization, and presenting an anisotropic strain thanks to magnetoelastic 
coupling. A neutron scattering study |^ is also confirmed, where, at variance with the prevalent 
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opinion, it is proposed that an electronic phase separation occurs for x < 0.02, with cluster spin 
glass domains coexisting with long range antiferromagnetic order. 

The interest in this finding is enhanced by the fact that, according to the prevalent interpre- 
tation, these domain walls coincide with the hole stripes, whose role in high-Tc superconductivity 
is much discussed, and whose fluctuation dynamics has not been probed by other experiments. 
Carrying on the analogy between the hole stripes interacting with the Sr'^^ dopants and other 
linear defect pinned by impurities, a relaxation process around 80 K at 1 kHz has been identi- 
fied as due to stripes overcoming the collective pinning barrier by thermal activation, with the 
measured barrier in very good agreement with a theoretical estimate j^. The assignment of 
such a relaxation process to the depinning stripe mechanism has been done after the analysis 
of a large number of anelastic spectra with Sr doping ranging from x = 0.008 to 0.20. In this 
manner, it has also been possible to evidence the effect of stripe-lattice coupling through the 
tilts of the octahedra, with clear changes when the prevalent stripe orientation passes from par- 
allel to diagonal with respect to the lattice modulation, and near the doping x = |, where the 
stripe spacing becomes commensurate with the lattice modulation. The locking of the stripes 
with domains of low-temperature tetragonal phase near x = | has also been seen as a strong 
depression of the absorption in the cluster spin glass phase in samples containing Ba instead 
of Sr, due to the reduction of the fraction of mobile stripes. I would like to point out that 
such phenomena are hardly detectable by other spectroscopies, which are dominated by the 
response of single charges and spins (dielectric and magnetic susceptibilities, NMR, /iSR) or 
may probe much shorter characteristic times (light and neutron spectroscopies). The anelastic 
spectroscopy is instead insensitive to individual charge and spin excitations, and probes only 
those that are coupled to strain. 

Even though it is commonly believed that anelastic experiments are too much affected 
by grain boundaries, "defects" and other vagaries, to be able to detect more fundamental 
excitations or processes, they are as reliable as any other spectroscopy, and the anelastic spectra 
are perfectly repeatable after years if the samples are good. I hope I showed that anelastic 
spectroscopy is able to provide reliable and often unique information not only on diffusive, 
reorientational and tunneling atomic motions, including collective modes like tilts of oxygen 
octahedra, but also on elusive charge and spin structures like the hole stripes in high-Tc cuprates. 
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Chapter 8 

Appendices 



8.1 Appendix A - analytical approximation of the order param- 
eter in the Bragg- Williams model of the CuOa; planes 

It is convenient to introduce the order parameter ^, such that 

Ci,2 = |(l±0 (8.1) 

with ^ = and 1 in the completely disordered and ordered states, respectively. Equation (|2.65|) 
becomes 

A solution ^ 7^ is possible when the slope of the logarithm in the left hand at small ^, 
is larger than the slope of the linear term, in the right hand equation; and this determines 
the critical temperature 

kBTc = - I) (8.3) 

Figure 15^ shows the order parameter ^{t), solution of eq. ()8.2|) computed with Mathematica, 
versus the reduced temperature t = T/Tq for some values of the concentration of dipoles c. For 
t > 0.2, the ^ (t) curves become steeper with c increasing up to 0.7, and then acquire again a 
shape similar to the low-c case. Note that for the case of oxygen in the CuOx planes of YBCO, 
there are two oxygen sites per Cu atom, so that in principle < x < 2 and the concentration 
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has to be defined as c = |; however, without doping in the Y or Ba sublattices, it is < x < 1, 
or c < 0.5. Instead, the RuO^ planes of Ru-1212 have c ~ 1. 




0.0 



0.0 0.2 GA^^.^ Q.Q 1.0 

Figure 8-1: Order parameter vs reduced temperature for various concentrations. 

The hmit t ^ or ,^ ^ 1 can be solved analytically by setting ^ = 1 — y and obtaining, to 
first order in y, 



In 



y(l-c)y {2-c)t 



■ exp 



(2-c)t 



(8.4) 



For fitting purposes, it is convenient to have an analytical approximation to (,{c,t). A 
tractable expression, valid for < c < 0.8 and > 10~'^, is the following: 



e(c,t) = (i-rr 



(8.5) 



with n and m depending on concentration as (see also Fig. I8-2|1 



n = 3.153 + 1.785c - 2.681c^ + 9.312c^ - 8.293c'' 



m = 0.51 + 0.02c 



(8.6) 

(8.7) 



For the case c > 0.8, relevant for Ru-1212, a simple expression like (|8.5jl is inadequate, but 
luckily the anelastic relaxation due to oxygen hopping in Ru-1212 at frequencies above 10^ Hz 
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Figure 8-2: Values of the m and n parameters in eq. (|8.5j) for fitting ^ (c, t). 
is well above Tq and there is no need of computing ^ (c, t). 

8.2 Appendix B - relaxation strength in a tetragonal polycrys- 
tal 

The Young's modulus of a tetragonal crystal along the direction n is (see e.g. |12j ) 

E'^ (n) = Sii {nf + 4) + 533n| + (^44 + 2S13) {nj + n^) (8.8) 
+ (566 + 2512) nlnl + 2S'i6nin2 (n^ - 71-2) 

where Siq = for various classes, including 4/mmm, which is appropriate for Ru-1212 |169j and 
tetragonal YBCO; in this case eq. 1)8. 8(1 coincides with eq. (6.2-3) of |H] and eq. (A5) of jl8Uj . 
The O atom or equivalently the O vacancy in the CuO^; or RUO2-5 planes can be considered 
as a (100) orthorhombic defect, which causes relaxation of the elastic constant |180| l8] 

5{Sn-S^,) = '^^{X^-X2f =6S^K (8.9) 
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Adopting the group theoretical comphances of a tetragonal crystal from Eq. (3.5), p. 125 and 
Table 5 of |18Uj . the elastic strain = Sijaj in matrix form may be written as 
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(8.10) 



In terms of these compliances eq. ()8.8() can be written as 

E-' (n) = i (5^'22 + S^^) (nf + 4) + + (^25^.12 ^ 25^) (n? + n^) 

566 +25l2 



In the hypothesis that stress is uniform over the polycrystalline sample (Reuss approximation) 
in order obtain the average ( E~^^ one has to average eq. (|8.11|) over the angles, with 



{E-') = is^'ii + ^S^''' + + Is^^ + Is^^ + ±S^ = (8.12) 

'5 15 15 15 15 15 

= I {2Su + 533) + 4 (2*^44 + See + 45i3 + 2^12) 
5 15 

The relaxation of the compliance due to jumps of the orthorhombic defect comes only from 
6S^^ / and therefore 

5 {E-^) = ^dS""' = ^5 {Sn - S12) (8.13) 
^ ' 15 15 

and the relaxation strength is 

_ 2 MS11-S12) _ 1 cvQ 2 ...4. 
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Porosity strongly affects the elastic moduli jl64j and also the relaxation strength in a non- 
trivial manner, and I simply propose the naive argument that a fraction p of void due to porosity 
should decrease both E and SE^"^ hy 1 — p (neglecting the effects of the shape of the pores on 
E); therefore it should be 

{5E-^) „ „ 5E-'^ 
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List of acronyms and symbols 



Acronyms 

AF = antiferromagnetic 
CSG = Cluster Spin Glass 
DW = Domain Wall 

EXAFS = Extended X-ray Absorption Fine Structure 

HTS = high-temperature superconductor/superconductivity 

HTT = High- Temperature Tetragonal 

LBCO = La2_^Ba^Cu04 

LSCO = La2-xSr^Cu04 

LTO = Low- Temperature Orthorhombic 

LTT = Low-Temperature Tetragonal 

/LtSR = muon spin relaxation 

NCO = Nd2Cu04+5 

NMR = Nuclear Magnetic Relaxation 

NQR = Nuclear Quadrupolar Relaxation 

Ru-1212 = RuSraGdCuaOg 

SG = Spin Glass 

TS = Tunnel System 

YBCO = YBaaCugOe+x 
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Symbols 



c = molar concentration 

Cij = elastic stiffness (matrix notation) 

A = relaxation strength 

AE = energy difference between two states in a relaxation process 

E = Young's modulus, activation energy 

Ep = pinning barrier for the stripes 

£ij = strain tensor (ej in matrix notation) 

J = spectral density (Eq. (|2.36|) ) 

Lc = collective pinning length of the stripes 

\ij = elastic dipole tensor (Aj in matrix notation) 

= oxygen 

01 = interstitial O atom 
p = hole density 

= ^ = ^ = elastic energy loss coefficient 
Sij = elastic compliance (matrix notation) 
aij = stress tensor (cjj in matrix notation) 
T = relaxation time 
Tc = superconducting temperature 
Tq = Curie- Weiss-like ordering temperature 
Td = temperature of the LTO/LTT transition 
Tg = temperature for freezing into the CSG phase 
Tn = Neel temperature for long range AF ordering 
Tt = temperature of the HTT/LTO transition 
Vo = oxygen vacancy 
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